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FOREWORD 

This  document  is  a  comprehensive  collection  of  equations, 
data,  and  supporting  information  for  use  in  processing  aircraft 
flight  data  and  correcting  it  to  standard  conditions.  It  does 
not  represent  a  final  effort;  as  n  w  material  comes  into  being, 
new  sections  and  revisions  to  existing  sections  will  be  issued. 

The  authors  are  indebted  to  Mr.  William  R.  Turley.,  Aerospace 
Engineer,  who  prepared  the  equations  on  the  dynamic  response  of 
vanes  which  appear  in  section  V,  and  to  Mrs.  Barbara  L.  Smith, 
Engineering  Technician.,  for  typing  all  of  the  equations. 

Prtpared  by:  Reviewed  and  appreved  by 
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EVERETT  W.  OUMLAP  - J.  ^IL 
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CtOBERT  M.  WHITE 
Brigadier  General,  USAF 
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4iBSTRA0t 

The  contents  of  this  report  wej^e  prepared  to  give  ^lose 
engaged  in  aircraft,  flight  test  an  xinderstanding;  of  the  analysis’ 
required  to  arrive  at.  standardized  flight;  data.  Toward  that  end, 
considerable  attention  was  given  to  the  derivation  of  equations. 

In  contrast  to  ear  Her  reports/  simplifying  assximptions  were  not 
made;  rather,  efforts. -were  made  to  Keep  the  derivations  of  all 
equations  as  nearly'  exact  as;  possible.  Emphasis  has  been; placed, 
on  climbs  and  levell  accelerations  since  these  tests,  particularly 
for  supersonic  aircra^,  consume  a  large  part,  of  a  test  progr^; 
and  require  calculaticjhs  which  are  much  .more  lengthy  than  -for 
other  tests.  The  information  in  this  document  was  the  basi.s  for 
the  development  of  uniform  digital  computer  programs  which  are 
being  constructed  for  use  in  processing  flight  data  and  .correcting 
it  to  standard  conditions.  These  programs  have  been  given  the 
name  Uniform  Flight  Test  Analysis  System  (UFTAS). 
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INTRODUCTION 

This  document  is  an  outgrowth  of  the  development  of  uniform 
digital  computer  programs  which  have  been  constructed  for  use  in 
processing  aircraft  flight  test  data  and  in  correcting  it  to  a  set 
of  standard  conditions.  These  programs  have  been  given  the  name 
Uniform  Plight  Test  Analysis  Systems  (UFTAS)  and  have  as  their 
basis  the  equations  developed  herein. 

In  contrast  to  earlier  reports  (e.g.,  Flight  Test  Engineering 
Handbook ,  AF  Technical  Report  No.  6273,  by  Russel  M.  Herrington, 
et  al.)  simplifying  assumptions  were  not  made;  rather  efforts  were 
made,  to  keep  the  derivations  of  all  equations  as  nearly  exact  as 
possible.  This  was  done  so  that  the  equations  would  remain 
valid  as  airplane  speeds  and  altitudes  increase  and/or  the  accuracy 
of  instrumentation  systems  is  improved.  In  general,  correction 
*  terms  in  the  final  equations  may  be  deleted  when  it  is  fo\uid  that 

their  magnitudes  are  small  compared  to  the  accuracy  desired  in 
the  end  results. 

Considerable  attention  was  given  to  detailing  the  derivations 
and  to  presenting  information  to  give  the  reader  an  understanding 
of  the  analysis  required  to  arrive  at  standardized  flight  data 
rather  than  incorporating  "cookbook"  procedures.  Although  most 
of  the  topics  covered  are  not  dependent  on  the  type  of  power 
plant,  performance  parameters  and  corrections  to  standard  condi¬ 
tions  have  been  included  for  jet  powered  aircraft  only.  Emphasis 
has  been  placed  on  climbs  and  level  accelerations  since  these 
t  tests  consume  the  major  portion  of  a  test  program,  particularly  for 

y 


supersonic  aircraft.  Further,  calculation  of  tent  parameters 
and  standardization  procedures  are  much  more  lengthy  than  for 
other  tests  (e.g.,  cruise  control,  turning  performance,  etc.) 
and,  hence,  the  difficulty  in  computing  standard  performance  for 
climbs  and  level  accelerations  is  much  greater.  The  basic 
approach  is  quite  simple,  however,  and  is  illustrated  in  figure  1. 

Briefly,  the  method  is  as  follows:  first,  the  performance 
parameter,  excess  thrust,  is  confuted  from  the  aircraft's  measured 
performance  via  the  equations  of  motion.  Next,  the  excess  thrust 
computed  from  test  conditions  is  corrected  to  standard  conditions. 
Finally,  standardized  climb  or  acceleration  parameters  are  com¬ 
puted  from  the  standardized  excess  thrust  by  means  of  equations 
of  motion. 

The  section  titled  Determination  of  Excess  Thrust  discusses 
in  detail  the  diverse  aircraft  flight  parameters  which  may  be  used 
in  these  calculations.  For  many  years  readings  of  airspeed  and 
altitude  from  conventional  airspeed  indicators  and  altimeters  have 
been  a  basis  for  measuring  test  performance.  The  accuracy  of 
excess  thrust  computed  from  such  readings  could  be  described  as, 
at  best,  adequate.  Data  obtained  at  low  altitudes  and  speeds  are 
satisfactory;  however,  as  altitudes  and  especially  speeds  are  in¬ 
creased,  the  use  of  airspeed  indicators  and  altimeters  as  the  sole 
Sources  of  performance  data  becomes  much  inferior  to  sensitive 
accelerometer  installations.  A  precision  of  less  than  +0.01g 
in  measured  flightpath  acceleration  has  been  obtained  during  flight 
tests  at  the  Air  Force  Flight  Test  Center  and  improvements  arv- 
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expected.  Although  not  in  general  use,  procedures  for  calculating 
excess  thrust  from  aircraft  position  measurements  (radar,  Askania 
camera,  etc.)  have  been  developed.  Accuracy  cf  excess  thrust  com¬ 
puted  with  data  from  these  systems  is  superior  to  that  from  the 
airspeed-altitude  method. 


The  most  convenient  parameter  with  which  to  work  in  standar¬ 
dizing  airplane  performance  data  is  excess  thrust.  Corrections, 
itt’  terms  of  excess  thrust,  are  most  easily  derived  from  equations 
of  motion.  If  accelerometers  are  used,  excess  thrn.st  can  be 
computed  from  measured  acceleration,  and  it  can  be  related  easily 
to  rate  of  climb>  turning  performance,  and  other  perfoannance 
parameters . 

In  writing  the  equations  of  motion  to  deduce  test  excess 
thrust >  (or  other  perfopaance  parameters)  the  assumption  of  a  flat, 
non-rotating  earth  with  a  constant  gravity  field  has  historically 
been  made.  (See  for  example,  AF  Technical  Report  6273,  Flight  Test 
Engineering  Handbook.)  To  ensure  that  test  performance  data  are 
not  degraded,  exact  equations  of  motion  have  been  written  which 
account  for  effects  from  the  earth's  rotation  and  for  variation 
in  its  gravity  field.  General  equations  in  terms  of  test  excess 
thrust  are  derived  in  the  section.  Determination  of  Excess  Thrust, 
for  various  flight  parameters.-  In  .preparation  for  the  de.rivatiqn 
of  these  equations,  the  coordinate  systems  employed  are  described 
in  the  following  section.  Coordinate  System?  and  Transformations, 
and  the  properties  of  the  earth  used  in  the  derivations  are  con¬ 
tained  in  the  section.  Geophysical  Properties.  In  the  next  two 
sections,  information  may  be  found  on  standard  atmospheres  and  on 
the  basic  measurements  needed  to  determine  airspeed,  altitude, 

Mach  number,  and  temperature. 

As  in  the  calculation  of  test  excess  thrust,  formulation  of 
excess  thrust  at  standard  conditions  has  been  carried  out  in  a 
more  nearly  exact  fashion  than  is  found  in  AP  Technical  Report 
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make  certain  that  accuracy  inherent  in  the  test  data  is  not  de¬ 
graded  in  the  standardization  process.  A  detailed  derivation 
of  the  equations  for  extrapolating  test  excess  thrust  to  standard 
conditions  appears  in  the  section.  Standardization  of  Excess 
Thrust. 

Once  excess  thrust  for  standard  conditions  has  been  deteimiined 
other  parameters,  (e.g.,  time  to  climb,  standard  weight,  etc.) 
follow  quite  easily  for  either  climbs  or  level  accelerations. 
.Equations  for  computing  these  parameters  appear  in  the  section,. 


SUMMARY 

General  information  about  coordinate  systems  is  presented:; 
this  is  followed  by  descriptions  of  t|ie  various  qoofdinate  sys^ 
terns  used  in  the  calculation  of  excess  ^riist.  Matrices  to 
make  trems formations  from  one  jixis  ^system  tq  another  Sire  de^ 
veloped.  This  is  done  by  first  treStin^g,  .the- genera|?  c^^ 

,  y  /  '  ' 

then  the  specific  cases  which  are;  needed  in  section  V,  Determina~ 
tion  of  Excess  Thrust. 
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SYfffiOLS  USED  IN  THIS  SECTION 


Symbol 

Definition, 

Units 

£ 

aircraft  roll  angle  about  the 
wind  x^raxis  ^airspeed  vector) 

rad 

eg 

ibenter  of  ^gravity 

pet  MAC 

geometric  altitude  of  the  radar 
coordinate  origin  _ 

=^ft 

^uhit  vector s  along  the  S/yi/z 
axes>  respectively 

Tsir- 

axes  transforniation  matrix  for 
the  rotation  through  the  angle 

0  about  the  x-axis ' 

MW. 

multiple  axis  transformation 
matrix. 

XfY/Z 

axes  labels  or  components  of  a 
vector  in^  a  rectangular  Cartesian 
coordinate  system  (appropriate" 
subscripts  denote  the  particular 
axes  system) 

d 

angle  of  attack 

rad 

3, 

sideslip  angle 

rad 

Y 

flightpath  climb  angle  measured 
f  rpiti)  the.  geocentric,  t^rizontal 
plane' 

rad 

r 

geodetic'  latitude  of  the  radar 
coordinate  origin 

rad 

«L 

aircraft  geocentric  latitude 

rad 

geocentric  latitude  of  the  radar 
coordinate  origin 

rad 

A6t 

meridian  tilt  angle  of  the  radar 
axes  horizontal  plane  (=6j^  -6^  ) 

rad 

«r 

=6_  f  the  angle  between  the 

radar  z  -axis  and  the  equatorial 
plane 

rad 

\ 
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6‘jr,0j»6s 

rotation,  angles  for  general 

axes,  transformation 

i^ad 

r 

longitude  of  the  radar 
coordinate  origin 

rad 

difference  between  the 
aircraft 'longitude-  and  the 
longitude  of  l^e  radar  coordi¬ 
nate  origin 

rad 

♦r 

the  angle  bctv;een  the  radar  x  -axis 
and  true  north 

rad 

IT 

the  angle  3 . 14159 . . .  radians 

rad 

a 

fiightpath  heading  angle 
measured  from  true  north 

rad 

[  3 

matrix 

; — 

transpose  matrix 

— 

t 

inverse  matrix 

— 

SiUsscripts 

--  ■  - 

b 

aircraft  b^y  ^es 

.^ef 

geocentric  axes 

— 

g 

north” » gebcgntricaiiy-direc ted 
nbrizpn  axes  (local  geocentric) 

-- 

T 

radar  axes 

— 

,w 

aircraft  wind  axes 

It 2, 3, 4 

general  axes  systems 

I-S 


irn'RODUCTIOK 

In  general,  aircraft  test  data  are  sensed  in  a  translating 
and  rotating  coordinate  system.  In  this  case,  to  coiapute  the 
forces  which  produce  the  Motion,  the  acceleration  data  i;.ust 
be  transforriied  to  an  inertial  coordinate  systeni  to  correct  for 
such  factors  as  centripetal  and  Coriolis  acceleration.  This 
section  describes  the  basic  coordinate  (or  axes)  syster.s  required 
for  the  expressions  and  equations  in  the  other  sections  of  the 
docjunent.  In  addition,  the  necessary  coordinate  transfomsation 
matrices  are  developed. 

COORDIHATb  SYf^TEMS- 

The  axes  systems  used  are  right-handed  rectangular  Cartesian 
coordinate  systems  v;ith  axes  denoted  respectively  by  the  symbols 
x,  y,  and  z  with  appropriate  subscripts  (See  figure  1) .  The 
three  axes  are  mutually  perpendicular,  and  the  direction  of  posi¬ 
tive  rotation  about  each  axis  is  indicated  by  the  curled  fingers 
if  the  axis  is  grasped  by  the  right  hand  with  the  thumb  pointing 
iii  the  positive  linear  direction  of  the  axis.  For  example,  rota¬ 
tion  ai^out  the  z-axis  v/ould  carry  the  x-axis  toward  the  y-axis. 
The  axes  system  is  uniquely  specified  if  the  position  of  the 
origin  and  the  directions  of  tv/o  of  the  three  axes  are  given,. 

The  third  axis  is  then  chosen  to  complete  the  right  hand  system. 

Associated  with  the  axes  system  is  a  set  of  three  unit  vec¬ 
tors,  i,  j,  and  k,  witli  the  same  subscript  as  x,  y,  and  z, 
colinear  respectively  with  x,  y,  and  z,  and  with  the  same  posi¬ 
tive  directions. '  A  vector  can  be  expressed  in  component  form 


as  three  scalars  multiplied  respectively  by  the  unit  vectors. 
GEOCENTRIC  COORDINATE  SYSTEM 

For  the  purpose  of  this  docuiaent  aircraft  trajectories 
will  be  considered  to  be  relatively  near  the  earth's  surface, 
and  the  center  of  the  earth  will  be  considered  an  inertial 
poijit.  The  geocentric  axes  system  (x^,  y^,  z^)  (shown  in 
figure  2)  will  then  be  fixed  in  the  earth  with  its  origin 
at  the  geocenter,  the  z^-axis  pointing  tov/ard  the  south  pole, 
and  the  x^-axis  in  the  equatorial  plane  pointing  toward  the 
earth's  surface  at  a  specified  longitude.  For  development  of 
the  radar  reduction  equations  this  longitude  will  be  chosen  as 
the  longitude  of  the  radar  coordinate  origin.  The  geocentric 
axes  system  is  then  inertial  except  that,  it  rotates  with  the 
earth . 

LOCAI.  GEOCENTRIC  COORDIIWi’L  SYSTEMS 

The  second  coordinate  system  to  be  defined  is  the  north-, 
geocentrically-directed  horizon  axes  system  (x^,  y^,  z^)  or 
local  geocentric-axes  system  for  short  (vSee  figure  3)  . 

Its  origin  can  be  located  at  either  the  surface  of  t}ie  earth 
(geoid)  on  the  radius  line  from  the  ccocenter  to  tlie  aircraft 
CO  or  at  the  aircraft  eg.  The  x  -axis  is  directed  north  and 
the  2  -axis  toward  the  geocenter.  The  y  -axis  is  then  pointing 
east.  These  local  geocentric  axes  arc  oriented  with  respect 
to  the  geocentric  axes  by  two  angles:  6^,  ,  the  aircraft  geo¬ 
centric  latitude,  and  tlie  difference  betv.’oen  the  aircraft 

and  radar  coordinate  oriain  longitudes. 
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AIRCRAFT  WIND  COORDINATE  SYSTEM 

The  third  coordinate  system  to  be  defined  is  the  aircraft 
wind»-axes  system  (x  ,  y  ,  2  )  (shown  in  figure  4) .  Its 

V/  W  W 

origin  is  located  at  the  aircraft  eg  v/ith  the  x^-axis  oriented 

in  the  direction  of  the  aircraft  airspeed  vector  and  the 

s^r-axis  directed  downv/ard  in  the  vertical  plane  of  symii\etry 

of  the  aircraft.  The  y  -axis  is  then  directed  out  of  the 

•'w 

right  side  of  the  aircraft.  The  wind  axes  are  oriented  with 
respect  to  the  local  geocentric  axes  by  three  angles;  the 
flightpath  heading  angle,  a,  measured  from  true  north  to  the 
geocentric  horizontal  projection  of  the  airspeed  vector;  the 
flightpath,  cli^  angle,  y ,i  measured  from  the  geocentric  hori¬ 
zontal  (xy  -plane)  to  the  airspeed  vector;  and  the  roll  angle 
about  the  airspeed  vector,  B.,  The  first  two  of  these  angles 
represents  the  aircraft,  flightpath  orientation  with  respect  to 
the  airmass.  The  angles  are  not  the  same  as  the  angles  of  the 
flightpath  with  respect  to- the  ground  if  the  wind  is  blowing. 
AIRCRAFT  BODY  COORDINATE  SYSTEM 

The  fourth  coordinate  system  to  be  defined  -is  the  aircraft 
body-axes  system  (Xj^,  yj^,  Zj^)  (shown  in  figure  5)  . 

Its  oriain  is  also  located  at  the  aircraft  eg  with  the 

X, -ctxis  directed  forward  along  the  longitudinal  axis  of  the 
b 

aircraft  and  the  Zj^-axis  directed  downward  in  the  vertical 

plane  of  sj^rametry  of  the  body.  The  body  axes  are  oriented 

! 

with  respect  to  the  wind  axes  by  the  sideslip  angle  and  angle 
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of  attack.  The  sideslip  angle,  3,  is  measured  from  the  air¬ 
speed  vector  (x^-axis)  to  its  projection  on  the  Xj^Zj^-plane. 

The  angle  of  attack,  a,  is  measured  from  the  projection  of 
the  Xj^-axis- on  the  x^y^-plane  to  the  Xj^-axis  itself. 

RADAR  COORDINATE  SYSTEM 

The  final  axes  system  to  be  described  is  the  racier  coordi- 

lia'te  system  (x^,  y^,  z^)  (shown  in  the  tv.’o  viev/s  of  figure  6), 

Its  origin  is  located  at  the  radar  site  at  an  altitude  h^.,  a 

longitude  and  a  geodetic  latitude  .  The  x  y  -plane 
r  Uj.  r  r 

is  oriented  v;ith  respect  to  the  geocentric  horizontal  by  a 

tilt  angle  A6_  in  the  meridian  plane.  Usually  the  tilt  angle 
■^r  - - - - 

has  a  negative  value  such  that  the  x^y^-plane  is  parallel  to  the 
geodetic  horizontal.  However,  in  general,  the  radar  plane  may 
not  be  exactly  geodetically  horizontal  due  to  local  anomolies 
in  the  earth's  gravitational  field.  Also,  if  there  is  any  error 
in  the  alinement  of  the  radar  axes  with  respect  to  true  north, 
the  angle  41^  between  the  Xj,-axis  and  the  meridan  plane  will 
be  non-zero.  If  the  radar  plane  is  geodetically  horizontal 

Xj 

r 

can  be  calculated  for  a  given  latitude  using  the  geodetic  to  geo¬ 
centric  latitude  transformation  equation  in  Section  III. 
COORDINATE  TRAi\^SFORMATIOK  MATRICES 

To  illustrate  the  method  by  which  a  transformation  matrix 
is  derived  a  general  three  axis  rotation  will  be  demonstrated, 
and  the  associated  transformation  matrices  will  be  developed 
for  this  rotation.  Succeeding  matrices  will  then  be  obtained 
by  analogy  with  the  general  matrices. 
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GEiNEPJ^  AXES  TRANSFORMATIONS 

From  the  geonietry  of  figure  7  /  which  shows  a  positive 
rotation  about  the  z-axis,  the  following  equations  for  the 
coordinates  (X2>  yj,  za).  in  terms  of  (xi,,  yi ,  zi)  can  be 
derived 

X2  w  >{cqs$^)Xj  +  (si^  +  (O)Zj^  (l)a 

y2  «  (-slnd^)Xj  +  (cos0j)yj  +  {0)z^  (l)b 

Z2  (O)Xj^  +  (0)yj  +  (l)Zj^  (l)c 

This  set  of  equations  is  alternately  expressed  as  the  matrix 
equation  X2  =  MXi,  or  in  component  form 


where  the  subscript  z\  on  M  d^pte?  the  ^is  of  rotation^  and 

0j  denotes  the  angle  of  rotation.  If  the  rotation  were  in  the 

opposite  direction  (through  a  negative  angle  61  which  is 

similar  to  the  rotation  of  (xa,  ya ,  za)  to  (xi,  yi ,  zi)  the 

signs  of  the  two  (sin  6)  terms  in  [M'.(0i)]  would  be  switched, 

z  1 

but  the  other  terms  in  the  matrix  would  be  unchanged.  From 
matrix  theory  the  rotation  matrix  from  (xa,  ya »  Zz)  to 
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(xi/  yw  2x)  i»  knov7h  to  be  the  inverse  of  [II  (6i)]  or 
tM_  (0x)]'’^.  Also,  these  transformation  matrices  can  be  shov/n 
to  be  orthogonal;,  which  by  definition  means  that  the  inverse  is 
>equal  to  the  transpose:  [M„  (9i)]^.  (The  orthog- 

onality  property  further  means  that  the  cofactors  of  each  matrix 
element  must  equal  the  element.  This  property  provides  a  con¬ 
venient  method  of  checking  any  matrix  for  errors  among  its 
eleKients.)  If  the  new  axes  (xj,  yi,  Zz)  are  now  rotated  through 
an  angle  02  about  tlie  ya-axis  to  forru  the  axes  systera  (xs,  ys,  zz)  , 
then  y3=yji,  and  z  v;6uld  replace  x  in  figure  7,  and  x  would 
replace  y;  so,  by  inspection  we  write  the  second  transformation 
matrix 


Similarly,  if  a  transformation  is  raade  to  (xi,,  y4,  z^)  by 
rotation  through  G3  about  the  X3-axis,  then  y  replaces  x  and 
z.  replaces  y  in  figure  7,  so  that  again  the  transformation 
written  by  inspection  is 


The  total  transformation  from  (xj,  yi,  zi)  to  (X4,  y^,  Zk) 
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is  then 


1 - 

X 

1 _ 

r  “1 

'"1 

^4 

-  [“X3<«3]l[«'y2‘«23[“zi‘«l3 

-  02'^3] 

N 

_ 1 

To  make  the  reverse  transformation  from  (xu,  y^f  Zk)  to 
(xi,  yi,  zx)  the  transpose  (inverse)  of  M(6i/  62/  63)  is  required 
The  transpose  is  obtained  by  use  of  the  reversal  law  for  the 
transpose  of  a  product  of  matrices 


(6) 


TRANSFORMATIONS  FROM  GEOCENTRIC  TO  LOCAL  GEOCENTRIC 

By  analysis  of  figure  3  it  can  be  determined  that  two  rota¬ 
tions  are  required  to  transform  the  geocentric  axes  (x^,  y^,  z^) 
to  the  local  geocentric  axes  (x^,  y^,  z^) .  First,  a  rotation 
about  the  z  -axis  throuah  the  angle  (tt+AAt)  yields  the  trans- 

formation  matrix  which  by  analogy  v/ith  (M  (Oi)]  is 

^  1 


( 77  +  AA 


COS(;r  +  AAj^) 

-sin(;r  +  AAj^) 
0 


sin(?r  +  AA^)  0 
COS(  n- -i  AA|^)  0 

0  1 


(7) 


This  is  simplified  by  trigonometric  identities  to  the  following 
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-cosAXj^  -slnAXj_  0 
sinAA|_  -cosAA|_  0 

0  0  1 


(8) 


Second,  a  rotation  about  the  y^-axis  thro?  gh  the  angle  {tt/2-6j^) 
yields  the  transformation  matrix  v/hich  by  analogy  v/ith  [I-Iy^CSa)! 
is 


cos(“-  S,  ) 
2  •- 

0 


0  -sin(f-S,  ) 
2  L 

1  0 

0  cos(4-S,  ) 

2  L 


which  simplifies  to 


slnS^^  0  -cosS|_ 


(9) 


(10) 


Expanding  the  matrix  product  we  obtain  the  matrix  of  the  trans¬ 
formation  from  geocentric  axes  to  local  geocentric  axes 


-cosAAj^sinSj^ 

sinAAj_ 

-cos  AA|_cosS^^ 


-sin  AAj^sinS^^  -cosSj^ 

-cosAAj^  0 

-sin  AA|_cosSj_  sinS^^ 


(11) 


The  inverse  of  this  transformation  matrix  for  the  reverse  trans- 
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formation  from  local  geocentric,  to  geocentric  coordinates,  is 


-cos  AXj^sInS^ 

s|nAX|_ 

rCOSAXj^COSS^ 

^!(5L/AXLj]  ■« 

-sin  AXj^slnSj^ 

-cosAX|_ 

-slnAXj_cosS|^ 

(12) 

-cosSj^ 

0 

TRANSFORMATIONS  FROM  LOCAL  GEOCENTRIC  TO  WIND  COORDINATES 

The  transformation  from  the  local  geocentric  coordinates  to 
the  aircraft  wind  coordinates  consists  of,  rotations  through  the 
flightpath  heading  angle,  o,  the  climb  angle,  y,  and  the  r611 
angle,  B,  in  the  same  sequence  and  about  the  same  axes  as  the 
rotations  illustrated,  through  the  angles  ,6i,  62,  and  63  in  the 
general  three  axes  transformation.  Consequently >  the  resulting 
transfoinnation  matrix  is. 


1 

0 . 

0  : 

cosy 

0 

-rsihy 

cos  {7 

Slha 

6 

- 

0 

cosB- 

sinB 

0 

j 

1 

0 

-sina 

COS  (7 

0 

r 

0 

-sin  B 

cos  6 

l^iny 

0 

cosy 

0 

0 

1 

(13) 


For  one  applies Sion  in  the  section  describing  the  radar  reduc¬ 
tion  equations  only  the  portion  involving  the  first  ty/o  trans¬ 
formations  through  0  and  y  to  an  unbanked  wind  axes  system  are 
required.  Expanding  this  product  first  we  obtain 


r. 


cosy,  cos  <7 


■‘.sina 

sinycosor 


cosy  sin  (7 
cos  a 
slnyslna 


-siny 

0 

cosy 


(14) 


l•*22 


Next,  forming  tte  product  to  rotate,  to  the  banked  wind  axes  we 

obtain  the  total  transfpjmi  matrix 

'  '  '  \ 


* 

- 

cosy  cos a 

cosy  sina 

-siny 

i 

[M(ff,y,B^  » 

-cos  Bsina 
+  sin  B  siny  cos(7 

cos B cos a 
+  sin  B  siny  sina 

sin  B  cosy 

1 

sinBsitkr 

+  cos  B  siny  cos o 

-sinBcosa 
+  cos  B  siny  sina 

cos B cosy 

(15) 


TRANSFOiy^lATIONS  FROM,. BODY  COORDINATES  TO  WIND  COORDINATES 


The  transformation  from  body  coordinates  to  wind  coordinates 
requires  tv;o  rotations  in  the  negative  directions  first  through 
a  and  then  6  about  the  and  z^- axes, respectively.  The  trans¬ 
formation  matrix  is  then 


cos/3  -sin^  0 

cos  a  0 

sina 

V  [M(a,/3^- 

sin/3  cos/3  0 

_  0  0  1 

0  1 

-sina  0 

0 

COSa 

(16) 


When  the  sideslip  is  neglected-  the  total  matrix  reduces  to  the 
second  one  involving  g  only.  When  3  is  included  the  total 
matrix  is 


cos /3  cos  a  -sln/3 
slnjScosa  cosjS 

-sin  a  C 


cos/3  sin  a 
sin/3  sina 
cos  a 


(17) 


I 

i 
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TRANSFORMATIONS  PROM  RADAR  TO  GEOCENTRIC  COOI^INATES 

To  transform  the  radar  measurements  t6  geocentric  cQordinates, 

transformations  through  the  two  orientation  angles  (Tr-<1)^)  about 

the  radar  2  -axis  and  {Tr/2-6^)  (where  6  =6-  -A6-  )  about  the 
r  r  r  Lj.  Lj. 

Yg-axis  are  required.  The  resulting  matrix  is  then 


|m  ( Y~  ■■ 

cos(|-S^)  0  -sMj-Sfl 

016 

cos(n-<f>f)  sinCff-^,,)  0 

-'sin(  n  -  <f)j)  cos(  17  -  ^)  0 

sln(-|-S^)  0  cosij-Sf) 

:  0  0  1 

(18) 


Performing  appropriate  trigonometric  substitutions 


sinSj.  0 

-cos  5, 

-cos<)5|. 

sin^f 

0 

0  1 

0 

-sin^f 

-COS^j 

,0 

cosSf  0 

;sin5f 

0 

0 

1 

Multiplying  the. 

two  matrices  together 

r 

' 

-cos^jislnSf 

sin^fSinSf 

rCOSSf 

-COS^fCOSSf 

-cos^,. 

sin^fCos$p 

0 

sinSf 

(19) 


(20) 
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SECTION  II 


GEOPHYSICAL  PROPERTIES 


II-l 


SUMMARY 


The  usual  assumptions  that  have  been  made  in  the  past 
(a  flat,  non-rotating  earth  with  a  constant  gravity  field) 
lead  to  significant  errors  when  performance  data  are  acquired 
with  inertial  navigation  systems  or  with  accelerometers  (such 
as  two-axis  flightpath  accelerometer  systems) .  Accelerations 
brought  about  by  the  earth's  rotation,  for  example,  can  be 
readily  sensed  by  current  installations  in  test  aircraft. 
Further,  the  magnitude  of  the  errors  caused  by  these  assump¬ 
tions  becomes  larger  as  speeds  and  altitudes  are  increased. 

The  more  nearly  exact  equations  derived  in  section  V,  Determi¬ 
nation  of  Excess  Thrust,  make  use  of  the  model  of  the  shape  of 
the  earth,  its  gravitational  field,  etc.,  set  down  in  this 
section . 


TABLE  OF  CONTENTS 


SYMBOLS  USED  IN  THIS  SECTION _  4 

INTRODUCTION -  6 

DIMENSIONAL  PROPERTIES _ 7 

Equatorial  and  Polar  Radii _  8  . 

Earth  Radius _ 9 

GEODETIC  AND  GEOCENTRIC  LATITUDES  AND  ALTITUDES - 10 

GRAVITATIONAL  FIELD  _ 13 

GEOMETRIC  AND  GEOPOTENTIAL  ALTITUDES _ 17 

GEOPHYSICAL  CONSTANTS _ 21 

REFERENCES  _ 23 


li~8 


Symbol 

f 


SYMBOLS  USED  IN  THIS  SECTION 


Definition 
earth's  flattening 


g,  local  effective  acceleration 

due  to  gravity 


gLv  »9Ly  geocentric  components  of  local 

*9  ^9  "  acceleration  due  to  gravity 

g  reference  acceleration  of 

gravity 


*SL 


sea  level  acceleration  of 
gravity 


H 


geocentric  components  of 
acceleration  due  to  gravita¬ 
tional  attraction  alone 

geopotential  altitude  in  geo- 
potential  units 


h 


geometric  (tapeline)  altitude 


hi 


geocentric  altitude  as  shown 
in  figure  2 


h2 


geodetic  altitude  as  shown  in 
figure  2 


ha  geometric  altitude  as  shown 

in  figure  2 

JzfJsfJ'*  coefficients  of  the  zonal  har¬ 

monics  of  the  earth's  gravita¬ 
tional  potential 


r 

r 


R. 


local  radius  of  the  earth 
polar  radius  of  the  earth 
equatorial  radius  of  the  earth 
effective  earth  radius 
time 


Units 

dimensionless 


ft  per  sec^ 
ft  per  sec? 
ft  per  sec^ 
ft  per  sec* 
ft  per  sec* 

ft 

f,t 

ft 

ft 

ft 

dimensionless 


ft 

ft 

ft 

ft 

sec 
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Symbol 

Definition 

Units 

XrY 

Cartesian  coordinates 

: - 

geocentric  position  coordi¬ 
nates 

ft 

geodetic  latitude 

rad 

geocentric  latitude 

rad 

X 

longitude 

rad 

product  of  universal  gravi¬ 
tational  constant  and  the 
mass  of  the  earth 

ft®  per  sec® 

$ 

gravitational  potential  func¬ 
tion 

f t^  per  sec® 

angular  rotation  rate  of  the 
earth 

rad  per  sec 
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INTRODUCTION.  ^5  -  " 

In  €his:  section  equations  for  the  earth's  geophysical 
properties  are  presented'.  The  geophysical  parameters 
which  are  necessary  for  atmospheric  trajectory  caicula- 
tipris  are  those  which  describe  the  earth ' s  mean  sea  level 
surface,,  the  aircraft's  position  relative  to  that  surface, 
the  gravitational  attraction  between  the  earth  and  the 

•  '  '  i>  -  .  « 

aircraft'^  and  the  earth's  rotation  rate,,  which  provides 
centrifugal  relief  from  the  force  of  gravity.  An  attempt, 
has  been  made  to  present  equations. .which  have  precision 
comparable  td>  the  general  aircraft  equations  of  motion 
?presented  in  the  section,  Detennination  of  Excess  Thrust. 

Various  approximations  to  the  earth's  shape'  and  gravi¬ 
tational  field  have  beerii  made  to  simplify  the  equations  of 
motion.  However^  with  the  introduction  of  high  perfomance 
aircraft  into  the  Air  Force;  inventory, improved  equations 
have  become  desirable.  High  accuracy  accelerometers  and 
precise  ihertiai,  fnayigation  systems  can  sense  accelerations 
such  as  t^ose  induoe<l‘  by  the  earth's  rotation,,  and  they  can 
sense  forces  caused;  by  variation!  in  gravity  .with  altitude 
and  latitude.  These  navigation,  systems  must  also  account 
for  the^  oblate  shape  of  t^©'  earth  in  order  to  provide  accu¬ 
rate  position  data.  The  old  assumptions  of  a  flat  earth 
and  constant  gravity  are  no  longer  adequate  in  most  cases. 
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Since  the  advent  of  artificial  earth  satellites 
a  number  of  refinements  in  the  measurements  of  the  earth's 
geophysical  properties  have  been  made.  These  precise  equa¬ 
tions  with  slight  simplifications  have  been  used  in  this 
section.  An  attempt  has  been  made  to  state  the  important 
assumptions  and  describe  the  simplifications  to  allow  for 
future  analysis  of  the  adequacy  of  the  equations. 

DIMENSIONAL  PROPERTIES 

The  strength  of  the  earth's  crust  is  low  enough  when  com¬ 
pared  to  total  earth  mass  and  rotational  energy  that  its  sur¬ 
face  shape  has  been  forced  to  assume  the  approximate  form 
of  a  ball  of  fluid  in  hydrodynamic  equilibrixam.  Such  a 
ball  of  fluid  in  space  (in  the  absence  of  surrounding 
bodies)  would  assume  a  spherical  shape  under  the  action  of 
internal  gravity;  however,  when  given  a  specific  rotational 
rate  about  a  "polar"  axis  it  would  develop  a  bulge  about 
its  equator  and  assume  the  shape  of  an  ellipsoid  of  rev¬ 
olution. 

In  actual  fact  the  earth  does  not  have  a  circular  equa¬ 
tor  and  its  density  is  not  uniform.  However,  for  reference 

purposes  an  ellipsoidal  surface  is  defined  which  approxi¬ 
mates  mean  sea  level.  Mean  sea  level  is  also  approximately 
the  geoid  surface  of  dynamical  balance  between  gravitational 
force  and  the  inertia  of  the  rotating  mass  (see  figure  1) . 


Figure  1  The  Reference  Ellipsoid 


EQUATORIAL  AND  POLAR  RADII 

The  dimensions  of  the  reference  ellipsoid  are  generally 
given  by  stating  an  equatorial  radius,  r^,  and  a  flattening, 
f,  defined  to  be  the  quotient  formed  by  dividing  the  equa¬ 
torial  radius  into  the  difference  between  the  equatorial 
and  polar  radii: 


This  equation  can  be  solved  to  yield  an  expression  for  the 
polar  radius  in  terms  of  the  equatorial  radius  and  earth 
flattening; 


EARTH  RADIUS 


Substitution  of  these  radii  into  the  polar  coordinate 
equation  for  an  ellipse  yields  a  general  equation  for  the 
radius  of  the  ellipsoidal  surface  as  a  function  of  the  geo¬ 
centric  latitude,  6^,  which  is  the  angle  between  the  radius 
vector  and  the  equatorial  plane: 


( 1  -  f)^  cos^Sj^  + 


One  alternate  form  of  this  equation  is  obtained  by  use  of 
tha  trigonometric  half-angle  formulas  and  by  dividing  by 
(l-f)"5 

9r  2 

^2.  _ _ tl2 _ 

[  1  +  (-X>.)2]  +  [  1  -  (-i_)2]cos25T 
l-f  1-f  L 

Another  computationally  convenient  form  of  this  equation 
can  be  obtained  by  expansion  using  a  Maclaurin  series  and 
appropriate  trigonometric  substitutions.  For  example,  for 
1/r  =  298.30  the  following  expansion  yields  the  same  radius 
as  equation  (4)  to  the  nearest  foot  or  better; 


r  -  ro(.99832172  +  .00167616  cos  2  +  .00000211  cos  4  5^)  (5) 
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GEODETIC  AND  GEOCENTRIC  LATITUDES  .AND  ALTITUDES 

Having  defined  the  reference  ellipsoidal  surface  of  the 
earth /  we  can  turn  to  the  problem  of  measuring  an  air¬ 
craft's  altitude  with  respect  to  that  surface ^  Pour  alti¬ 
tudes  will  be  defined.  An  exaggerate' d  illustration; of 
these  altitudes  is  shown  in.  figure  2.  The  first  is  the 

geocentric .altitude  which  is  merely  the  distance  measured 

% 

along  an  extension  of  the  radius  vector  from  the  ellipsoid 
surface  up  to  the  aircraft.  This  distance  wken  added  to 
the  earth  radius  at  that  point  gives  the  straight  line  disr 
tarice  between  the  earth's  center  and  ;:he  aircraft i  The,  geo¬ 
centric  altitude  is  required  in  the  general  equations  which 
describ  ^  acce? erations  in  the  geocentric  reference  system. 


Figure  2  Altitudes  Related  to  the  Geoid 
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The  second  altitude  to  be  defined  is  the  geodetic  alti¬ 
tude.  This  is  the  straight  line  distance  from  the  aircraft 
tb  the  nearest  point  .on  the  geoid  (sea  level)  surface.  The 
angle  between  this  line  and  the  equatorial  plane  is  defined 
to  be  the  geodetic  latitude,  The  geodetic  altitude  is 

the  true  altitude  or  straight  line  height  of  a  point  above 
the  earth's  surface.  By  observation  of  figure  2  one  can 
see  that  the  differences  between  the  geocentric  and  geodetic 
latitudes  and  altitudes  vary  and  are  a  maximxp  near  a  lati¬ 
tude  of  45  degrees.  These  differences  also  vary  with  alti¬ 
tude  . 

A  precise  equation  for  transforming  between  geodetic  and 
geocentric  latitudes  is  presented  in  reference  1,  pp.  96  ff. 
However,  for  aircraft  trajectories  which  generally  lie  below 
150,000  feet,  the  precise  equations  can  be  approximated  with 
little  error  by  the  sea  level  expression. 

This  relation  betv;een  the  two  latitudes  at  sea  level  is 
easily  derived.  The  slope  of  the  ellipsoid  surface  is 
obtained  by  differentiation  of  its  rectangular  coordinate 
equation : 


(6) 


Then  since  a  geodetic  latitude  line  is  by  definition  every 
whe-  normal  to  the  surface,  its  slope  is  the  negative 
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reciprocal  of  the  surface  slope  : 


tanS^ 


( -52.)^  tan  5f 
■■p 


(7) 


Substituting  for  the  radii  ratio  from  equation  (2)  and 
rearranging  we  obtain  the  simple  exact  relation  between 
the  latitudes  of  a  point  at  sea  level : 

tanSj^  -  (1  — O^tanSp  (8) 

The  maximum  difference  between  the  two  latitudes  is  about 
eleven  and  one-half  minutes  of  arc  or  eleven  and  one-half 
nautical  miles  error  in  position  on  the  earth's  surface  if 
the  incorrect  latitude  were  to  be  used. 

The  relation  of  equation  (8)  can  be  approximated  by  the 
first  term  of  a  series  (reference  2,  p.  485) : 


.  ~(. 19323889) sin 25p  (9) 

Equations  could  also  be  presented  relating  the  geo¬ 
centric  and  geodetic  altitudes,  but  little  error  results 
from  assuming  that  they  are  equal.  The  difference  is  one 
foot  or  less  belov/  150,000  feet  altitude  and  less  than  60 
feet  at  an  altitude  of  about  one  earth  radius  (reference 
1,  p.  102) .  The  remaining  two  altitudes  as  shown  in  figure 
2  are  described  later. 


GRAVITATIONAL  FIELD 


By  definition  the  gravitational  potential  function  of 
the  earth,  is  a  scalar  quantity  such  that  the  partial 
derivative  of  $  in  a  given  direction  yields  the  accelera¬ 
tion  due  to  gravity  (reference  3,  p.  173) .  For  example 
in  the  x  direction  by  definition  (e) 


ii.  B 

"  dt2 


(10) 


If  the  earth  were  a  homogeneous  sphere,  it  would  have  a 
central  force  field  with  a  gravitational  potential  function 
given  by 

^  "  m  (11) 

where  (r+h)  is  the  distance  from  the  geocenter  to  a  given 
point.  Differentiation  of  this  function  in  the  (r+h) 
direction  yields 

5(r  +  h)  d(r  +  h)  r  +  h  .  ,,2 

U  +  n; 

or 


d^(r  +  h)  ^ 
dt^  "  (r  +  h)^ 


(12) 


which  is  the  inverse  square  law  gravity  equation. 
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since  the  earth  is  not  a  homogeneous  sphere  nor  even 

an  ellipsoid  of  revolution,,  its  true  potential  function 
must  be  obtained  by  an  integration  over  its  entire  vol¬ 
ume  or  alternately  as  a  series  solution  to  the  differen¬ 
tial  equation  of  Laplace  (reference  4,  p.  141) .  Neglect¬ 
ing  the  earth's  asphericity  with  longitude  (assuming  the 
shape  is  an  ellipsoid  of  revolution)  and  only  presenting 
the  first  four  terms,  the  series  for  the  earth's  potential 
is  (reference  5,  p.  2) 


4) 


(r  +  h)  2 


(•^)^(l-3sin^5L^ 


+  ~(“~-)^l3  -  S  sin^Sj  )  sin  St 
2  r  +  h  ^ 

Jo^)4(3  _  aosin^s  +  35sin%  )] 
8  r  +  h  ^  *-■ 


The  components  of  gravitational  acceleration  in  the 
local  geocentric  coordinate  system  can  be  obtained  by  par¬ 
tial  differentiation  of  the  potential  function  in  the  three 

coordinate  directions  x  ,  y  .  and  z  .  To  perform  this 

y  y  9 

differentiation  the  following  relations  are  required: 


cr  .  .  JL.J±.  (14a) 

SXg  r  +  h 


1  21 
(r  +  h)sin5j^ 


(14b) 


ll-M 


d<if 

d{r  +  h) 


(14c) 


Performing  the  differentiation  indicated  by  equation 
(14a)  we  obtain  the  north  component  of  gravitational  accel¬ 


eration  : 


•  gx  -  -(-%)('-^)'^[3J2sin5j 

6Xg  p^.},  2  L 

-  J.J3(J^.)(l-5sin25j  ) 

2  ^  r  +  h  L. 

“  f  (  3  -  7  sin^S^ )  sin  ^  J cosSl  (13) 

Since  the  variation  of  the  potential  function  with  longi 
tude  has  been  neglected  the  eastx^/ard  component  of  the  gravi¬ 
tational  acceleration  is  zero ; 


gv_  ■  0 


(15b) 


Finally,  evaluating  equation  (14c) ,  we  obtain  the  geo¬ 
centric  vertical  component  of  the  gravitational  acceleration 


(-^)(-^)^[  1  +  1  J2(-^)(l  -  3sin2s  ) 

r2r+h  ^ 

o 


+  2J3(-~-j-)^(3  -  5sin^5j^)sinSj^ 


r„  ,4 


—  2.  J.(— 9.— )  (3-30sin  5r  +  35 sin  Sr  )] 
ft  ^  r  +  h  L.  L 


(15a: 


(15c 


These  components  yield  the  acceleration  due  to  gravita¬ 
tional  attraction  alone.  To  determine  the  apparent  acceler¬ 
ation  of  a  particle  at  a  fixed  point  with  respect  to  the 


earth's  surface  the  so  caileia  '" 


which" 


is  contributed  to  the  particlfe  by  thCjA^ar^-'s' 
must  be  included.  The  components;  of  hh^  ".^1^ 
relief"  actually  arise  as  a  result  of  the  ,p,articl;e 
motion  with  respect  to  inertial  spao.e,r  and  the>de'rivation 
of  these  components  appears  in  the  section^  Determination 
of  Excess  Thrust.  For  the  present  purpose,  the  "cehtrifur 
gal  relief"  components,  will  be  added  to  the  gravitational 
attraction  acceleration  components  without  further  cpniraent. 
The  resulting  centrifugally  relieved  geocentric  components 
of  the  acceleration  are  then 


^(-^)(-52--)'^[3J2sinSj  -  1  J3(  JL.)(1  _  5sin25.  ) 

2-  r  +  h  ^  'b  2  ^  f  +  h  b 


5  T  (  *0  _  7  \ 


”  “JaI— (3  -  7sih  Sj  )sinSr  ]cos5r 
2  ^  r  +  h  b  L  L 


->  0)^  ( r  +  h )  cos  sin 


(16a) 


(16b) 


(i^)(JL-)^[l  +  3  j  (Jh  )2(i.33i„^g  ) 
2r+h  2r+h  ^ 


+  2J3(-^^)^(3 -»  5sin^Sj^)sin5£^ 

)'^(3-30sin^gr  +35sin'^Sy  )] 
8‘*r  +  h  '  b  L 

-  (a*2(  r  + h)cos^5f 


(160 


I 


^  The  magnitude:  of  ^e  resultant  accel^fatifonid^  t© 
gravity  at  a  point  ^ove  the  earthvls  thexi: 


GEOMETRIC  AND  GEOPOTENTIAL  ALTITUDES 

The  total  gravity  cpinponents  of  equation  (16)  have  been 
obtained  by  suininiug  ^e  final  acceleration  contributions  s 
from  gravitational  attraction  and  "centrifugal  relief."  A 
different  approach  is  taken  in.  reference  6,  p.  ‘5.  The  total 
potential,  or  geopotential  at  a  point  is  defined  as^  the.  sum 
of  the  gravitational  potential  and  the  fictitious  "centrif¬ 
ugal  force"  potential.  The  lines  of  gravity  force  which 
result  from  this  potential  are  by  def ini tioh  evei^where  per¬ 
pendicular  to  the  surfaces,  of  constant  pctential,  and  since 
the  potential  surfaces  are  ellisoids  the  lines>  of  force  are 
pwved  .(reference  6.,.  figure  I,.2v4(a)i)...  The.  curvature,  ijs, 
such  that  latitude  increases  with  altitude  along  the  line 
of-  forces 

Geometric  altitude  is  defined  to  be  the  curved  distance 
measured  along  a  line  of  force  from  the  zero  potential  level 
(approximately  sea  level)  to  the  given  altitude  point. 

It  is  convenient  to  define  another  altitude  parameter 
which  is  especially  useful  in  the  model  atmosphere  equations 
of  the  section.  Atmospheric  Environment.  This  parameter  is 
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physically  equivalent  to  the  geopqtehtiai  but  has  the^ units 
of  length  and  is  called  the  geppbtential  altitude.  It  is 
defined  by  the  equation 

h  -  /  ikdh  (18) 

O'  gr 

where  g^^  is  the  standard  value  of:  the  reference'  sea  level 
acceleration  due  to  gravity,  in  accordance  with  its  defin¬ 
ition  the  geopotential  altitude  is  also  measured  along  the 
curved  lines  of  force;  however,  the  physical  or  geometric 
length  of  a  geopotential  foot  is  not  cpnstantii  The  length 

i 

is  a  Junction  of  the  local  gravity  and  consequently  increases 
with*  altitude  because  of  the  reduction  in  gravity. 

The  most  precise  technique  for  calculating  geopotential 
altitude  requires  a  simultaneous  numerical  Integration  of 
the  differential  equations  relating  the  acceleration  of 
gravity -and  Jhe  latitude  and  radius  of  each^  point  along 
the  curved  line  of  force  to  the  geometric  altitude  of  that 
ppint  (reference  6>  pp.  6.,  7)  . 

The  curvature  of  the  line  of  force  is  exaggerated  in 
reference  6,  emd  for  the  range  of  altitudes  considered  in 
this  document  the  geometric  and  geopotential  altitudes  can 
be  considered  to  be  measured,  along  straight  lines  coinciding 
with  the  geodetic  altitude  lines.  The  error  associated  with 
this  assumption  is  negligible.  Consequently,  the  three 
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can 


different  altitudes,  geocentric,  geodetic  and  geometric, 
all  be  considered  equal,  and  we  can  speak  of  the  geometric 
altitude  as  being  representative  of  all  three. 


By  an  appropriate  assumption  an  equation  can  be  derived 
which  allows  direct  calculation  of  geopotential  altitude 
from  geometric  altitude  and  vice  versa  without  the  necessity 
of  the  numerical  integration  from  sea  level  up  to  the  given 
altitude  (reference  7,  pp.  217,  218  and  488).  The  assump¬ 
tion  which  is  necessary  is  that  the  centrifugal  relief  is 
not  a  function  of  altitude.  By  analysis  of  equation  (16) 
it  can  be  detemnined  that  this  assumption  amounts  to  less 
than  one  percent  error  in  magnitude  of  the  centrifugal  relief 
term  for  altitudes  below  200,000  feet.  Since  the  term  is  a 
small  fraction  of  the  total  gravitational  acceleration  the 
assumption  introduces  negligible  error. 

To  develop  the  conversion  equation  an  inverse  square  law 
gravity  field  is  assuirted,  but  the  earth's  effective  radius, 

R  ,  and«sea  level  acceleration  due  to  gravity  are  expressed 
as  a  function  of  latitude  to  make  the  resulting  equation  fit 
the  actual  gravity  field  more  closely.  For  an  inverse  square 
gravity  field  the  acceleration  due  to  gravity  is 

8L  ■  (19) 
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Substituting  this  gravity  expression  into  equation  (18)  and 
integrating,  we  obtain  the  conversion  equation. 


H 


^SL  ^  'j 

Sr  Rg  +  h 


(20) 


Equation  (20)  can  easily  be  solved  to  obtain  the  inverse 
equation  for  h  as  a  function  of  H. 

To  obtain  the  expression  for  we  first  differentiate 
equation  (19) . 


dh 


(Rg  +  h)^ 


Evaluated  at  sea  level  (h  *=  0) ,  this  becoires 


IfLl  - 

'h-0 


and  solving  for 


R 


e 


5h  h-0 


(21) 


11-20 


IQ 


Equations  for  the  sea  level  gravity  and  the  partial  deriva¬ 
tive  at  sea  level  are  presented  in  reference  (7).  Similar 
equations  can  also  be  derived  by  appropriate  substitution 
and  differentiation  of  equation  (16) . 

The  value  of  geopotential  altitude  at  each  geometric 
altitude  as  given  by  the  more  nearly  exact  numerical  inte¬ 
gration  of  reference  6  and  as  given  by  equation  (20)  are 
identical  to  the  nearest  foot  up  through  188,000  feet. 

GEOPHYSICAL  CONSTANTS 

In  order  to  select  a  consistent  set  of  geophysical  con¬ 
stants  a  brief  study  was  made  to  compare  the  values  of 
acceleration  of  gravity  as  computed  by  equations  (17)  and 
(19).  Equation  (19)  was  evaluated  by  using  the  following 
equation  for  the  effective  earth  radius  obtained  from 
reference  7; 

64.344882(1  -  .0026373 cos25n  +  .0000059 cos ^25^) 

R  =  _ y _ D 

3.085462x10’^  +  2.27xl0'^cos25p  -  2xlO~^^cos4Sp 

In  order  to  evaluate  equation  (17)  the  constants  in 
equations  (4)  anc  (16)  had  to  be  supplied.  The  values  of 
these  constants  as  shown  in  references  3,  6,  and  9  were 
substituted,  and  the  accelerations  from  equation  (17)  pro¬ 
duced  by  each  set  of  constants  were  compared  with  those 
from  equation  (19) .  The  set  of  constants  from  reference 


(22) 


11-21 


9  produced  the  closest  agreement  between  equations  (17)  and 
(19) .  These  values  were  the  ones  "most  used"  for  orbital 
calculations . 

The  values  from  reference  6  and  those  from  reference  3 
which  are  more  recent,  produce  accelerations  which  differ 
from  those  of  equation  (19).  If  the  values  from  reference 
3  or  later  values  are  to  be  used  in  equation  (4)  and  (16) , 
then  a  new  expression  for  the  effective  earth  radius  should 
be  derived  by  evaluation  and  differentiation  of  equation 
(17)  in  accordance  with  equation  (21) .  However,  it  is 

believed  that  the  constants  from  reference  9  are  adequate 
for  aircraft  trajectory  analysis.  These  constants  are  pre¬ 
sented  in  the  following  table. 


Constant 

Value 

Units 

r_ 

6,378,165 

meters 

0 

20,925,738 

feet 

3.986032x10“* 

3  : 

meters  per  sec 

1.4077768x10^® 

ft®  per  sec^ 

1/f 

298.30 

nondimens ional 

J2 

1082.30x10“® 

nondimens ional 

J3 

-2.3x10“® 

nondimens ional 

Jh 

-1.8x10“® 

nondimensional 

% 

7.2921x10"'* 

radians  per  sec 
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SUMMARY 


The  basic  assumptions,  definitions,  and  constants  which 
have  been  used  in  generating  model  atmospheres  are  presented. 
Such  an  atmosphere  provides  the  norm  to  which  all  test  data 
are  corrected.  Information  about  the  various  models  is  also 
presented,  together  with  the  concepts  of  geometric  and  geo- 
potential  altitudes. 
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Subscripts 


SYMBOLS  USED  IN  THIS  SECTION 
Definition 

acceleration  due  to  gravity 

local  effective  acceleration 
due  to  gravity 

reference  acceleration  of 
grav'd  ty 

sea  level  acceleration  of 
gravity 

geometric  altitude 
qeopotential  altitude 
temperature  gradient,  dT^/dK 
molecular  weight  of  air  or 
Mach  number 

molecular  weight  of  air  at 
sea  level 

ambient  pressure 

local  radius  of  the  earth 

universal  gas  constant 

ambient  temperature 

temperature  of  the  ice  point 
(273.15°K) 

m.olecular  scale  teiriperature 
air  density 

air  density  ratio,  p/Pq-^ 
geopotential 

base  of  atmospheric  layer 
sea  level 


Units 
f t/sec^ 
ft/sec^ 

f t/soc^ 

ft/sec^ 

ft 

ft 

“K/ft 

dimensionless 

dimensionless 

dimensionless 

Ib/ft^ 

ft 

ttVsec^^K 
deg  K 
deg  K 

deg  K 
slugs/ft^ 
dimiensionless 
f t^/sec^ 
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ATMOSPHERIC  ENVIRONMENT 


INTRODUCTION 

The  physical  characteristics  of  the  earth's  atmosphere 
vary  greatly,  changing  from,  day  to  day  and  v;ith  seasons  of 
the  year.  The  performance  of  an  aircraft  is  dependent  on 
the  nature  of  the  airmass  through  which  it  flies.  For 
example,  the  thrust  of  a  turbojet  engine  increases  appre¬ 
ciably  with  a  decrease  in  air  temperature.  Therefore  som.e 
set  of  standard  conditions  must  be  established  in  order 
for  performance  data  to  have  some  meaning  when  correlating 
data  from  one  flight  to  another  or  comparing  the  perform¬ 
ance  c£  one  aircraft  to  that  of  another.  V/hen  flight  test 
results  are  reduced  to  standard  conditions,  ideally  the 
corrections  applied  ought  to  be  as  small  as  possible  to 
minimize  errors  due  to  linear  approximations.  (Reference 
the  section.  Standardization  of  Excess  Thrust.)  This 
might  imply  the  use  of  a  reference  atmosphere  which  accu¬ 
rately  represf;.ited  the  mean  atmospheric  properties  at  the 
test  site  except  for  the  fact  that  reduction  of  test  data 
acquired  elsev.'hsre  in  the  world  for  purposes  of  comparison 
could  then  entail  large  corrections.  A  comprom.ise  between 
small  corrections  and  universal  applicability  requires  the 
use  of  an  idealized  middle-latitude,  year-round  mean  atmio- 
spheric  model. 
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MODEL  ENVIRONMENT 


BASIC  ASSUMPTIONS 


Since  the  earliest  models  of  the  atmosphere  were  pub¬ 
lished  (circa  1920)  the  same  assumptions,  with  minor  va' 
iations,  have  bev-n  made  in  the  range  of  0  to  20,000  meters 
(*65,600  feet) .  These  have  included; 

(1)  Sea  level  temperature  is  15*0 

(2)  A  constant  temperature  gradient  from  sea 
level  to  about  11,000  meters  (36,089  feet;)  ,  and  then 


(3)  A  constant  temperature  from  11,000  to  20,000 


meters 


Additionally,  it  has  been  assumed  that 

(4)  The  air  is  dry 

(5)  The  atmosphere  is  a  perfect  gas  so  that  the 
equation  of  state 


applies 


(6)  Hydrostatic  equilibrium  exists,  assuming  that 


the  atmosphere  is  static  with  respect  to  the  earth; 


-  pg  dh 


(2) 


By  co?nbinir.§  equations  (1)  and.  {2)  the  usual  differen- 
ti.c  ,  form  of  the  barometric  equatior-.  is  obtained 


If  g  is  taken  co  be  constant  and  T  is  replaced  by  a  lin- 

d. 

ear  function  of  h,  equation  (3)  can  be  integrated  quite 
simply  to  calculate  pressures.  This  procedure  was  fol¬ 
lowed  in  tabulating  the  older  atmospheres;  however.,  the 
assumption  of  constant  g  becomes  inadequate  as  altitudes 
are  extended  and  new  techniques  (discussed  in  following 
paragraphs)  were  devised. 

PRIMARY  CONSTANTS 

In  order  to  compute  properties  of  a  model  atmosphere, 
it  is  necessary  to  establi^.h  values  for  basic  constants 
appropriate  to  the  earth's  atmosphere.  To  illustrate, 
constants  which  have  been  used  to  form  tabular  values  of 
concern  in  aircraft  flight  test  have  been  extracted  from 
reference  4  and  appear  in  table  1.  For  a  more  complete 
list  of  constants  together  with  a  discussi.on  of  their 
origin,  see  reference  4,  pages  4  and  5. 
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Symbol 

Units 

P 

29,92126 

in,  Hg 

^SL 

0.076474 

lb  ft“^ 

15 

®C 

gf* 

32.1741 

ft  sec  ^ 

Ti 

273.15 

®K 

y 

1.40 

dimensionless 

R 

3089.80 

.2  -2  OI/-1 

ft  sec  °K 

UNIVERSAL  GAS  CONSTANT 

The  universal  gas  constant  in  the  perfect  gas  law  has 

the  dimensions  of  energy  mole*”^  (degree  absolute)  When 

divided  by  the  mass  of  one  mole,  the  gas  constant  of  air 

has  dimensions  of  energy  (unit  mass) (degree  absolute)"^. 

In  the  English  system  of  units  used  in  flight  test  work, 

the  energy  is  expressed  in  ft-lb^,  the  mass  in  Ib^^^'  and 

temperature  unit  in  *K  so  that  the  gas  constant  becomes 

ftrlb£].b_“'“K“*  or  ft  ft-sec"'^®K”* .  If  the  gas  law  is 
f  m  ^ 

written 

Pa  = 

R  has  tile  dijitiensions  ft  If  the  law  is  written 


*  The  a<3optecl  value  of  g^.  is  9.80665  meters  sec~^  which 

converged  to  the  English  system  of  units  is  32.jL7’4649  feet 
sec"’^.  The  above  value  from  the  UiS.  1962  Atmosphere, 
reference  4,  has  been  rounded  incorrectly. 
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Pa 

as  was  done  in  this  report  in  equation  (1)  then  R  has  the 
dimensions  ft^sec“^®K“ ' . 

SPEED  OF  SOUND 

The  speed  of  sound  is  by  definition  the  speed  of  propa¬ 
gation  of  the  wave  formed  by  an  infinitesimal  pressure  dis¬ 
turbance.  Such  a  disturbance  very  closely  approximates  an 
adiabatic  reversible  (isentropic)  process.  Considering  this, 
the  conservation  equations  for  mass  and  momentum  can  be  com¬ 
bined  to  give 


2 

a 


(iP) 

dp  s 


and  for  an  ideal  gas 


In  the  ARDC  1959  and  the  U.S.  1962  atmospheres  {references 
4  and  6)  y  has  been  taken  to  be  1.40  exact  (to  an  altitude 
of  90  kilometers) . 

RELATIONSHIPS  USED  TO  DEFINE  MODEL  ATMOSPHERES 

As  was  previously  pointed  out,  the  barometric  equation 
(equation  3)  can  be  integrated  easily  provided  that  a  con¬ 
stant  value  of  g  is  assumed.  This  was  done  in  the  older 
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atmospheres,  and  data  v/ere  tabulated  as  a  function  of  tape- 
line  altitude,  h.  At  the  high  altitudes  to  which  more 
recent  models  have  been  computed  this  assumption  is  no 
longer  valid.  Further,  variation  in  molecular  weight,  M, 
has  been  accounted  for  by  writing  the  perfect  gas  law  in 
the  form 

PgM 

^  "  RTg  (5) 

See,  for  example,  reference  4,  page  5,  and  reference  6, 
page  4.  Integration  of  the  resulting  barometric  equation 
becomes  quite  complex,  even  when  g  and  M  are  replaced  by 
very  simple  functions  of  h.  To  retain  the  mathematical 
simplicity  of  the  equations  used  for  low  altitude  calcu¬ 
lations,  two  transformations  of  variables  have  been  made. 
The  two  new  parameters  are;  geopotential  altitude,  H,  from 
combining  g  and  h,  and  molecular-scale  temperature,  Tj^^, 
from  combining  T^  and  M.  By  defining  1!^  as  a  series  of 
linear  functions  of  H,  integration  may  be  carried  out 
exactly  as  for  equations  using  geometric  altitude  and  con¬ 
stant  g. 
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Molecular-Scale  Temperature 

By  definition  the  molecular-scaJ.e  temperature  is 


Tm 


Mot 

M  a 


(6) 


At  extreme  altitudes  is  widely  different  from  the 
kinetic  temperature  (reference  4,  page  8);  however,  since 
in  the  altitude  range  of  interest  in  aircraft 
flight  test,  no  further  consideration  is  given  to  this 
parameter. 

Geopotential  Altitude* 

The  geopotential  at  an  altitude  h  is  the  potential 
energy  of  a  unit  mass  at  that  altitude  relative  to  the 
potential  energy  of  the  same  mass  at  mean  sea  level.  In 
differential  form  geopotential,  if  is  related  to  geometric 
(tapeline)  altitude  by 

=  gj^dh  (7) 

The  force  of  gravity  in  equation  (7)  is  the  resultant  of 
two  forces:  (1)  the  gravitational  attraction  which  derives 
from  Newton's  universal  law  of  gravitation,  and  (2)  the 
centrifugal  force  (commonly  called  centrifugal  relief)  in 

*  Tbe  following  information  on  geopotential  altitude  per¬ 
tains  to  model  atmospheres.  Additional  information  may  be 
found  in  the  section.  Geophysical  Properties. 
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a  refexvince  frame  attached  to  the  earth,  assuming  that  the 
atmosphere  rotates  with  the  earth.  Hence,  a  local  value 
of  g  is  a  function  of  both  vertical  displacement*  above 
the  earth's  surface  and  latitude. 

Integrating  equation  (7)  the  geopotential  at  an  alti¬ 
tude  of  h  is 

h 

^  =  /  gi  (8) 

0 


In  accordance  with  equation  (8)  the  geopotential  alti¬ 
tude  is  defined  as 


H 


/"fkdh 

0  gf 


(9) 


By  introducing  the  reference  gr^avity,  g^,  geopoten¬ 
tial  altitude  has  the  dimensions  of  length  (geopotential 
feet  in  the  English  system)  ,  and  is  equivalent  to  '!>,  which 
is  bhe  amount  of  work  done  in  raising  a  unit  mass  from 


*'  In'  the  U.S.  1962  Atmosphere  a  refinement  is  made  in  that 
h  is  measured  along  the  line  of  force  through  the  point, 
from  the  equipotential  surface  for  which  |J>  =  0  to  the 
point  in  question,  causing  the  lines  of  force  to  curve 
toward  the  poles.  The  difference  in  distance  according 
to  this  definition  and  a  straight  line  distance  is-  negli¬ 
gible  within  the  sensible  atmosphere  and  is,  therefore,  of  j 
no  consequence  in  aircraft  fliaht  testing.  ■ 
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mean  sea  level  to  a  geometric  altitude  of  h,  A  value  of 
9.80665  meters  sec*"^  (32. 17405  feet  sec"^)  was  adopted  for 
the  ICAO  Standard  Atmosphere  and  for  the  ARDC  1959  Atmo- 
sphere.  The  same  value  was  also  adopted  for  the  U.S. 
Standard  Atmosphere,  1962  for  a  geographic  latitude  of 
exactly  45®.  It  should  be  noted  that  physical  displace¬ 
ments  between  equipotential  surfaces  separated  by  a  con¬ 
stant  amount  in  terms  of  geopotential  is  not  constant. 
Rather,  the  physical  displacement  increases  v/ith  altitude 
because  of  the  decreasing  values  of  g. 

In  differential  form  equation  (9)  is 

gfdH  = 

Substituting  equation  (10)  ,  equation  (3)  may  be  written 
as 


<1(1"  ly  -  ,11 

J.  Ig 

In  order  to  integrate  this  equation  tc  find  geopoten¬ 
tial  altitude  as  a  function  of  pressure,  H(T  )  is  substi- 
tuted.  To  simplify  the  procedure,  it  is  assumed  that  the 
atmosphere  is  made  up  of  layers  in  which  the  temperature 
varies  linearly  with  geopotential  altitude  (constant 
temperature  gradient) . 
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The  general  temperature-altitude  relationship  is  then 


I'a  ■  <Ta)b+ 


(12) 


Substituting  equation  (12)  in  equation  (11)  and  integra¬ 
ting  produces  the  following  equation: 


InP,  .  Ln|(H-"i)l 

M 


(13) 


Taking  (P_)v,  to  be  the  pressure  at  the  base  of  the  layer 

o  ID 

we  have 


Pa  “  (Pa)b[ 


<Ta)b 


Sr./RL 


(T3)b..^j(H-Hb) 


M 


(14) 


when  the  temperature  gradient,  is  not  zero.  Following 
the  same  steps  with  Lj^  =  0  leads  to 


Pa  =  (Pa)b  exp[--_— - ] 


P^Ta)b 


(15) 


Properties  which  appear  in  the  above  equations  have 
been  taken  from  reference  4  and,  after  converting  units, 
are  presented  in  table  2. 
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Table  2  -  Properties  Defining  the  U,S 
Standard  Atmosphere/  1962 


Altitude,  H 

Gradient,  L 

Atmospheric 
Pressure,  Pg 

Atmospheric 
Temperature,  Tg 

km 

feet 

“K/foot 

in.  Hg 

“K 

0.000 

-1.98120x10"^ 

29.92126 

288.15 

11.000 

36,039.24  ■ 

0 

6.68321 

216.65 

20.000 

65,616.80 

0.30480x10"^ 

1.61671 

216.65 

32.000 

104,986.88 

0.85344x10"^ 

0.25632 

228.65 

47.000 

154,199.48 

0 

0.032750 

270.65 

52.000 

170,603.67 

-0.60960x10"^ 

0.017423 

270.65 

61.000 

200,131.23 

0.0053773 

252.65 

Pressure,  Temperature,  and  Density  Ratios 

From  the  preceding  equations,  general  equations  defin¬ 
ing  pressure,  temperature,  and  density  ratios  may  be  found 
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Prom  equation  (14) 


(Pa).  (TA  i 


when  0  and,  from  equation  (15) 


(Pa)b  ^-gr(H-Hb), 

exp[ - - ) 


^SL 


R^Ta)b 


when  L„  =  0, 

M 

The  general  equation  for  temperature  ratio 
tion  (12)  is 


^SL 


From  the  perfect  gas  law  6 
(16)  and  (18) 


00  so  that  from 

a 


o 


^'^a)b-^Hl(H-Hb) 


+  I 


when  Lj^  7^  0  and,  from  equations  (17)  and  (18) 


(16) 


(17) 


from  equa- 


(18) 


equations 


(19) 
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a 


(20) 


=  ff,  exp[ 


-  gr  (H-Hb) 


when  L„  =  0 
M 


Geometric  Altitude 


In  flight  test  applications  the  calculation  of  geomet¬ 
ric  (tapeline)  altitude  is  frequently  required.  Geometric 
altitudes  cannot  be  computed  directly  for  off-standard 
conditions;  an  integration  procedure  must  be  resorted  to 
using  continuous  profiles  of  pressure  and  temperature. 

From  dP  /P  =  d(lnP  )  equation  (3)  becomes 

ci  Gi  H 


ddnPg) 


_!k 

RT, 


dh 


(21) 


From  the  inverse-square  law  of  gravitation  (reference 
equation  (19)  in  the  section.  Geophysical  Properties) 


^SL 


R, 


(Re  +  h)- 


(22) 
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Substituting  equation  (22)  in  equation  (21)  and  rearranging 


dh 


JL 

SSL 


Ta( 


£e. 

R 


+”h  2 

— )  ddnP^) 


(23) 


Equation  (23)  may  be  written  in  integral  form  as 


h, 


hn 

/  dh 
n-1 


R  ,  j!k 

SSL  2  Re  Pa 


(24) 


with  the  assumptions  that  the  temperature  is  constant  at 


=*11-1 


and  that 


^  s,  ^n-1 

Fg  Rg 


Integrating  equation  (24) 


h 


n 


Cl 


*n-l 


+  ^  1  ^n-1  ^  ^^n-1 

^SL  2  R3  P 


(25) 
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Geometric  altitudes  are  found  by  summing  increments  as 
indicated  by  equation  (25) . 

STRUCTURE  OF  THE  ATMOSPHERE 

The  atmposphere  has  been  broken  up  into  four  major 
regions  which  are  associated  with  physical  characteristics 
listed  in  table  3.  The  names  of  atmospheric  shells  and 
boundaries  used  in  the  table  have  been  adopted  by  the 
World  Meteorological  Organization. 


Table  3  -  Description  of  Atmospheric  Shells 


Name 

Description 

Troposphere 

The  region  nearest  the  earth's  surface  having  a 
uniform  decrease  in  temperature  with  altitude. 

The  troposphere  is  the  domain  of  weather  where 
turbulence  caused  by  convection  occurs.  At  the 
tropopause  (top  of  the  troposphere)  high  winds 
are  common  and  the  highest  cirrus  clouds  are 
found. 

Stratosphere 

The  region  above  the  troposphere  having  a  con¬ 
stant  temperature  followed  by  increasing  tempera¬ 
tures,  reaching  a  maximum  at  the  stratopause. 
Maximum  of  atmospheric  ozone  is  found  near  the 
top  of  this  region.  Turbulence  is  infrequent. 

Mesosphere 

Temperature  remains  constant  with  altitude  above 
the  stratopause  and  then  decreases.  This  region 
is  in  radiative  equilibrium  between  ultraviolet 
ozone  heating  by  the  upper  fringe  of  ozone  region 
and  the  infrared  ozone  and  carbon  dioxide  cooling 
by  radiation  to  space. 

Thermosphere 

The  region  of  rising  temperature  above  the  major 
temperature  minimum.  No  upper  altitude  limit. 
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COMPARISON  OF  ATMOSPHERIC’ ^MODELS  li;! -  ^CENT  iUp 
Standard,  atmospheres  -which  have  :beeh  used 
years  are:  - 

(;bj  International  Civil  Aviation  Ori^hizatign 
Standard  A tmosphere^  1952  adopted  by  ihe ^National  Advisory- 
Committee  for  Aeronautics  and  pubiishedt  in  NACA  Report  1:235 
(.2)  Air  Research  and  Development  iConim'and  M 
Atmosphere,-  1956 

(3) .  U.S.  Extension  to,  the  ICAO  Standard' Atmosphere , 

1958 

(4)  ARpC  Model  Atmosphere^  1959  «  * 

i5)  U.SV  Standard  Ai^dspher^^  1962 

The  International  Ci'yi'l  Aviatidjas-Organizati^^n;^  Staiid^d 
Atmosphere  was-  adopted: ^by  .thejiNACA  ^pn  .Ndyembe'r  :2Q>.  i952  and 
is  contained  iii  NA£a  Re'ppife;  .i23"5> :  Standard'-^A'ttiosphare  ~ 
Tables  and  Data  jor  Alii-^tudes.io  65>i80b-  Feet.  .The;  eguations_^ 
.pi  this  rjeppr.t  v;ere,- used  to  extend  SpVpOO  feet 

andappehr  ,,in  .refP^ipnce:  2  To  an;  alt^^^^  ‘86>6oo  $‘ept> 

'tha  iproplpties.  ;s'uch'  as  e>  a  irPm;  NACa  *RepPrt.  123S 

.are  idehtical  with  ihosje  in  the  .‘J^C  ^Hpdel;  Atmosphere ,  19’56;. 
the  UiS.  ExtensiPn  ^to  the  iCAO  . Standard  Atmosphere ,  ,195i^8:  and; 
the  ARDC  Model  .Atmosphere/  19.6.9;.  '  .  - 

Some  important  changes -were:  imade  in  the-'-yvS'i  'Standard;'. 
Atmosphere,  1962.  The  principal 'difference- between,  the  ;y,S«, 


l9.65^and^  previous  models  is  that  the  altitude  at  the  top 
"Of:  the  stratosphere  is  65,616.8  geopotential  feet  in  the 
fonnei"and"^2, 021.0  feet  in  (2),  (3),  and  (4)  from  the 
list,  above..  Above  65,616.8  feet,  the  U.S.  1962  Atmosphere 
has  a  -temperature  gradient  of  0.3048®K/1000  feet  as 
depicted  in  figure  1. 


comparison  of  Temperature  Profiles  in 
Standard  Atmospheres 

Figure  1 


H!»21 


The  temperature  in  the  U.S.  1962  atmosphere  is  4.4* 
Kelvin  higher  at  an  altitude  of  80/000  feet;  thus,  an 
airspeed  of  700  knots  (for  example)  at  this  altitude 
would  represent  M  -  1.208  in  the  1962  atmosphere  and 
M  =1.220  in  the  others.  The  difference  in  temperature 
produces  differences  in  tabulated  properties  at  altitudes 
over  65,000  feet  (e.g.,  6,  6,  etc). 
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SECTION  IV 


PLIGHT  PARAMETERS 
FROM  SENSED 
ENVIRONMENT  (ON-BOARD 
STATIC  AND  TOTAL  PRESSURE, 
AND  TOTAL  TEMPERATURE) 


SUMMARY 


The  information  in  this  section  on  airspeeds,  tempera¬ 
tures,  etc.,  is  similar  to  that  found  in  other  documents  (e.g., 
AF  Technical  Report  6273,  Flight  Test  Engineering  Handbook) , 
but  has  been  included  for  the  sake  of  completeness  and  to  have 
a  source  of  basic  equations  for  use  in  subsequent  sections. 

In  addition,  an  examination  has  been  made  of  the  effects  of 
high  speed  on  the  usual  assumptions  that  air  obeys  the  equation 
of  state  P  =  pRT  and  has  a  constant  ratio  of  specific  heats. 
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SYMBOLS  USED  IN  THIS  SECTION 


Symbol 

Definition 

Units 

a 

speed  of  sound 

ft/ sec 

Cl 

airplane  lift  coefficient 

dimensionless 

g 

acceleration  due  to  gravity 

ft/sec ^ 

h 

enthalpy 

BTU/lb 

j 

mechanical  equivalent  of  heat 

ft-lb/BTU 

K 

temperature  probe  recovery 
factor 

dimensionless 

M 

flight  Mach  number 

dimensionless 

pc 

correction  to  Mach  number  for 
position  error 

dimensionless 

load  factor  along  the  Z“axis 

dimensionless 

Pa 

ambient  pressure 

Ib/ft^ 

P. 

1 

Pg 

indicated  pressure 

Ib/ft^ 

static  pressure 

Ib/ft^ 

Pt 

total  pressure 

Ib/ft'- 

^Pp 

static  pressure  source  position 
error 

Ib/ft^ 

impact  pressure 

Ib/ft^ 

R 

universal  gas  constant 

ft^/sec^'^K 

S 

wing  area 

ft* 

""a 

ambient  temperature 

deg  K 

total  temperature 

deg  K 

calibrated  airspeed 

ft/sec 

com.prcssibility  correction  to 
calibrated  airspeed 

ft/sec 
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Units 


Symbol 


V 


V. 


AV_ 


pc 


W 

z 

a 

Y 

'a 

X 

y 

p 

Subscripts 

ic 


s 

SL 

t 


Definition 
equivalent  airspeed 
true  airspeed 

correction  for  airspeed  position 
error 

airplane  gross  weight 
compressibj lity  factor 
angle  of  attack 

ratio  of  specific  heats:  Cp/C^ 
ambient  pressure  ratio 
pressure  lag  constant 
absolute  viscosity 
air  density 

indicated  corrected  for  instrument 
error 

corrected  for  lag 
static 
sea  level 
total 


ft/sec 

ft/sec 

ft/sec 

lb 

dimensionless 

rad 

dimensionless 

dimensionless 

sec 

Ib-sec/ft^ 

slugs/ft^ 


IV-5 


INTRODUCTION 


Pressure  altitude,  airspeed,  Mach  number,  and  free  air 
temperature  are  basic  parameters  in  the  performance  of  air¬ 
craft.  Conventional  instruments  used  to  measure  these 
quantities  are  the  altimeter,  the  airspeed  indicator,  the 
machmeter,  and  the  free  air  temperature  probe.  (Mach  num¬ 
bers  deduced  from  altimeter  and  airspeed  readings  are  pre¬ 
ferred  to  those  from  machmeters.)  Relationships  of  the 
basic  parameters  to  environmental  conditions  sensed  on 
board  an  aircraft  (static  and  total  pressures,  and  total 
temperature)  are  developed  in  this  section.  More  compre¬ 
hensive  derivations  of  equations  and  descriptions  of  tlie 
construction  and  calibration  of  instruments  may  be  found 
in  reference  2.  It  should  be  noted  that  in  this  reference 
the  usual  simplifying  assun\ptions  are  made  (e.g.,  constant 
ratio  of  specific  heats) .  These  assumptions  lead  to 
errors  in  calibrated  air  data  when  Mach  numbers  get  much 
above  two.  Real  gas  effects  as  they  influence  calibrated 
air  data  together  with  equations  which  may  be  used  at  high 
supersonic  and  hypersonic  speeds  are  presented  in  this 
section. 
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pbessure  altitude 

Fror,  the  general  equations  relating  pressure  and  alti¬ 
tude  (equations  (18)  and  (19)  in  the  section,  Atmosphor^ 
Environment)  we  have,  after  substituting  constants  appro¬ 
priate  to  the  U.S.  Standard  1962  atmosphere  from  table  2,  for 

«a  “ 


.  (1-6.87558x10-‘H)'^®® 

Tor  - 16,404.20  S  H  S  36,089.24  gcopotenda)  feel 


.  0.223360  exp(- 4.80637  rl0'^(  II -36059.M)  I 

for  36,089.24  <  H  <  65,616.80  geopotential  feet 

Sjj  -  0.0540322[1  +  1.40688 xl0~^’(n- 65.61 6.80) 

for  65,616.80  <,  H  <  104,986.88  geopotential  feet 


0.0085664911  +  3.732S’xlO‘^(H  -  104,986.88)]' 


-12.2012  (4) 


for  104,986.88  <  H  <  154,199.48  geopolential  feet 

Altimeters  are  built  and  calibrated  accoroing  to  these 
relationships  (or  perhaps  to  another  model  atmosphere) . 
Differences  in  altitude  between  the-  various  model  atmo¬ 
spheres  are  fairly  small  (rcferenc'^  1,  pages  13  -  15); 
however,  the  standard  used  for  instrument  calibration 
should  be  known.  If  laboratory  calibrations  are  made  using 
a  model  atmosphere  other  than  the  one  desired,  corrections 


be  easil'i  imade;  through  the  gerierat  pres'sxixieHai^"^^ 
equations  and  ;prp'perties'  def ihing=  -the  i|wo  tatfixpspiveres . 

The  static  pressure  s^ehsed  at  the  spatic 'sou^^ce  the 
altimeter.,  P  may  differ  slightly  from,  the  atmospheric 
pressure,  P.  because  of  pressure'  Mg  arid  positiph  error. 
These  errors  are  discussed  in  later  ^paragraph . 

-.GhLIBPATED  AIRSPEED 

The  airspeed  indicator  senses  the  difference  between 
•two  pressures;  total  pressure,  and- static  pressure, 

P  .  The  difference  in  pressure  .is  'conyerted  to  a  speed 

O.  '  V 

through  Bernoulli 's  compressible  equation  for  frictioniess 
adiabatic  (isentropic)  flow  in  which  airspeed  is  expressed 
as  the  difference 'between  total  and;  static,  pressures . 
Assxoming.  that.  P^  -  P^,  Bernoulli '>s  equation  -may  be 
expressed  as 

^  .  2^  .  [1*  _  1  (5) 

2,  '  a 

for  subsonic  speeds; 

VJith  Y  =  1.40,  equation  (5)  becomes 


u  +  o.2(hh^'^ 


(6) 


At  supersonic  speeds  a  detached  shock  wave  forms  in 
front  of  the  total  pressure  probe,  and  equation  (5)  is  no 


longer  valid.  In  this  case  the  Rayleigh  supersonic  pitot 


formula 


y  1 

ry-t-l/\)2]^r  y-hl  ^  J 

1-2 


relating  total  pressure  behind  the  shock  to  the  free  stream 
ambient  (atmospheric)  pressure  is  used.  It  should  be  noted 
that  =  P^'  -  where  the  total  pressure,  sensed 

at  the  pitot  head  is  behind  the  shock  and  is  not  equal  to 
the  free  stream  total  pressure,  P^,  at  supersonic  speeds. 
Substituting  1.40  for  y  and  simplifying,  equation  (7)  becomes 

_  '  166.921 

2c  .  — - — J - 1  .  (8) 

’'"a  [7(1)2-11^*^ 

a 

Examination  of  the  above  equations  shows  that  true  sp'eed, 
V^,  is  dependent  on  the  speed  of  sound,  a,  and  atm.ospheric 
pressure,  P^,  in  addition  to  the  differential  pressure,  q^. 
Therefore,  an  airspeed  indicator  measuring  differential 
pressure  can  be  made  to  read  true  airspeed  at  only  one  set  of 
atmospheric  conditions.  Sea  level  standard  has  been  selected 
arbitrarily,  and  the  dials  of  airspeed  indicators  are  scaled 
so  that  a  given  differential  pressure  will  indicate  a  speed 
in  accordance  v;ith  equations  (6)  and  (8)  in  v'hich  sea  level 
standard  and  P  are  inserted.  This  sea  level  standard  value 
of  V.  is  defined  as  calibrated  airspeed,  V  ,  and  is  found 
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from  airspeed  indicator  readings  after  corrections  have 
been  rcade  for  instrument  and  position  errors. 

Accordingly,  equations  (6)  and  (8)  may  be  rewritten  as 


il-  -  [1  +  0.2(]^)2]^*^  ~  1 


^SL 


^SL 


and 


^SL 


166.921 

t7(lL)=- 

®SL 


-  1 


(9) 


(10) 


EQUIVALENT  AIRSPEED 

The  equivalent  airspeed,  V^,  is  the  result  of  correct¬ 
ing  the  calibrated  airspeed  for  compressibility  effects. 

The  airspeed  indicator  is  calibrated  to  read  correctly  at 
standard  sea  level  conditions  and  thus  has  a  compressibil¬ 
ity  correction  appropriate  for  these  conditions.  VJhen  an 
aircraft  is  operated  at  altitude,  however,  the  compensation 
is  inadequate  and  an  additional  correction  must  be  applied. 

Hence,  V  and  V  are  related  by 
c  e 

=  Vg  +  AVp  (11) 


Equivalent  airspeed  coupled  with  standard  sea  level 
density  produces  t).ie  same  dynamic  pressure  as  the  true  air¬ 
speed  (speed  relative  to  the  airmass)  coupled  with  the 


AYj.  (knot*) 


actual  air  density  at  flight  conditions.  Therefore 


i;  2  V  ^ 

\P  "  \  PS'L 


(12) 


or 


(13) 


For  subsonic  flight,  solving  equation  (5)  for  , 


-  5a2[(^  ~  1] 

‘  Pa 


(14) 


Combining  equations  (4)  and  (5)  from  the  section,  Atmo¬ 
spheric  Environment,  the  speed  of  sound  in  a  perfect  gas 
may  be  expressed  as 


a  .  (I^k)'/^ 

P 


(15) 


Replacino  V. ^  by  V^^/a  and  a^o  by  yV  /p  . ; 

u  6  SL  oJu 


^SL  Pa 


(16) 


Solving  equation  (6)  for  calibrated  airspeed: 


7R 


«SLr/ic.  ,i2/7  ,,,1/2 


V,  -  +  nr 

•  a 


^S'L 


(17) 
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j  The  difference  between  equations  (16)  and  (17)  is  the  com¬ 

pressibility  correction,  AV^,  in  equation  (11) .  If  = 

*  P  calibrated  and  equivalent  airspeeds  are  the  same  and 

^SL 

AV  is  zero  rcaardless  of  the  macnitude  of  q_.  Also,  AV 
[  .  c  •'  •  -*0  c 

can  be  calculated  for  any  altitude  and  calibrated  airspeed 

since  a  value  of  V,  determines  q  . 

c  c 

MACK  NUMbER 

Mach  number  is  defined  as  the  ratio  of  true  airspeed 
to  idle  local  speed  of  sound: 

V 

M  -  -1  (18) 

a 

From  Bernoulli’s  equation  with  isentropic  flow  of  a  perfect 
gas 

y 

A  -  (I+IzlIm^)^"^  (19) 

This  equation  relates  Mach  nuirber  to  free  strearr.  total 
and  static  pressures  and  is  applical>le  for  supersonic  as 
well  as  subsonic  flight.  It  should  be  remendjered,  however, 

'  that  P^'  rather  than  P^  is  sensed  by  a  pitot  probe  in  super¬ 

sonic  flight. 

Substituting  1.40  for  y,  equation  (19)  becomes 

A  .  (1+0.2M^)^‘^  (20) 
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iiaaBM»gaa«aatiiaaig»a 


The  machmeter  equation  for  subsonic  flight  is  found 
by  substituting  the  definition  for  Mach  number  in  equation 


A  1=  [1  +  «  1 


Solvina  for  M 


y-i 


With  Y  =  1.40 


For  supersonic  flight  from  equation  (7) 

y  _L 

-5c  _  ^ 5  _ 


)  -  1 


1~  >'  +  2yM" 


This  equation  cannot  be  solved  explicitly  for  Mach  number. 
It  can,  however,  be  put  in  the  form 


(3  +l)(y +  1)  +  y  -  1 


9y(y  +  lM^  y  +  1 

^  l-y  +  2yM^ 


which  may  be  solved  by  iteration. 
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Substituting  1.40  for  y  and  rearranging 


[ 0.7766628 (^+  1)(1  - 


.sviaL  temperaturl 

If  the  air  surrounding  a  temperature  probe  is  brought 
to  a  complete  stop  adiabaticnlly  the  resulting  temperature/ 
if  sensed  correctly  is 

T,  =  +  (27) 


For  various  reasons,  such  as  radiation  or  heat  leakage, 
most  probes  do  not  register  the  full  adiabatic  temperature 
rise.  Introducing  a  probe  recovery  factor,  K,  the  equation 

Tt  =  Ta(l+K21^M^)  (28) 

may  be  v/ritten  to  account  for  a  lack  of  complete  adiabatic 
temperature  rise.  The  magnitude  of  K  is  between  0.95  and 
1.00  for  most  installations  and  can  often  be  assup’ed  con¬ 
stant.  Variations  v/ith  altitude  and  Mach  number  should  be 
expected,  howev-'^r.  particularly  at  supersonic  speeds. 
Methods  for  determining  K  for  a  given  installation  are  dis¬ 
cussed  in  reference  2. 


ri'1.1 BRTVTFD  AIR  DATA  AT  JilGH  SPEEDS 

For  flight  at  speeds  below  Mach  2  it  is  generally  ade¬ 
quate  to  assume  that  air  behaves  as  an  ideal  gas  obeying 
the  equation  of  state.  P  =  pRT  and  having  a  constant  ratio 
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t 


of  specific  heats  taken  to  be  y  =  1.40.  However,  at  speeds 
slightly  above  Mach  2  and  for  higher  Mach  numbers  this 
assumption  becomes  invalid.  In  general  it  is  necessary 
to  account  for  real-gas  effects  v;here  the  specific  heats 
are  functions  of  temperature  and  pressure,  and  the  equa¬ 
tion  of  state  must  be  expanded  to  include  the  effects  of 
intermolecular  attraction  at  high  densities  and  dissoci¬ 
ation  and  ionization  at  high  temperatures  and  low  pres¬ 
sures.  A  general  form  of  the  equation  of  state  includes 
tile  compressibility  factor,  Z,  as  follows 

P  =  pZRT  (29) 

Viithout  the  assumption  of  constant  specific  heats  it  is 
impractical  to  obtain  an  analytic  expression  for  param¬ 
eters  such  as  the  total  temperatures  and  pressures  in 
front  of  and  behind  the  shock.  A  method  has  been  employed 
to  compute  these  parameters  by  interpolating  air  thermo¬ 
dynamic  properties  from  a  comprehensive  set  of  tables 
(reference  3)  and  computing  the  total  and  behind  the  shock 
conditions  using  the  conservation  equations  of  mass, 
momentum,  and  energy  as  follows: 

92%  =  PiVi  (30) 

P2  +  P2'^2  =  (31) 
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(32) 


h 


t 


The  free  stream  temperature,  pressure,  and  density  are 
defined  by  the  am.bient  conditions  at  a  given  altitude,  and 
the  ambient  enthalpy  and  entropy  can  be  interpolated  from 
the  thermodynciir.ic  tables  of  reference  3.  The  total  enthalpy 
at  each  Mach  number  is  given  by  equation  (32) ,  and  by  defi¬ 
nition  the  total  (isentropic)  temperature  and  pressure,  T^ 
and  P^,  can  be  interpolated  at  this  total  enthalpy  and  the 
ambient  entropy. 

The  static  properties  behind  the  shock  can  be  obtained 
by  an  iterative  solution  of  equations  (29) ,  (30)  ,  (31) ,  and 
{32),  and  then  the  entropy  behind  the  shock  can  be  interpo¬ 
lated  from  the  thermodynamic  ta}:)les.  The  total  temperature 
and  pressure  behind  the  shock,  T^'  and  can  then  be 
interpolated  from  the  same  tables  as  a  function  of  the 
total  enthalpy  and  the  behind-the-shock  entropy. 

The  values  of  T^,  P^,  T^' ,  and  P^'  obtained  by  this 
method  will  generally  differ  from  the  values  given  by  equa¬ 
tions  (20),  (24)  and  (27).  However,  these  differences  in 
no  way  affect  the  physical  operation  of  the  airspeed  indi¬ 
cator.  Consequently,  the  definition  of  differential  pres¬ 
sure  the  equations  for  calibrated  air¬ 

speed,  equations  (9)  and  (10)  are  still  applicable.  As  an 
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illustrafcion  of  the  real  gas  effects  on  the  values  of 
and  a  comparison  with  the  values  from  the  ideal  gas 
equations  is  presented  in  figure  2.  The  real  gas  tempera¬ 
ture  ratio  diverges  from  the  real  gas  value  immediately 
above  Mach  2.  The  pressure  remains  within  plus  or  minus 
one  percent  up  past  Mach  3  but  diverges  beyond  this  point. 
Tlie  differences  in  calibrated  airspeed  arc  not  as 
largo,  being  much  less  than  one  percent  up  to  Mach  5. 

Figure  2  can  be  used  as  a  tool  to  decide  when  the 
real  gas  method  described  here  should  be  used  for  data 
reduction.  Computer  programs  have  been  developed  to  per¬ 
form  these  calculations,  and  a  set  of  tables  of  tempera¬ 
ture  and  pressure  ratios  and  calibrated  airspeeds  has 
been  produced.  These  arc  described  in  reference  4. 

ERRORS  IH  PRESSURE  MEAgUREMEMT 

In  addition  to  the  usual  instrument  error,  altimeters 
and  airspeed  indicators  are  subject  to  two  additional 
errors.  They  are  position  error  and  pressure  lag  error 
and  are  discussed  in  the  following  paragraphs. 

POSITION  ERROR 

To  determine  the  speed  and  altitude  at  which  an  air¬ 
craft  is  flying,  values  of  dynamic  and  amljient  pressure 
are  required  as  indicated  by  the  preceding  equations. 


The  pressures  rectistered  By  a.  prtot-static  system  will, 
in  general,  differ  from  the  desired  values.  Distortions 
in  flow  are  caused  by  the  presence  of  t),ie  aircraft  and  by  the 
system  itself;  pressure  sensors  are  located  in  a  flov.'  field 
which  is  different  from  the  flov;  field  distant  from  the  air¬ 
craft.  The  resultino  errors  are  called  position  errors. 

Total  Pressure  Error 

At  subsonic  speeds  the  flov*;  perturbations  due  to  the  air¬ 
craft  static  pressure  field  are  very  nearly  is en tropic-  and 
hence  do  not  affect  the  total  pressure.  Therefore;,  as  long 
as  tlie  to’ al  pressure  source  is  not  located  in  a  v;irig  v/ake, 

in  a  boundary  layer,  or  in  a  region  of  local  superson5.c  flov’, 

..  ** 

the  total  pressure  error  due  to  the  position  of  the  total  pres¬ 
sure  head  v/ill  usually  be  negligible; 

Kose  boom  pitot-static  systems  are  installed  on  supersonic 
aircraft  so  that  the  total,  pressure  pickup  v;iil  be  located 
ahead  of  any  shock  v;aves  formed  by  the  aircraft.  The  shock 
wave  due  to  the  pickup  itself  is  accounted  for  by  the*  equation 
against  which  airspeed  indicators  are  calibrated  (reference 
equation  \10)). 

Failure  of  the  total  pressure  sensor  to  register  the  local 
pressure  may  result  from  the  shape  of  the  head,  inclination  of 
the^  flow,  or  burred  or  misshapen "pitot  lips.  Since  pitot-static 
probes  are  tested  in  wind  tunnels  before  intallation  to  assure 
good  design  and  commonly  used  probes  produce  no  significant 
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errors  due  to  inclination  to  the  relative  wind  up  to  approxi¬ 
mately  20  degrees,  there  should  be  no  significant  errors  in 
total  pressure  measurement.  In  flight  test  applications  it 
is  usually  presumed  that  all  of  the  position  error  originates 
in  the  static  pressure,  source.  The  possibility  of  a  total 
pressure  -error  must,  however,  be  considered,  and  airspeed 
calibrations  should  be  investigated  to  find  if  position  errors 
in  total  pressure  do  exist. 

Static  Pressure  Error 

The  static  pressure  field  surrounding  an  aircraft  in 
flight  is  a  function  of  speed  and  altitude  as  well  as  the 
secondary  parameters;  angle  of  attack,  Mach  nurri>er,  and 
Reynolds  number.  Kence,  it  is  seldom  possible  to  find  a, 
location  for  the  static  .pressure  source  where  the  free  strear.''. 
pressure  will  be  sensed  under  all  flight  conditions.  There¬ 
fore,  an  error  in  the  measurement  of  the  static  pressure  due 
to  the  position  of  the  static  pressure  source  in  the  aircraft 
pressure  field  will  generally  exist. 

At  subsonic  speeds  it  is  often  possible  to  find  som.e 
position  on  the  aircraft  fuselage  v/here  the  static  pressure 
error  is  small  under  all  flight  conditions,  and  flush  static 
ports  are  usually  installed  on  the  fuselages  of  subsonic  air¬ 
craft.  Nose  booms  are  generally  installed  on  supersonic 
aircraft  to  minim.ize  the  possibility  of  total  pressure  error; 
static  pressure  sensed  on  a  boom  has  the  advantage  that  in 
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supersonic  flight  the  bov.’  x^ave  formed  by  the  aircraft  is.  down¬ 
stream  of  the  static  pressure  ports  so  that  the  pressure 
measured  is  unaffected  by  pressure  field  of  the  aircraft. 

Not  considering  Reynolds  number  and  assximing  that  side¬ 
slip  angles  are  small,  the  functional  statement 


f(M,  a) 


(33) 


may  be  written. 


Defining  the  position  error. 


as 


(34) 


equation  (33)  may  be  modified  to 


AP 

-LE  «  f(M,  a) 


(35) 


Since  q  /P  is  related  to  Mach  number  only,  equation  (35) , 
in  terms  of  indicated  values  corrected  for  instrument  error, 
can  be  changed  to 


AP„ 


fljc) 


(36) 


1C 


or 


=  f(Mic,€L.  ) 


(37) 
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For  constant  n„VJ 
z 


-  fCVjc)  only 


(43) 


1C 


Prom  subsequent  derivations  of  equations  (67)  and  (70)  it 
can  be  seen  that 


AV 


pc 


(44) 


Since  AP  /a  =  f(V.  )  only  in  the  lov?  I-3ach  number  ranre 
P  ic' 

from  equation  (43) 


^Vpc  -  f(Vic)  only 


(45) 


for  constant  n„W  in  the  absence  of  Mach  num])er  effects, 
z 

At  higher  sneeds,  when  there  may  be  both  M,.  and  Cb 

"  1 C  Ij  » 

1C 

effects,  airspeed  calibrations  may  be  put  in  the  form 

AP  /q^  or  am  versus  Data  from  nose  booms  or  wing 

10 

tip  booms  will  usually  form  a  single  line  for  Mach  numbers 
greater  than  about  0.6.  At  lower  speeds  variation  in 
may  be  of  consequence  especially  for  large  airplanes  which 
have  large  changes  in  weight. 

For  aircraft  equipped  with  nose  booms  the  static  pressure 
error  increases  very  rapidly  with  Mach  numiber  in  the  vicinity 
of  Mach  one.  The  bow  wave  passes  behind  the  static  ports  at 
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=  1.03  or  so,  and  the  pressure  error  becomes  quite  small 
(perhaps  zero) . 

In  summary,  the  form  in  which  data  from  airspeed  cali¬ 
brations  is  put  depends  on  the  type  of  airspeed  system,  speed 
range,  and  importance  of  Kach  nusT^ber  and  weight  effects. 


Calculation  of  Calibrated  Airspeed  and  Calibrated  Altitude 

The  three  cases  noted  belov;  are  generally  used.  Equations 
for  making  corrections  for  position  error  to  both  airspeed  and 
altitude  are  presented  for  these  three  cases  assuming  that  no 
error  exists  in  the  total  pressure  source.  With  this  assump¬ 
tion  a  common  value  of  AP^  is  applied  to  airspeed  and  to 
altitude. 

Case  I; 


APp 


(46) 


nc 


or  its  alternate 


Sc 


(47) 


is  readily  found  at  the  test  conditions  and  AP  from 


p'  -*0 


1C 


P 


AP. 


■^qc. 

Cf.  ‘c 


(48) 
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where  q  is  cowputed  from  the  modified  form  of  equation  (9) 
^ic 

1c:„  -  II  (49) 


Case  II:  AM  =  f(M.  ) 
pc  1C 

By  definition 


which  can  be  v/ritten  as 


1-^ 


From  equation  (20)  for  subsonic  flight 


1  - 

(1+0.2m2)^*® 


Expansion  of  equation  (52)  in  a  Taylor  series  about 
through  the  first  tv/o  terms  produces 


APp  1.4MipAM  0.7(1- 1.6Mip^)AMpp=^ 


1  +0.2M;, 


(l+0.2Mij,)^ 


From  the  machmeter  equation  for  supersonic  flight 


i4;i;q9im.  ‘ 


(7M:^^-1)^*^ 
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and 


Pt' 


166.921 
(7M^-  1)2*5 


(55) 


Following  the  same  procedure  that  v;as  used  with  the  subsonic 
equations 


'^Pp  .  7(2Mi^-’-l)4Mp^ 


7  ( 21  Mjc'^  -  23.5  Mj/  +  4)AMg^ 


‘1C 


Mic"(7Mi,2_i)2 


(56) 


Rev;riting  equation  (6)  as 


(57) 


and  dividing  equation  (53)  by  equation  (57) 


l-4MieAMp^  +  0-7(1  ~l-6Mic^)AMpc^ 
(W0.2M|c^)  (1  +  0.2Mic^)^ 

[(1  +  0.2 -  1] 


(58) 


for  <  1.0.  The  equivalent  expression  for  the  supersonic 
case  is  obtained  by  dividing  equation  (56)  by 


’’Cjc  _  166.921  MjJ 

*^3  (7M;/-1)^’6 


v/hich  results  from  equation  (8)  ,  to  give 


(59) 
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7(2Mie^-l)AMp^  7(21M|p'*-23.5Mip^.t  4)AMpp^ 

Mip(7Mip^-l) _ Mip^7Mi/-l)^ 

,  166.921  _  1, 

2  2.5 

(7Mjp2-l) 


As  for  Case  I,  AP  is  determined  from 

P 


APn 

he 


requiring  equations  (59)  and  (49)  for  subsonic  flight.  At 
supersonic  speeds  equation  (6C)  and  the  modified  form  of 
equation  (10) 


166.921  (-5^)'^ 

; _ a 

V-  2  2.5 

[7(Acf  -  11 

a 


are  used. 


Case  III:  =  f(V.  ) 

pc  1C 


From  the  relationship 


APp  .  (Pt'  -  P3)  -  (P,'  -  P3) 
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reproducible 

AP  is 
P 


AP. 


q.  -  qc; 


IC 


f(VJ  -  f(V;c) 


Froir,  equations  (9)  and  (49) 


V  2  3.5  y.  2  3.5 

APn  -  Paor  Itl  +  0.2(-^)  ]  ~  [1  +  0.2(-i£.)  ]  ! 


Expanding  ecmation  (65)  in  a  Taylor  series  through  th-c 
tvo  terRs  and  evaluatina  at  V  =  V. ^ 

C  1C 


AP. 


1  4P  V- 

nc  V-  2  2.5 

®SL 


0.7  P 

®SL  V-  2  1*5  V  2  9 

+  - ^[l+0.2(^)]  [l  +  1.2(-i£)  lAV  2 

^SL  ®SL  ®SL  ^ 


or  dividing  by  AV 


AP. 


AV, 


pc 


pc 
1 4P  V 

asLhc 

2 

^SL 


V:  2  2.5 
[l+0.2(-i£.)  ] 

^SL 


0.7  P. 


aSL 


V  2  1*3 


V;.  .2, 


+  - ^[l+0.2(-^)]  [1 +  1.2(-i^)  ]AV 


“SL 


"SL 


a 

®SL 


Follovrinc  the  same  rrocedx'.re  for  supersonic  flight, 
fron  equation  (10) 


AP. 


166.921  (^)^ 

p  I _ "SL 

^SL  V  2  2.5 

[7(-.£-)  ~1] 

®SL 


166.921  (-iiS-)^ 

_ ^SL 

V-  2  2.5 

[7(Jk.)  -1] 

^SL 


(64) 

(65) 

first 

(66) 

(67) 

(69) 
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Expanding  ecuation  (68)  in  a  Taylor  series  through  the  first 
two  terras  and  evaluating  at 


AP„ 


7  X  166.921  Pg,,  2(^)^  -  1 

a-,  ap 


“SL 


{ 


"SL 


'SL 


V*  2  3.5 

[7(J^)  -1] 

^SL 


!AV, 


pc 


7  X  166.921  R 


V-  4  V*  2  9 

[14(^)  -9(-i^)%6]AV  2 

fSL  ( ^SL  ^SL _ 


*SL 


V>  2  ‘1--5 

[7(J1^)  -  1] 

®SL 


(69) 


or  dividing  by  AV 

P"- 


AP„ 


AV. 


pc 


1168.45P  (if 
_ ^SL 

®SL 


V-  2 

2(Jk.)^  1 

^SL _ 

V-  2  3.5 

[7(-^)  -  1] 

®SL 


+  584.224  P. 


^SL 


V  4  V  2 

[14(-i£-)^-  9(-AS-)  +  6]AV_ 

^SL  “51  ^ 


"SL 


"SL 


V*  2  4.5 

[7(^)-'  -  1] 

®SL 


(70) 


AP  for  Case  III  is  found  froir. 
F 


AV. 


AV 


pc 


pc 


(71) 


AP  can  be  determined  as  described  above  for  each  of  the 
P 

three  cases  and  then  used  to  find  both  H  and  V  from  U .  and 

c  C  1C 
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V.  .  To  calculate  ii  .  P  is  first  found  from  the  appropriate 
X  o  c  s  * 

pressure  -  altitude  relationship  (equations  (1)  through  (4)) 
at  Ii.  .  Then  ambient  pressure  is  computed  from 

•L  ^ 

Pa  -  Pa  -  APp  («) 

II  ,  corresponding  to  P  ,  raav  then  be  found  through  the  sam.e 

pressure  -  altitude  relationships. 

V  may  be  found  in  a  similar  fashion  by  computing  q  , 

O  C  t 

1C 

corresponding  to  from  equation  (40)  or  equation  (62) ; 

then  impact  pressure  from 

%  “  %.  +  (73) 

1C 

V  /  corresponding  to  may  be  computed  from  equation  (0)  or 
c  c 

from  equation  (10)  through  an  iterative  method. 

PPXSSURE  LAC  ERROR 

Pressure  gages  such  as  airspeed  indicators  and  altimeters 
are  subject  to  pres.suro  lag  errors  when  airspeed  or  altitude 
are  changing.  Pressure  at  the  source  differs  from  that  regis¬ 
tered  by  a  pressure  measuring  device  because  of; 

(1)  pressure  drops  in  tubing  resulting  from  viscous 
friction  betv/een  the  moving  air  and  the  walls  of  the 
tubing 

(2)  inertia  of  the  air  in  the  tubincr 

NOT  REPRODUCIBLE 
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(3)  pressure  drops  across  orifices  and  restrictions 

(4)  acoustic  lag,  i.e.,  the  time  required  for  a  pressure 
disturbance  to  travel  the  length  of  the  tubing,  and 

(5)  instrument  damping,  inertia,  anci  ''riction 

Since  an  alti’ meter  is  an  absolute  pressure  measuring 
instrument  the  effect  of  lag  on  indicated  altitude  readings 
is  fairly  obvious.  For  example,  during  a  climb  the  indicated- 
altitude  tends  to  be  less  than  the  actual  altitude;  an  air¬ 
speed  indicator,  being  a  differential  pressure  gage  is  affected 
by  lags  in  both  total  and  static  pressures  so  that  the  .error 
may  be  either  positive  or  negative. 

Since  the  error  in  both  pressures  is  in  the  sarue  direc¬ 
tion  the  net  effec  on  impact  pressure  and  hence  calibrated 
airspeed  is  compensating.  Corrections  to  altitude  are,  in 
general,  of  more  consequence  then  corrections,  to  airspeed, 

From  the.  above  list  of  factors  which  affect  the  indicated 
pressure  in  an  airspeed  system  it  is  apparent  that  a  complete 
mathematical  treatment  of’  the- response  to  varying-  pressure 
would  be  prohibitly  complex.  It  has  been  found,  however, 
that  lag  can  be  attributed  largely  to  viscous  friction  and 

be  adequately  described  by  the  equation 

-  P(t)  -  (74) 

by  which  the  indicated  pressure  lags 


that  rhe  system  can 


A-i-  +  R 
dt  * 


where  X.  is  the  time 


behind  the  source  pressure.  Use  of  this  equation  assumes 
that 

(1)  the  mass  of  the  air  in  the  system  is  zero 

"(2)  the  rate  of  chance  of  the  applied  pressure  is  nearly 
constant 

(3:)  laminar  flov  exists  (For  this  to  be  true  it  is  neces¬ 
sary  that  the  Reynolds  number  be  less  than  2000;  in 
typical  airspeed  systems  a  valtie  of  500  is  seldom 
exceeded  in  flight.) 

(4)  the  pressure  lag  is  .small  compared  with  the  applied 
pressure  (This  is  generally  the  case;  hov/ever,  at 
very  high  altitudes  this  assumption  becomes  critical.) 

(5)  the  acoustic  lag  is  negligible  (This  assurption  can  be 
easily  checked  by  computing  acoustic  lag  from  the 
length  of  tubing  and  the  speed  of  sound  in  the  tubing 
and  comparing  it  to  X  from  equation  (74)). 

(6)  the  pressure  drop  across  the  orifices  and  restrictions 
is  negligible  (This  is  true  only  if  a  minimum,  of  such 
restrictions  exist  so  that  the  tubing  is  nearly  a  smooth, 
straight  "pipe”  of  uniform,  diameter.) 

For  a  constant  rate  of  change  of  applied  pressure, 

P(t)  =  (dP/dt)t  and  equation  (74)  may  be  solved  as 


A -yd.  +  p. 

/^f  i 
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From  dimensional  analysis  it  can  be  shown  that,  for  a 
particular  installation, 


A 


(76) 


Hence  if  the  lac  constant  is  obtained  at  one  value  of  p/P 
it  may  be  extrapolated  to  other  values  by  the  expressit)n 


h.  =  ^1^2 
^2 


(77) 


Usually,  experimental  data  are  used  to  Compute  a  sea  level 
lag  constant  from,  v^hich  the  lag  constant  at  any  values  of  y 
and  P  can  be  obtained,  using  the  equation 


X  *• 


(78) 


With  the  lag  constants  for  the  static  and  total  pressure  sys~ 
terns  known,  the  error  in  altimeter  and  airspeed  indicator  read¬ 
ings  due  to  pressure  lag  can  be  calculated  for  any  set  of  test 
conditions . 


Corrections  to  Altitude 

The  lag  constant  for  the  static  pressure  system,  X  ,  can 

s 

be  defined  from  equation  (75)  as 


V7ith  the  approximation  that  dP  /dt  =  dP  /dt  equation  (79) 

s 

becomes 


AP. 


n 


dPg/dt 


(80) 


From  the  assumption  Of  hydrostatic  equilibrium 


dP_  -  -pgdH: 


1C 


(81) 


Then  from  equations  (80)  and  (81)  the  altimeter  lag  correction 
is 


AH 


ICg 


£ic 

dt 


(82) 


Replacing  P_/P-  with  equation  (78)  may  be  written  as 

S  cLrtr  S 


^Hic  1 

f^SL  h 


(83) 


The  altimeter  lag  correction  can  be  evaluated  from  the  experi¬ 
mentally  determined  value  of  /^SL  function  of 

temperature  from  figure  3) ,  5  corresponding  to  H.  .  and 

w  xc 

dH^^/dt  the  rate  of  change  of  observed  altitude  corrected  for 
instrument  error. 
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Corrections  to  Airspeed 

Lag  in  both  the  total  and  static  pressure  systems  may  be 
accounted  for  by 


^C: 


1C( 


P  *  -  P 
no 


(85) 


The  error  in  impact  pressure  due  to  lag  is,  by  definition 


which  may  be  stated  as 


■  ‘Pu-P.’)-(Ps,-Ps) 


It  follows  from  equation  (75)  that 


dR’  dP„ 


^iC£ 


^  dt  ^ 


dt 


Differentiating  equation  (85)  with  respect  to  time 


-ICJi,  . 


■iP^  ■iPs, 


dt 


dt  dt 


Substituting  equation  (89)  in  equation  (88) 


dqc.  dP_ 


dt 


dt 


(86) 


(87) 


(88) 


(89) 


(90) 


f 
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with  the  approximations 


and 


dt  dt 


equation  (90)  becomes 


-  (.3  -  A,)^ 


(91) 


Differentiating  equation  (81) 


dPs 


-pgdH 


1C 


(92) 


and  differentiating  the  subsonic  airspeed  equation  (equation 
(49)) 


Dividing  by  dt 


dq, 


^ic 


I’^Paqr  '^ic  V.  2  2.5 

- !SLif  [1  +  0.2(-i£-)  ]  dVj 

2  ^CT 

®SL 


1C 


(93) 


l-'^PacT  ^ic  V-  2  2.5  dV- 

- JL21[1  +  0.2  (^)] 


(94) 


Differentiating  equation  (88)  and  substituting  together 
with  equations  (89)  and  (90)  in  equation  (87)  the  airspeed 
indicator  correction  expressed  as  a  finite  difference  is 


AV, 


ic^ 


A 


/  ^  dHic 

(Ag  A{)pg 


l-4Paci  V-  2  2.5 

- +  0.2(-^)  ] 


(95) 


*SL 


"SL 


'^ic  1  ^SL 

With  the  supersonic  airspeed  equation,  equation  (59) 


AV 


ICfl 


‘at  V-  6.. .  V; 


(96) 


3738.11(-i^)  [2(-i^)  -  1] 
^SL  ^SL _ 

[7(Ji£)^  -  1]®-^ 

^SL 


Indicated  airspeed  corrected  for  instrument  and  lag  error i 
V.  ,  may  be  found  from  the  above  equations  and 

At 


ICfl 


Vic  + 


(97) 


As  in  the  case  of  the  altimeter,  /Vg^)  (l/<Sg)  ; 

ic 

similarly 


^Hic  aSL 


^SL 


f'sL  Ps  +  q„. 


(98) 


nc 
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SECTION  V 

DETERMINATION  OF 
EXCESS  THRUST 


SUMMARY 

Excess  thrust  has  a  number  of  advantages  to  recommend  it 
as  a  basic  parameter  in  defining  the  performance  of  an  aircraft. 
For  example,  corrections  in  terms  of  excess  thrust,  are  most 
easily  derived  from  equations  of  motion,  and  procedures  common 
to  both  climbs  and  level  accelerations  may  be  set  up  using 
excess  thrust  as  a  basis  for  arriving  at  the  desired  parameters 
in  both  cases.  Various  means  of  computing  excess  thrust  are 
available.  The  advantages  and  disadvantages  of  each  are  dis¬ 
cussed,  and  the  equations  for  computing  it  are  derived  in  detail. 
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SYMBOLS  USED  IN  THIS  SECTION 


Symbol 


Definition 


Units 


components  of  accel¬ 
eration 


ft  per  sec’ 


acceleration  factor 


dimensionless 


acceleration  factor 
associated  with 
energy  height 

roll  angle  about  the 
airspeed  vector 

aCj^/3a 

drao 


dimensionless 


per  rad 


^  •r*  r- 


net  thrust 

local  effective  acnelefa- 
tion  due  to  gravity 

reference  acceleration 
■due  to  gravity 

components  of  accel¬ 
eration  of  gravity 

tapeline  altitude 

pressure  altitude^ 

energy  height 

uhi^:  vectcj  3 

thrust  angle  pf  inci.*? 
dence 

moment  of  inertia  about 
y  axis 


'f  t  per  sec" 


ft  per  sec’ 


ft  per  sec’ 


slug.rft® 


dimensionless 


J2  fda  r <^4 


coefficients  of  the 
zonal  harmonics  of  the 
earth's  gravitational 
potential 

distance  from  vane  to 
aircraft  eg 


dimensionless 


T 

lift 

lb 

M 

flight  Mach  number 

dimensionless 

components  of  lead 
factor 

dimensionless 

P 

roll  rate 

rad  per  sec 

ambient  pressure 

lb  per  ft^ 

q 

pitch  rate 

rad  per  sec 

q 

dynamic  px  S't^nure 

lb  per  ft^ 

r 

local  radius  cf  Ihe. 
earth 

ft 

r 

ya^'j  rate 

rad  per  sec 

equatorial  radius  of 
the  earth 

ft 

s 

Laplace  operator 

_  -  - 

characteristic  area  of 
vane 

t 

time 

sec 

hprizoiital  compbnent  of 
adrcrd^t  speed  relative 
to-  -^e  -ground 

ft  per  sec 

inertial  speed 

ft  per  sec 

true  airspeed 

ft  per  sec 

wind  speed 

ft  per  sec 

$ 

edrcraft  velocity 

Induced  -by  the  earth ' s 
rotation 

ft  per  sec 

w 

airplane  gross  weight 

lb 

x,y,z 

Cartesian  coordinate 
system  (s'lbscripts 
denote  particular  axes 
system) 

a 

angle  of  attack 

rad 
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I 


^%com  bending 

change  in  sensed  angle 
of  attack  due  to  boom 
bending 

rad 

change  in  sensed  angle 
of  attack  due  to  pitch 
rate 

rad 

change  ih  sensed  angle* 
of  attack  due  to  upwash 

rad 

Aa 

(0 

change  in  sensed  angle 
of  attack  due  to  dynamic 
lag 

rad 

Y 

flightpath  climb  angle 
measured  from  the  geo¬ 
centric  horj.2ontal 
plane 

rad 

aircraft  geodetic  lati- 
__tude 

rad 

-  '■ 

aircraft  geocentric  lat¬ 
itude 

rad 

misaiinement,  angle 

rad 

? 

damping  ratio 

dimensipnless 

AX^ 

yL. 

difference  petweeh  -^e 
"  aircraft  longitude  and 
th'e  longitude  ;pf  the 
iadar  coordinate  origin 

xad 

% 

product  of  the  universal 
gravitational  constant 

mass  of  tdie  .earth 

f .per  sec^ 

a 

f 3.i<'htpa^  heading  angle 
measured'  from  true  north 

rad 

vehicle  ground- track 
heading  angle  measured 
from  true  north 

rad 

bank  angle 

rad' 

<f>r 

angle  between  radar 
X^-akis  and  true  north 

rad 

direction  from  which  the 
wind  blows  (from  true 
north) 

V-8 

(0 


frftquency 


rad  per  sec 
rad  per  sec 


undamped  natural 
frequency 

rotational  of.  rad  per  see 

wind- axes  sys®"2rs 

Wq  angular  velocity  of  the  rad  per  sec 

earth 


\ 


Subscripts 


b 

eg 

e 

g 

ic 

m 

6 

■p:  ■ 

y 

'll! 


body  axes 

center  of  gravity 

geocentric  aices 

local-geocentric  axes 

indicated  corrected  tor 
instrument'  b'rrdr 

.measured 

maximum  value 

radar  axes 

vane 

wind  axes 
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GENERAL  DISCUSSION  OF  TUE  VARIOUS  xMETIiODS 


Airplane  performance  data  can  be  determined  by  several 
different  methods  requiring  different  instrumenraticn. 

These  methods,  which  have  received  varying  degrees  of  use, 
are; 

a.  Airspeed-altitude 

b.  Accelerometer 

c.  Position  measurement 

d.  Time  history  of  dynamic  pressure 

e.  Rate  of  climb  indicator  eohiiec ted  to  total  pressure 

f.  Air  temperature  thermometer 

Methods  d,.  e^,.  and  :f  are’-dascribadiiiT  reference  1 -.aiid-  rn 
Chapter  7  of  reference  2  b.ut.  -are  excluded  from  further  discus-  . 
sion  since  they  have  never  received  popular  acceptance .  ,0f  the 

-remaining  Ihtesir  the  airspeed-radtitude  metJiod  Has  been  widely 
=used  for  rihany  y,ears.>;=  accelercmeters  have  ,beeir  used  during: 
numerous  flight  test  programs;  position  measurement  (e.g., 
radar)  has  been  employed  in  .some  isolated  instances i 

The  most  -convenient  parameter  with  which  to  work  in 
standardizing  airplane  performance  data  is  excess  thrust. 
Corrections,  in  terms  of  excess  thrust,  are  most  easily 
derived  from  equations  of  motion.  Excess  thrust  can  be 
related  easily  to  rate  of  climb,  turning  performance,  and 
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other  performance  parameters.  Standardization  procedures 
common  to  both  climbs  and  level  accelerations  may  be  set 
up  using  test  excess  thrust  as  a  basis  for  arriving  at  the 
desired  standard  parameters  in  both  cases.  Equations 
defining  test  excess  thrust  are  derived  in  subsequent  para¬ 
graphs  for  the  three  methods,  a,  b,  and  c,  starting  from 
the  same  basic  equations  of  motion. 

Since  the  quality  of  standarized  data  can  be  expected 
to  be  no  better  than  that  of  the  test  data',  some  observa¬ 
tions  in  regard  to  the  calculation  of  excess  thrust  are  set 
down  for  the  three  methods.  Advantages  and  disadvantages 
of  each  are  considered. 


AI'RS^■SD:?ALTiTUpEf'^^ETHO.p 

The-  bulk  of  tlie  peffprinancS  test  programs  which  have 

been  conducted  up  to- the  present  time  have  made  use  of 

measurements  cif  airspeed,  altitude,  and  time  {.usually  an 

airspeed  -iidi Ciator,  and  altimeter,  and  a  clock  mounted  on 

a  £Jho tor  panel)  tc  gather  performance  data.  This  '.^thcid 

! 

has  been  in  roU”c.;.ne  use  since  the  early  1950 's.  Several 
schemes  for  processing  the  data  have  been  tried.  The  one 
in  most  comiflon  use  in  later  years  for  both  climbs  and  accel¬ 
erations  has  been  to  compute  total  energy  as  a  function  of 
time,  then  curve  fit  and  differentiate  to  find  rate  of  change 
of  energy  (or  excess  thrust) . 
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The  chief  merits  of  the  airspeed-altitude  method  are  (1) 
the  instruments  are  alv.’ays  available  since  they  are  needed 
on  board  any  test  aircraft  as  a  part  of  its  test  instrumenta¬ 
tion  system.  Hence  "spec  al"  instrumentation  {■  iz  sensitive 
accelero  eters)  or  ground  based  positioning  equipment  is  not 
needed,  and  (2)  the  instruments  are  very  reliable  and  usually 
quite  consistent  in  behavior  so  that  relatively  little  data 
compared  to  sensitive  accelerometer  system.s  should  be  lost 
because  of  instrument  malfunction.  Their  reliability-  is  also 
generally  superior  to  that  of  position  measuring  equipment 
since  contact  v/ith  an  aircraft  is  occasionally  lost  due  to 
excessive  range,  clouds,  etc. 

The  airspeed-altitude  method  of  calculating  excess  thrust 
has  been  found  to  be  the  le'ast  accurate  of  the  methods 
described.  (See  reference  3v  example)  .  There,  .are  several 
factors  which  contribute  to  the  inaccuracy  of  this  method,  the 
most  important  being  the  errors  in  determining  true  airspeed 
and  tapeline  rate  of  climb  due  to  errors  in  ambient  tempera¬ 
ture,  pressure  lag  during  climbs,  readability  of  instruments, 
etc.  The  reason  these  errors  are  so  important  is  because 
excess  thrust  is  computed  by  differentiating  specific  energy, 
E/W  (energy  height,  H^) 


R. 

w 


S3 


H 


+ 


2g 
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and  any  errors  in  these  parameters  are  amplified  by  the 
differentiation  process.  Also,  the  magnitude  of  che  time 
interval  used  for  curve  fitting  and  differentiating  to 
obtain  rates  of  change  of  specific  energy  (and  consequently 
excess  thrust)  may  have  a  decided  bearing  on  the  results. 

If  small  intervals  are  used,  spurious  variations-  are  intro¬ 
duced  into  the  calculation  of  excess  thrust.  On  the  other 

/ 

hand,  when  excessively  large  intervals  are  used,  variations 
which  actually  existed  are  eliminated,  or  at  least  reduced 
in  size. 

Another  disadvantage  of  the  airspeed-altitude  method  is 
that  corrections  for  wind  gradients  are  needed  which  are 
frequently  quite  sxabstantial  for  climbs.  In  contrast  only 
relatively  minor  corrections  for  variations  in  normal  load 
factor  are  required  to  standardize  climb  performance  when 
using  any  of  the  other  methods. 

As  flight  speeds  and  altitudes  increase,  the  accuracy 
inherent  in  the  airspeed-altitude  method  becomes  worse. 
Above  a  Mach  number  of,  say,  two,  other  means  of  obtaining 
excess  thrust  should  be  used,  if  possible. 

ACCELEROMETER  METHOD 

Flight  tests  to  evaluate  the  practicability  of  measur¬ 
ing  longitudinal  accelerations  by  means  of  a  sensitive 
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accelerometer  were  conducted  in  the  mid-1950's.  (See 
reference  4.)  Sensitive  accelerometers  were  installed  in 
category  test  aircraft  for  use  as  a  prime  source  of  per¬ 
formance  data  starting  in  the  early  1960 ‘,s.  Accuracy  of 
the  data  recorded  during  these  days  was  comparable  to  that 
acquired  with  airspeed-altimeter  measurements;  however, 
with  improvements  in  accelerometers,  installation  design, 
and  means  of  recording,  data  from  accelerometers  are  at 
present,  decidedly  more  accurate  than  that  from  the  air-r? 
speed-altitude  method. 


Two  different  types  of  installation  have,  been  tried. 

In  the  first,  a  sensitive  accelerometer  has  been  mounted 
oh  a  vane  (similar  to  that  used  to  measure  ingle  of  attack) 
so  that  the  sensitive  axis  of  the  accelerometer  remains 
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ali'ned  with  the  Ideal  airflow.  Unfortunately,  the  upwash 
angle  caused  by  the  presence  of  the  pitot  head  and  perhaps 
the  nose  of  the  aircraft  is  significant,  and  accurate  cal¬ 
ibrations  must  be  made  in  order  for  the  desired  component 
of  acceleration  along  the  airplane's  velocity  vector  to 
be  found.  Secondary  effects  which  may  be  of  consequence 
if  the  pitch  angle  is  changing  rapidly  are  due  to  pitch 
rate,  pitch  accelerations,  and  dynamic  effects  causing 
a  lag  in  the  position  of  the  vane.  These  effects  can  be 
expected  to  be  negligible  during  level  accelerations  and 
pirobabiy  during  climbs  but  may  require  minor  corrections 
to  measured  accelerations  during  roller  coaster  maneuvers 
or  qther  tests  when  pitch  rates  are  high.  In  the  second 
type  of  installation,  a  sensitive  accelerometer  is  hardr 
mounted  at  or  near  the  airplane's  center  of  gravity.  In 
this  case  the  location  of  the  sensitive  axis  -of  the  accel¬ 
erometer  must  be  very  well  kjnown.  Further,  the  accelera- 
t-iqn^ririust  -be  -resolved  through^  the  angle  of  attack  and  cor¬ 
rections  have  to  be  made  to  it  similar  to  those  described 
fqr  the  vane  mounted  accelerometer. 

The  most  attractive  features  of  the  accelerometer 
method  are  that  load  factor,  n  ,  along  the  flight  path  is 
obtained,  after  making  the  corrections  and  axis  transfor¬ 
mation  described  above.  The  resultant  force  is  the  engine 
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thrust  minus  the  airplane,  drag-  neither,  of  'which  are 

affected  greatly  by  wind  gradients.,  >/^lso;:,  ihstantaneous 
values  are  recorded ,  and  uncertainti^es  incurred  by-  data 
smoothing  and  differentiation  required  , by  the  other  methods 
are  avoided. 

The  reliability  of  the  flightpath  accelerometer  is 
better  than  that  of  position  measuring  equipment  (radar 
and  Askania  camera)  since  to'  iaccelerometer  is  carried 
with  the  aircraft  and  there  is  no  danger  of  incomplete 
flight  coverage;  however,  experience  has  shown  the  flight- 
path  accelerometer  to  be  less  reliable  than  the  airspeed- 
altitude  instruments. 

The-  most  recent  accelerbmeter  installations  have  been 
a  two-axis  accelerometer  systemrwhich  measures  normal 
load  factor  as  wdll  as  load  factor  along  the  flight  path; 
This  makes  for  a  consider^le  improvement  in  the  accuracy 
of  the  resultant  acceleration;  along  the;  flight,  path 
together  with  an  improved  normal  load  factor  for  use  in 
the  standarizatiqn  equations. 

In  addition  to  one-  and.  two-axis  .accelerometer  is.ysteKS , 
three-axis  accelerometer  data  from  an  inertial  navigation 
system  have  been  used  with  marked  success  as  a  prime  source 
of  performance  information.  Data  from  a  threer-axis  system 
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mounted’  on  a  stable,  platform  has  the  considerable  advantage 
that  neither  upWash  angle  nor  angle  of  attack  need  he 
known  if;  orientation  angles  to  locate  the  aeceleratioh 
components  are^  known;  they  are  generally  available  from  pn- 
bp^d;  computers  in  aircraft  which  have  inertial  navigation 
systems  as  standard  equipment.  Drawbacks  to  the  use  of; 
these  systiemsi  are:  additional  complexity  with  attendant 
maintenance  requirements#  emd  need  for  pre-flight  and  post¬ 
flight  checks.  Also#  operable  navigation  systems  have  not 
been,  installed#  in  general#  in  aircraft  used  for  category 
performance  testing. 

POSITION  IffiASUREMENT  METHOD 

Radar  eind  Askania  camera  data  have  been  .used  to  com¬ 
pute  perfomahce  information  in  only  isolated  instances, 
although  both  have  yielded  satisfactory  results  (reference 
3hi  since  the  data  from  both  sources  are  used  in  the 
spie  way#  extrept  for  processing,  of  raw'  ;data  to.  tind  posi¬ 
tion  and  velocity  components:,  they  are  discussed  together. 

The  accuracy  in  both  cases  depends  primarily  on  the 
quality  of  the  tracking  data,  which  in  turn  depends  on 
such  factors  sas  nixmber  of  recording  stations,  range,  ele¬ 
vation  angle,  etc.  As  in  the  airspeedf-altitude  method, 
any  errors  in  tracking  data  are  ^plified  in  the  differ¬ 
entiation;  process  required  to  compute  excess  thrust. 


The  accuracy  inherent  in  these  ;systeriis>  however,  will 
produce  excess  thrust  ot  higher  quality  than  will  the 
airs  peed- altitude  method.  A -serious  disadvantage  may  be 
that  coverage  is-  not  always  available/  especially  ifpr 
high  .performance  aircraft-  Also,^  data  turnaround  time' 
may  be  excessive. 

Like  the  accelerometer  method,  phiy  slight  correcr 
tions  for  wind  gradients  (for  normal,  load',  factor)?  are 
heeded . 

Table  i  shows  a  rating  of  each  method  from  the  stand¬ 
point.  of  accuracy,  reliability,  aircraft  equipment  requirec' 
(the  least  ainpunt  being  considered  best)-,,  and  data  pro¬ 
cessing  effo.rt;  (least  required  by  engineering,  personnel 
considered  best) . 


Table  1  -  Comparison  of  Methods  for  Caiculatihg, 

Excess  Thrust' 


Rating 

Accuracy 

Reliability  ' 

.  .Aircr^t  Equipment 
i  „  Required  . 

'Data:Process5ng 

Efiort 

"  > 

;  1 

Plight  Path 
Accelerometer 

Airspeed- 

Altitude 

Askania 

iRadar  ; 

2 

Askania 

Plight-'Path 

Accelerometer 

:Radar 

/ 

Plight  Path  = 

AccelerometVr  < 

3 

Radar 

Radar 

Plight^Path 

Accelerometer 

' 

Askania 

4. 

Airspeed-  > 
Altitude 

Askania 

Airspeed- 

Altitucie 

‘  Airspeed-  ' 

«  Altitude 
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EQUATIONS  .OF  MOTION 

Components,  of  the  earth's  gravi tational.^at traction  taken 


from  the  section/ 


Prbpertielt ,  are 


[sJgslnfiL 

*0 

“  1^4  <  7^  A-3  ~  ^sin?  Sl  )  sin  5  J-Cos  8^ 


0 

+  2  J3  ( A  3  -  5  sin2  Sl  )  sin  Sl 
~  i  >^<  3  -  30  sin2  +  35  sin^ 


in  a  north-/  gepcentrically-directed  system.  Because  of 
the  earth^s  equatorial  bulge/  is  always  pointed  toward 

. ,  ■  g  ; 

the  equator  (i.e,/  toward^  the  soutn  in*  the  northern  hemi¬ 
sphere  and  toward  the  north  in  the  southern  hemisphere) . 
g  is  assumed  to  be  zero  and  g  is  directed  toward  the 

'  g 

earth's  center.  By  transforming  these  components  of 
gravitational  .attraction'  through  the  angles  0 ,  y,  and  B/ 
the  resultant  along  the  velocity  vector  is  found.  The 
longitudinal  equation  of  motion  in  the  wind-axes  system 


L  ,+  P„  8in(jtt  +  iij, ).—  ^  (  KgCOS  y  +  g^sitnycos aj 


(7) 


It  is  convenient  to  cqitibine  gravitatibrial,  attraction 
and-  aircraft  acceiieration  by  defining  load-  fabtpfs  along 
^e;  airspeed  vector  and  along  the  lift  vector.  Acceler-r 
ometers  sense  these  ;ioad  facitors  directly  .>  Along'  the 
airspeed  vector  the'  load  factor  is 


n-  '«  g.cosycosc  +  a_  i 

*w  Sf  " ’S;  *g  "  '  *w 


(8) 


In  the  vertical  direction  along  the  lift  vector,  the 
load  factor  is 


“  -LCgL  Coiy  +  g- sin  y cose- a„  ) 
®r  ■■  w 


(9) 


Equations  (6)  and  (7)  iiiay  then  be  rewritten  as 


5'nC0s(a,+  lT)-D  - 


(10) 


and 


Li  '+  sin  ( (z  4*-  i*p  ^  **  Hg  yi 

w 


(11) 
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AIRSPEfcD-ALTITUDi:.  METHOD 

Precise  derivations  of  equations  are  made  which  may  be 
used  ;with  on-boafd  measurements  of  airspeed  and  altitude 
to  give  excess  thrust.  The  magnitude  of  some  of  the  terms: 
is.  quite  small,  and>  considering  the  inaccuracies  in  con- 
ventiphai  instruments.,  they  may  (be  safely  eliminated.  If,, 
howevefv  improved  instrumentation  becomes  available,  use 
of  the.  precise  equations  may  be  found  desirable. 

DERIVATION  OF  GENERAL  EQUATIONS 
Inertial.  Velocity  in  Wind-Axes  System 

Referring  to  figure  4  in  the  section.  Coordinate  Sys¬ 
tems  and  Transformations,  the  inertial  velocity  may  be 
expressed*  as 

+  (12) 


In  the  wind-axes  system  the  airspeed,  V^,  is  V^i.  Per¬ 
forming  transformations,  first  through  the  angle  a,  then 
through  y  and  B,  the  velocity  due  to  the  earth’s  rotation 


where 
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(21) 


(22) 


The  acceleration  becomes 


iSi 


I  I.  - 


which  may  be  found  by  differentiating  equation  (18) : 


m  \^-\^cosycoti^~o)  +  V,,(  y8tnycos(0-<7)  +  (0  -  (7)cosy8iB(^^' -  a)] 

+  OiJje  +  UcosS^^cosysina  +  ajir  +  hX-  S^sinS^^cosysina  —  ycosS^sinyslnff 
+  acosS^cosYcosa ) 


(23) 


-  V^^[s'"y8iiiBcos(0  -  a)' +  cosBsin(0- a)]  -  V^[ycosyslnBcos(0  -  a) 

•  •  * 

+.  B  ainyco3Bcos(0  -  <r)  -  (0- a)8lny8inB-8in(0  -  ct)  -  B  sinBsin(0  -  a) 

^  • 

+  (0-O')co8Bcos(0-<t)]  +  6)^(r-i- h)coaS|^(siny8ina8inB  +  cosffco8B 

•  # 

—  +  h)S^8iiid£^(8mysiiK7sinB  +  cosacosB)  +  h)  co8S£^(ycosy8inff  sinB 

•  *  •  * 

■f  asitiyooscrsinB  -f  B  sinysinacosB  -  asinacosB  -  BcosasinB) 


(24) 


Yg  “  —  V^^[rinyco8Boos(0  -  ff)  -  sioBsiii(0-(7)l  -  V.^[yco8ycosBcos(0  -  a) 

I 

•  ^9 

-  B'5lny8inBco8(0  -  ff)  -  (0— a)slnyfcosBsiB(0  -  (?)  -  BcosB  8lii(0  -  (?) 

-  (0  -  (7)8iaB  coe(.^  -  (?)]  +  <i»^(r  +  k)  cosS^Csinysina cosB  —  cosc? sinB) 

-  6)Jit  +  h)5j  siiiS£^(8iny8in(7cosB  —  coscrsinB)  +  WcosSj^CycosysinacosB 

•  •  •  • 

+  asiayco8<?(»8B - B8iny8in<7sinB  +  aslnasinB  -  Bcos(7CosB) 


(25) 


If  no  rate  of  roll  exists,*  as  may  be  expected  during  per¬ 
formance  flight  testing,  B,  may  be  set  equal  to  zero  and 
equations  (24)  and  (25)  reduced  to  the  following; 

Vy  n  -  V^JsinysinBcosC^  -  +  cosB  sinC^i  -  a))  -  V^^tycosysinBcosC^  -  a) 

•  #  • 

-  (0  -  CT)sin)'sinB  sin(^4  -  o)  +  {<l>  -  a)eo8Bcosi.\{i  ~  a)]  +  (u^(r  +  b)cosSj^(siny8inasinB 

•  • 

+  cosocosB)  -  &j^{t  +  h)Sj^sinSj^(sinysliiff8inB  +  cosctcosB)  +  +  b)coa5|^(>'co8>'8ino'sinB 

+  ff s;-:yros(7slnB  -  <;sinacosB) 

(26) 

•  •  « 

Vj,  (sin}'cosBcos((/(  -  o)  -  sinBsin(^&  -  <t)1  -  V^[ycosyco8Bco8(^i' -  tr) 

«  •  # 

-  (0  -  ff)sinycosBsin{0  -  a)  -  (^  -  <7)sinBcos(0  -  a)  J  +  cjJp  +  h)cosS£^(8iny8iri(7cosB 

♦  • 

-  vosasinB)  -  <j^(r  f  h)Sj^sinS|^(siny8in(7cosB  ~  cosffslnB)  +  &)J[r  +  h)cosSj^(ycosysiiiacosB 

+  asinycosacosB  +  asinvsinB)  (27) 

Next,  referring  bac)c  to  equation  (20)  ,  it  is  necessary  to 
define  the  rotational  velocity  of  the  wind-axes  system. 

This  is  first  determined  in  the  north,  east,  down  system, 
denoted  by  the  subscript  g; 

bi^,  n  +  A)cos5l  —  ysincrj  ig  +  (ycosa  —  jg  +  [(7  -  (w^+ A)8in5L]bg  (28) 


*B  has  no  influence  on  a^  ,  a  ,  or  a^  .  This  can  be 

shown  by  retaining  terms  which  contain  B  in  equations  (23) 
and  (25)  as  well  as  in  subsequent  equations  which  define 
w  .  All  terms  containing  B  then  vanish  due  to  cancella¬ 
tion  of  terms. 
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wi^  the  rotational  velocity  due  to  a  rate  of  roll  set 


egu^l  to  zero. 


Transforming,  as  before,  to  the  wind-axes  system 


<UWx 

“  •  •  * 

(£t)g+ A)cos£^-  ysinff 

Wy 

Im] 

•  ; 

ycosa  -  5. 

o  -  {<0^  + A)  sinSj^ 

(29) 


Performing  the  matrix  multiplication 

61^  “  [(6)^  AjcosSj^  —  y8m<T]cosyco8a+ [ycosa  —  Sj^lcosyslnff 

-  [c\-  (<u^+ A)  8iiiS£]8iiiy 


(30) 


0)  t  [(<u  +  A)co^,  —  ysinalisinycosasinB  —  8inaco8B)  +  {ycooa  —  S.  )(siny8inasinB 

y  ^ 

+  cos(tco8B)\^+  [o-  -  (<i;^+  A)  siii5^]co8y8inB  (31) 

B  [(<u^+  A)cosSl  —  ysinaKcosCTsinycosB  +  sioasinB)  +  (yco8(T  —  5j^)(sinysin(7cosB 
Z 

-  C08<7sinB)  +  [a  —  ‘(i>.  + A)8iii5i ]co8yco8B 

•  ^  (32) 

simplifying  equations  (30)  ,  (31)  and  (32)  we  have 


^)'(«>85LC08yco8<7  +  slnS|^8lny) -5£^co8y8ina -CTsiny  (33) 
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•  • 

(o)^  +  A)'(cosS^sinyco8(7ainB  —  cosSi^aincrcosB  —  sinS^cosVsinB)  +  ycosB 

-  S£^(sinysiQasinB  +  cosacosB) -forcosysinB  (34) 


! 


“  ((i>^ -I- A)(cos£^cosasinyco8B +cosS£^«ina8inB  —  sinSjjCosycosB)  —  ysinB 

•  . 

_  Sj^CsinyainacoaB  —  cosoainB)  4-(7co8yco8B  (35) 


Again  referring  back  to  equation  (20) ,  to  evaluate  the 
acceleration  due  to  the  rotation  of  the  wind- axes  system 
we  hav,’ 

S^xVj  - 


i 


j 


<u„  <a„  6}^ 

'^x  '^y  ”2 

V  V  V 


(36) 


Components  of  acceleration  are 


<a. 


(37) 


(38) 


a 


(39) 
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Substituting  terus  from  equations  (18),  (30),  (31),  and 
(32)  we  have 


*=  <“6^=^*+  rt)(r  +  hXcosSj^ainSj^cosycosa  -  coJ^S^sinV)  -VJw^+  A)[co8ysl«l(^^i  -  a)sia5£^ 

—  sinycos5|^aiii\2>]  +  y<oJit  +  UcosS^^sin^'sina  —  oujj  +  h)cosS^cosyco8a 

-  y V^8jByoo8(0 -  ff)  +  SV,^sSnycos0  +  V^^ffcosy8lii(v4- :r)  (40) 

•<Uy  "  A)(co8S£^cosff8inycoaB  +  co8Sj^eiii<78inB  —  8iiiS^co8ycosB)  —  V^XsinB 

•  ^ 

~  V|Sj^(8liiy8m(7co8B  —  co8a8inB)  +  V^ffcosycosB  +  A)[co8(^&  -  a)8iii5£^co8B 

-  coaVcoaS^^ain^sinB  —  8iiiy8in(0  —  a)8inS£^sinB]  +  V^yco8ycos(^4  -  a)  sinB 

-  V^.5LC08yc0B\&8liiB  -  V^ff[co8(0  -  a')cosB  —  siny 8li»((4  -  a)  sinB] 

f  2 

+  A)(r  +  h)  [- sina  cosS^sinSj^cosB  +  cosycos  5|^8lnB  +  8inycosacoeSj^slnSL8uiB] 

-  +  h)co8y  8iiiaoo8dj^8lnB  4  ^  hHsinacosd^^cosB  —  8iiiycoe(7co8S^8iiiB)(41) 


^-HcosSLstaycosasIiiB  -  coeSLsJnaoosB  -  sinSj  co8y8liiB)  +  ycosB 
-  5L(8lnysiii<r8inB  +  oosacjsB)  +  CTco8y8ijiB]  —  VJa>^+  A) [cosysiin^ cobSlCosB 
+  col^<l/  -  a)  r  Indj^sinB  +  siiiy8iii(^  -a)8ta5j^co8  B]4  +  A){r  +  h)co85|^(cosycoe5j^co8B 

+  sinasinS^^sinB  +  alnycosvsinS^cosB)  —  \^S£^co8yco8^o[wB  +  \^ff[co8(^  —  0')»'inB 
+  8inysin(0-er)co8B]  +  V^ycosy  coa((^  -  a)oo8B  -  yaij[r  •)■  hkosysina  cosSj^cosB 
”  +  Wcosdj^CsInffsinB  +  sliiyoos«Tco8B) 


(42) 


To  Obtain  total  accelerations  in  the  I,  J,  and  ic,  direc¬ 
tions,  as  indicated  by  equation  (20)  ,  we  combine  terms  as 
follows : 
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(43) 
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Equations  (25)  and  (42)  yield 


a  «  -V.yeosB  +  V.SiCsinysInbslBB  +  cosffcosB) —  V|0cos>'siiiB  -  V|(f;;  +AKBinyoo8ffco85,8inB 

—  siaffcosSj^coaBx-  cosyainSjLsinB)  —  V^^lsinycosC^  -  0)casB  -  8ln(0  -  0'  aiiiB) 

^  V^^.^ii3lny8ia(;4 -0)oosB  4-  cO!rf0-0)s1iiBj  -  V^^(£i»^+ A)[co8y  cosS^cosB 
4-  fct»(0  ~  0)8inSj^8inB  +  sinysiuC^  —  a)8inS£^cosB]  +  +  h)  co8S£^(8lny8in0cosB 

-  coi^ainB)  -  w^{r  +  h)S£^8in5|^{sin>'sin0co8B  —  co808lnB)  AKr  +  WcobSj^* 

•(cosycoaSj^cosB  +  sin08inS|^sinB  4-  8inyco808inSj^co8B)  —  V^Sj^cosycos^cosB  (48) 


The  lateral  acceleration,  a  ,  will  tend  to  produce  a 
sideslip,  the  magnitude  of  which  will  depend  on  the 
aerodynamic  characteristics  of  the  airplane.  The  resulting 
increase  in  drag,  even  in  extreme  cases,  will  be  slight, 
and  the  effect  of  a  on  aircraft  performance  can  be 
safely  ignored. 

Equations  (46)  and  (48)  may  be  modified  to  more  easily 
compute  accelerations  in  the  x^  and  directions  based  on 
on-board  measurements  of  airspeed  and  altitude.  In  doing 
this  the  parameters  V  ,  4,  and  X  will  be  replaced  by 

W 

more  basic  quantities. 

In  the  preceding  derivations  accelerations  have  been 
expressed  in  terms  of  parameters  in  a  geocentric  reference 
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systeni.  The  flightpath  clindj  angle  in  a  geodetic  reference 
system,  however,  is  desired.  Prom  equation  (9)  of  the  sec¬ 
tion,  Geophysical  Properties, 

K  “  Sd  -  0.19323889 sIn2  3D  (49) 

with  angles  expressed  in  degrees.  Equation  (49)  can  then 
be  used  to  relate  the  geocentric  and  geodetic  flightpath 
climb  angles  by  an  extension  of  the  relation  illustrated 
in  figure  2. 


DOWN 


Figure  2  GEOCENTRIC  FLIGHTPATH  ANGLE  FOR  GEODETIC 
HORIZONTAL  FLIGHT 
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y 


Xd  +  Sd)cosct 


(50) 


or  substituting  equation 


(40) 


y  m  yjj  —  0.19323889  sln2  3jj  cos  a 


(51) 


The  inverse  equation  for  geodetic  climb  angle  is 


yjj  ■  y  +  0.19323889  sln2Sj3  cos  a  (52) 

where  the  geodetic  latitude,  can  be  replaced  by  the 
geocentric  latitude,  6^/  with  little  loss  in  accuracy. 

From  the  definitions  of  geocentric  and  geodetic 
climb  angles  the  following  equations  for  rate  of  climb 
are 


h  -  VjSinyu 


(53) 


and 


9  0 

h  +  r  »  Vjsiny 


Next,  the  ground  speed  may  be  expressed  as 


(54) 
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(VjCOSyCOSff  —  V^C08^)l 


+  ( V^cosyslna-' V^sin0)jg 


(55) 


The  heading  cingle  of  an  airplane's  gro\md  track, 


is  then 


a 


V,  cos  y  sin  a  —  V„  sin  ti¬ 
tan  rf - ^ ^ ^ 

Vj.  cos  y  cos  o  —  cos  l/l 


(56) 


which  may  be  used  to  evaluate  the  rates  of  chemge  of 
latitude  and  longitude. 


Vgcoso 
r  +  h 


(57) 


and 


A 


(r  +  h)cosgL 


(58) 


In  the  preceding  derivations  leading  to  a  and  a 

A  z 

w  w 

(equations  (46)  and  (48))  winds  have  been  specified  in  a 


V-34 


geocentric  reference  system.  As  an  aid  in  the  following 
development  it  will  be  assumed  that  the  wind  can  be  con¬ 
sidered  either  geocentrically  or  geodetically  horizontal 
without  error =  The  wind  speed  and  wind  direction  deriva¬ 
tives  may  then  be  expressed  through  chain  differentiation 
as 


’w 


dh 


dVjj, 

dh 


V^sinyu 


(59) 


and 


(60) 


The  geodetic  climb  angle  can  be  replaced  by  a  geocentric 
climb  angle  in  keeping  with  the  assumption  noted  above. 

Introducing  these  cjxpressions  into  equations  (46)  and 
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(48)  we  obtain  directly 


'w 


•  d  Vu,  aw 

Vj-  sin yjj cosy oos<  ^  -  a)  +  Vj. sin y^cosy sin ( 0  -a) 

V_  cosa» 

+  sin  y  cosy  cos  Sl  sin  ff  —  ct)^{r  +  h)Ji — j-yS  sinSj^  cosyslno 

V-cos<T_  VpSina_  , 

■*■  \  - — Bsinycos^  -.V„[w*  +  .— _i - i — ]• 

r  +  h  ®  (r  +  h)cosSL 

•[cosysin(\&  —  ff)sInSL  ~  sinycosSj^  sin^  ] 

V_  sina„  , 

+  <uJ6j^+«— JL - 2-— ](r  +  h)(cos5i  sSrvSf  cosycosa 

®  *  (r  +  h)cosSL  ''LL 


d  0 


—  cos^Sj^  siny). 


(61) 


and 


“W 


«  ~  V^ycosB  +  Vj 


VLCOS0r_ 

-2 - 2  (  sin  y  sin  a  sinB  +  cosa  cosB) 

r  +  h 


„  *  V„  sino„ 

*- VjffcosysinB  -Vj(6)«+  - s - e — >(sinycosacos&  sinB 

(r  +  h)cosSj^ 

—  sinacos^LCO^B  —  cosysin  Sj^sinB) 


dV, 


-  -j^Vj.sinyn[sinycos{^  -  ct)cosB  -  sin(^  -  <7)sinB] 


+  V, 


^  Vj  sin  yjj  [  sin  y  sin ( ^  -  ff)cosB  +  cos(  -  <7)sinB] 


V  sin(7 

—  i  w®  + - ^  ^ —  1  [  cos  y  sin  cosSlCosB  —  cos(\&  ~  o)sinB 

(r  +  WcosS^ 


+  siny  sln(  tr)  sinS^^cosB]  +  w^Vj-slny  cosSL(slnysinacosB 

V„sinff„ 


V  ef 

—  COS  a  sin  B ) -'  <a^(  r  +  h  )  — a— j-s 


sm^(sinysinacosB  —  cosasinB) 


Vp  sin<7„ 

+  - - - ^)(r  +  h)cos5T  (csyf^osSr  cosB  +  sinosinSr  sinB 

(r  +  WcosS^ 


V_COS<T» 

+  sin  ycosasinSr  cosB)  —  V„ -2 - 2. 

^  (r  +  h) 


(62) 
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By  combining  terms,  sxibstituting  trigonometric  identities, 
and  rearranging,  equation  (50)  is  reduced  to 

-  ■^Sln2y[-^  C0S(V/  -  a)  -  sin(  -  a>  ] 

V  V 

+  —y  g  [slnycosf^ff  ~a  )  -  cos  y  tan  5t  sin  (  0  -<T)slnff_] 
r  +  n  ®  ° 

~  -  Vj.cosyslna)sJnycos5L  +  Vg  sin  5l  cos  y  sin(  <7g  -  (,) 

-  <u^V„[cosysin5Lsin(^&  -  <,)  sin  y cos  5^,  sin  ] 

+  (r  +  h )( cosSj^  sin cos ycos (7  —  cos^ sin y)  (63) 

where  was  replaced  by  y  iii  the  wind  gradient  terms. 


In  most  cases  the  bank  angle  is  kept  small,  and  it 
can  be  assumed  that  sin  y  =  0  and  cos  y  =  1.  With  this 
assumption  equation  (62)  can  be  reduced  by  combining 
terms,  making  substitutions  for  trigonometric  identities, 


and  rearranging  to 


d  V 

“  -Vfcy  -  Vtsin2y[-^cos(./.  -  a)  -  v„:^sin(^  -n)J 


dtp 


V,V-  V  V  , 

+  - - ^  COS(<T  -  ff)  —  ^  c  cosy  cos  (  Ip  ^  \ 

r  +  h  6  -rug/ 

+  sinytanS^sinCi^  -tf)sin0g] 

+  0)^  cos  I  VjSln  c  +  ( Vj  sin^  y  sin  a  +  Vg  cos  >'  sin  Og )  ] 

-  w©V^[cosysin0  cos^L  +  sinysinC^Ji  -a)sin5Ll 

+  Wg,  VgSin  ysinSj^  sin(ag  -  o) 

+  (r  +  hXcosycos^Sj^  +  sin ycosS^  sin  5^0080) 


V..V„ 


(64) 
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where  was  again  replaced  by  y  in  the  wind  gradient 
terms . 

The  terms  containing  the  products  (u^V.  1  /  Iu)-,(V^sino 

w  w  w  9  ^ 

-  V.cosysino)]  and  [w^V  sin{o„  -a)]  in  equations  (20)  and 
u  ^  g  g 

(21)  can  be  neglected  with  little  error.  If  they  are 
dropped  the  equations  reduce  to 

.  V.  f  dVjj,  dtp 

a  ®  U - 5  sln2y[’^  cos( v!r  -  a)  -  V„-~sin(0  -  a)] 

*w  ^2  a'''  dh 

V  V 

+  — S-l sin yco£( 0  —  ff.)  +  cosytanSt  sln(iit  —  a)3inff„] 

r  +  h  ®  u  r  g 

(65) 

+  tu^^(r  ^  h )(cos cosycosa  —  cos^S^siny) 

and 

X  ”  -  Vtsln2y[^cos(0  -  it)  -  v^^sln(V/  -a)] 

Vf  V  V  V 

+  — i^cos(o„  -  it)  -  [cosy cos (i&  —  a_) 

r+h  ®  r+h  ^ 

+  sinytanS^slnC^  -ff)siniTg] 

+  w^cos  5^  [  Vj  sin  IT  +  ( Vj  sln^  y  sin  it  +  Vg  cos  y  sin  iTg )  ] 
f  <u^(r  +  h)(cosycos25j__  +  sin ycosSj^ sin Sl cost*)  (66) 


ENERGY  HEIGHT 

In  the  early  days  of  aviation,  a  climb  amounted  chiefly 
to  increasing  an  airplane’s  potential  energy  with  changes 
in  kinetic  energy  being  quite  small  in  comparison.  As 
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maxlmxuti  speeds  of  airplanes  increased,  particularly  after 
turbojet  engines  came  into  being,  it  became  common  prac¬ 
tice  to  account  for  changes  in  kinetic  energy  when  com¬ 
puting  instantaneous  rates  of  climb  and  determining  "best 
climb"  schedules.  The  concept  of  energy  height  (frequently 
called  specific  energy)  was  introduced  by  German  engineers 
during  World  War  II  and  gained  world-wide  acceptance  in 
dealing  with  the  performance  of  aircraft  powered  by  turbo¬ 
jet  engines. 

Energy  height  is  found  quite  simply  by  considering  the 
total  energy  to  be  the  sum  of  the  potential  and  kinetic 
energies .  That  is 


wVf2  (67) 

E  =  Total  energy  *  WH„  +  - 

®  2gf 

where  total  energy  is  arbitrarily  referenced  to  H  =  0 
and  V.  =  0.  H  is  used  to  compute  potential  energy, 
recalling  that  geopotential  altitude  is  equivalent  to  the 
amount  of  work  done  in  raising  a  unit  mass  from  mean  sea 
level  to  a  geomrtric  altitude  of  h.  Energy  height  is, 
then,  the  total  energy  per  pound  (in  the  English  system) 
of  weight.  Thus 

v,2 

He  “  ~  (68) 
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which  is  independent  of  aircraft  mass.  In  its  classical 
form  energy  height  has  been  defined  as  h  +  V^/2g;  however, 
when  using  a  model  atmosphere  in  which  geopotential  rather 
than  geometric  altitude  is  employed,  II  in  geopotential 
feet  should  be  used  as  a  measure  of  potential  energy. 

Energy  height  is  useful  in  optimizing  climb  performance. 

It  is  energy  height  which  must  be  gained  mostly  rapidly 
(rather  than  altitude)  to  minimize  time  to  climb  since 
potential  and  kinetic  energies  are  readily  interchangeal-ile. 
Also,  energy  height  has  been  used  (reference  the  section 
Standard  Climb  Schedules)  as  an  independent  variable  in 
"optimum"  schedules  during  v/hich  altitude  is  not  monotonic 
increasing.  Differentiating  equation  (68)  we  have 

♦  «  Vf  • 

He  ‘  He  +  (69) 

It  should  be  noted  that  accelerations  arising  from  varia¬ 
tions  in  wind  and  from  the  earth's  rotation  as  seen  in 
equation  (46) ,  for  example,  are  not  accounted  for  in  equa¬ 
tion  (69) . 

As  in  equation  (10)  of  the  section  Atmospheric 
Environment 


gj.dHe  -  gLdh  (70) 
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so  that  equation  (69)  may  be  rewritten  as 


h  +  -iV. 


LOAD  FACTORS 

The  accelerations  from  equations  (65)  and  (66)  are 
substituted  into  equations  (8)  and  (9)  to  obtain  general 
equations  to  describe  load  factors  in  terms  of  geocentric 
parameters.  The  gravitational  and  centrifugal  relief 
terns  are  grouped  together: 


— 1[  gg  -  +  h ) cos2 ]  sin y 

“  +  h)cosSLSingL]cosycosa 

•  V.  ,  dV„  dip 

+  Vj - ^sin2y[-^cos(^  "a)  -  -^sln(ip  -*&)] 

V  V  1 

+  -2L^[sinycos(  -  org)  +  cosytan  sin(  ^  -(r)sini7g]| 


(U^^Cr  +  h)cos2gL3cosy 
+  [  gjj^-  0)^*^  ( r  +  h  >  cos  5^  ^  ^  ^ 

+  Vxy  +  V. sin2y[^^cos(^  —  a)  —  Vy.  -^sin(p  -  a)] 
^  ^  dh  "  db 

y.  V  V  V 

- s_Scos(oc  -  a)  +  [cosy cos -  <7g) 

r+h  ^  r+h  ^ 

+  sinytanSL  sin(\ir  —  <7)sinag]  —  (W^cosSl  [  Vjslno 
+  ( Vj  sln^  y  sin  cr  +  Vg  cosy  sin  Og )  )| 


V~»l 


where  the  geocentric  components  of  the  acceleration 

due  to  gravitational  attraction  (g  and  g_  )  are  computed 

9  9 

from  equations  (1)  and  (2) . 


With  zero  wind  equations  (72)  and  (73)  reduce  to 

Ujj  -  — “  <u^^(r  +  h)co825L]8iny 

w  V  H 

"  [g*  “  6)J^(r +  h)co8SLsln5Llco8yco8a  + Vj| 


— |[gz  -  ^.ii)co82SLlco8y 

gf  I  g 

+  [g^  -  <a^^(r +  h)co8SL8i“5L^®^“^®°®® 

/  V^2  I 

+  V^y  - cosy  —  2a»^Vj.co8  5jj sin (t| 


Since  the  local  centrifugally  relieved  acceleration 
due  to  gravity  is  approximately  normal  to  the  geoid, 
equations  (72)  through  (75)  can  be  rewritten  including 
the  total  resultant  gravity  vector,  gj^; 


f  .  V,  dV* 

J-jgLsinyjj  +  Vt  -  -^8in2y[-^c08(f  -  a) 
gf  I 

-  V^~sln(^  ~a)]  +  [slnycos(\ir  -  a  ) 

dh  rvh  ^ 

f  cosytan 5l 0  — a)8lnag]| 


and 


if  *  2  dV 

%  -  +  Vfcy  +  VjSin  Yl  ~-50s(^  -  <r) 

dd,  V^V 

■"  -<^)] - —  C0S(<jr_  -  or) 

dil  r  +  h  * 

V  .V 

+  [cosycos(^  +  sinytanSj^sinC^  -a)sinffg] 

-  oj^cos  S|^  t  Vj.  sin  CT  +  (Vjsin^yslnff  +  Vg  cosy  sin  ffg)]|  (77) 

For  zero  wind 


•  Vi.'' 

l.(gLCosyjj  +  V^y  -  — s~cosy  -  2(i>g>V^cosfLSino)  (79) 


The  load  factors  of  equations  (72)  through  (75)  are  not 
exactly  equal  to  those  computed  with  equations  (76) 
through  (79)  because  of  the  approximations  in  the 
derivations  including  neglect  of  the  differences  in  the 
roll  and  heading  angles  of  an  aircraft  when  referenced 
to  the  geocentric  instead  of  the  geodetic  horizontal 
plane. 


It  is  frequently  convenient  to  relate  the  longitu~ 
dinal  load  factor,  n^^  ,  to  energy  height.  This  may  be 
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f 

V 


done  as  follows:  by  rearranging  equation  (69}  and  sub** 
stituting  equation  (33)  we  have 


(80) 


Since  the  right  hand  side  of  this  equation  appears  in 
equations  (76)  amd  (78) ,  the  left  hand  side  can  be  intro¬ 
duced  into  these  equations  to  produce 


ACCELERATION  FACTORS 

Equation  (71)  can  be  rearranged  to  produce  equations 
for  geometric  rate  of  climb,  h,  in  terms  of  either  alti¬ 
tude  or  energy  height.  Two  different  acceleration 
factors  result.  First,  expanding  by  the  chain  rule, 
assuming  h  0,  equation  (71)  becomes 


V-44 


(83) 


% 


gr  •  V*  dVj. 

-ih  +  -Ji _ 5h 

dh 


Solving  for  h 


% 


V*  dV, 


«L  dh 


(84) 


The  denominator  is  defined  as  the  acceleration  factor, 

Aj,  and  is  used  with  any  of  the  continuous  climbs  described 
in  the  section,  Standard  Climb  Schedules. 


«T 


V,  dV, 

+  -i — i) 
Kl  dh' 


(85) 


With  no  wind,  this  provides  the  following  equation  for 
rate  of  climb  using  equation  (82) 


h 


(86) 


The  second  acceleration  factor  is  useful  when  climb 
data  are  standardized  at  constant  energy  height  (reference 
option  4,  page  9,  Appendix  I),  Again  expanding  in 
equation  (71)  but  assuming  H^,  0 
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I 


«r  Sr  dHfe  ^ 


Solving  for  h 


8l  d% 


The  acceleration  factor  is,  then 


iLn  -  Zl£XL) 

8L  '8r 


With  no  wind  we  have  from  equation  (82) 
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ACCELEROMETER  METHODS 


In  the  preceding  paragraphs  the  general  equations  for 
computing  inertie-.l  accelerations  and  excess  thrust  of  an 
aircraft  were  developed  for  the  airspeed-altitude  method. 
Using  this  method  a  number  of  individual  terms  must  be 
evaluated  in  order  to  obtain  the  total  accelerations  and 
resulting  forces  (e.g.#  airspeed  derivative,  rate  of  climb, 
wind  gradients,  Coriolis  accelerations,  centrifugal  relief). 
Evaluation  or  some  of  the  terms  involves  numerical  differ¬ 
entiations,  and  the  results  are  dependent  on  the  numerical 
data  editing  and  differentiation  methods.  In  addition  the 
basic  parameters  airspeed  and  altitude,  obtained  with  con¬ 
ventional  instruments,  have  limited  accuracies  which  are 
degraded  due  to  complications  of  position  error  and  pressure 
lag.  Also,  it  is  difficult  to  evaluate  wind  gradients  pre¬ 
cisely  because  of  uncertainties  about  fluctuations  of  atmo¬ 
spheric  conditions  with  time  of  day  and  with  distance  from 
the  aircraft  flight  corridor.  The  same  numerical  problems 
of  data  editing  and  differentiation  are  also  encountered 
when  weather  balloon  position  coordinates  are  successively 
transformed  to  wind  velocity  and  gradients. 

In  light  of  these  uncertainties,  methods  have  been 
develOt>ed  to  measure  an  aircraft's  accelerations  more 
directly  using  accelerometers.  These  methods  are  not  with¬ 
out  their  own  problems,  but  they  do,  in  general,  provide 
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much  superior  values  of  excess  thrust.  Accelerometers  sense 
the  inertial  or  total  accelerations  acting  on  an  aircraft 
and  their  readings  can  be  converted  directly  into  forces 
by  multiplying  by  the  aircraft  weight.  Consequently ,  two 
of  the  major  problems  of  the  airspeed-altitude  method  can 
be  eliminated;  data  editing  and  differentiation,  as  well  as 
measurement  of  atmospheric  winds.,  are  not  required. 

ACCELEROMETER  IJISTALLATIOD^S 

Sensitive  accelerometers  have  been  installed  on  test 
aircraft  and  have  produced  excess  thrust  data  of  signifi¬ 
cantly  better  precision  than  similar  data  from  measurements 
of  airspeed  cind  altitude.  (See  reference  5,  for  example.) 


Vane  Mounted 

With  an  accelerometer  mounted  on  a  vane  (similar  to  an 
angle  of  attack  vane  to  keep  it  alined  with  tiie  local  flow)  , 
excess  thrust  could  be  found  immediately,  knowing  aircraft 
weight ,  from 

F  =  n  W 

*w  (93.) 

The  local  flow,  however,  is  not  coincident  with  the  airplane's 
velocity  vector  because  of  upwash,  pitch  rate,  etc.  (Correc¬ 
tions  are  derived  in  subsequent  paragraphs.)  Further, 
mechanical  misalinement  of  the  sensitive  axis  of  the  acceler¬ 
ometer  relative  to  the  vane  may  be  expected.  Hence,  the 
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sensitive  axis  of  the  accelerometer  is  displaced  from  the 
local  flow  at  the  vane  (and,  therefore ,  the  longitudinal 
axis  of  the  vane)  by  the  misalinement  angle,  Since 

accelerations  in  the  wind-axes  must  be  found  in  order  to 
compute  excess  thrust,  transformations  through  the  angle 
Ema  other  correction  angles  are  required.  These 
transformations  are  carried  out  in  a  later  paragraph  under 
Errors  in  Measured  Accelerations. 


Fixed  eg  Mounted 

Rather  than  installing  a  two-axis  accelerometer  system 
on  a  vane,  it  may  be  hard  mounted  near  the  aircraft's  cen¬ 
ter  of  gravity.  This  has  the  advantage  that  errors  in 
measured  accelerations  caused  by  changes  in  attitude  arc 
made  negligible;  however,  since  those  corrections  are  gen¬ 
erally  small  and  can  be  made  easily,  thi^  advantage  is  a 
slight  one.  If  load  factors’  are  measurea  wildi  acceler¬ 
ometers  located  near  the  center  of  gravity,  they  must  be 
transformed  through  any  misalinement  angle  representing 
the  displacement  of  the  sensitive  axis  of  the  longitudinal 
accelerometer  from  the  airplane  body  axis  ana  through  the 
angle  of  attack.  This  permits  loads  in  the  direction  of 

and  normal  to  the  velocity^ vector  to  be  found. 

r  ■ 

I  ; 

The  sources  of  error  in  the  position  of  a  vane  to  which 
an  accelerometer  system  is  attached  also  exist  in  measured 
angle  of  attack  when  sensed  by  a  vane.  They  appear  also 
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in  values  obtained  from  a  pressure  sensing  device  with  the 
exception  of  dynamic  response.  To  avoid  degrading  the 
accuracy  inherent  in  the  accelerometer  system/  correc¬ 
tions  for  the  errors  described  in  the  following  para¬ 
graphs  should  be  made  and  angle  of  attack  should  be 
accurately  known  (to  within,  say,  0.1  degrees)  when 
using  an  accelerometer  system  mounted  near  the  airplane's 
center  of  gravity. 

ERRORS  IN  MEASURED  ANGLE  OF  ATTACK 

Angle  of  attack  is  generally  sensed  with  a  vane 
mounted  on  a  nO;3e  boom  well  ahead  of  the  aircraft.  The 
accuracy  of  the  sensed  value  (angular  displacement  of  the 
relative  wind  from  the  airplane  body  axis)  is  adequate. 
Substantial  corrections,  however,  must  be  applied  in 
order  to  find  true  angle  of  attack.  These  corrections 
arise  xrom:  bending  of  the  boom,  pitching  velocity  (which 
adds  a  component  of  velocity  to  the  vane  not  experienced 
by  the  airplane's  center  of  gravity),  upwash  (created  by 
the  presence  of  the  boom,  nose  of  the  airplane,  and  its 
wing),  and  lag  in  the  vane  position  caused  by  rapid  motions 
abo It  the  pitch  axis. 

Angle  of  attack  has  been  determined  less  frequently 
through  r  differential  pressure  sensing  device  attached  to 
the  boom  to  obtain  AP/q^,  which  can  be  related  to  indicated 
angle  of  attack.  Corrections  similar  to  those  described 
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for  the  angle  of  attack  vane  need  to  bis  made  for  this  sort 
of  installation  with  the  excepcion  of  dynamics  response. 

Boom  Bending 

Loads  from  aerodynamic  forces  and  from  inertia  cause 
bending  of  the  boom  which  results  in  errors  in  the 
measured  angle  of  attack  since  the  positron  of  the  vane 
is  referenced  to  the  axis  of  the  boom.  Bending  from  aero¬ 
dynamic  loads  is  usually  negligible  although  it  may  be 
estimated  from  data  contained  in  reference  6.  Adjustments 
to  boom  bending  for  loads  due  to  inertia  may  be  made  from 
a  calibration  of  static  deflections  of  the  boom  when  loaded 
with  weights  to  represent  inertial  forces  experienced  in 
flight. 

Upwash 

The  largest  correction  to  be  made  in  finding  true  angle 
of  attack  generally  stems  from  upwash.  As  previously 
pointed  out,  upwash  is  generated  by  the  boom,  the  nose  of 
the  airplane,  and  the  wing.  The  upwash  is  most  pronounced 
at  high  angles  of  attack  and  decreases  with  lift  coefficient. 
At  supersonic  speeds  effects  of  fuselage  ana  wing  disappear, 
of  course.  Upwash  generated  by  the  boom  may  be  measured  in 
a  wind  tunnel. 

Two-dimensional  incompressible  flow  theory  has  been  used 
to  compute  the  effects  of  upwash  at  low  speeds.  Figure  3, 
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Thitrttical  Efffef  ts  on  Angle  of  Attack  Measurements  of  Upwash  frem  the  Nose 
Boom  and  Fuselage  at  Low  Speeds 


eiWKiwn 


taken  from  reference  7  page  16/  shows  the  results  of  such 
computations.  Indicated  in  the  figure  are  variations  in 
angle  of  attack  from  the  boom  itself  and  from  the  fuselage, 
assumed  to  be  a  blunt  circular  cylinder. 

In  reference  8  equations  are  presented  which  may  be 
used  to  estimate  the  effects  on  angle  of  attack  of  the 
fuselage  and  wing.  Effects  of  the  fuselage  are  found  using 
potential  theory,  considering  the  fuselage  to  be  a  half 
rotational  body;  influence  of  the  wing  is  computed  assuming 
a  bound  vortex  in  the  quarter-chord  line  of  the  wing  sec¬ 
tion  and  solving  for  the  induced  vertical  velocity  with  a 
Biot-Savart  equation.  Estimated  upwash  angles  have  been 
computed  using  references  7  and  8  and  compared  to  in-flight 
calibrations  obtained  with  an  A-37B  (reference  9).  Corre¬ 
lation  of  the  estimated  angles  with  those  from  the  calibra¬ 
tions  was  reasonably  good.  In  flight  test  applications, 
however,  it  is  desirable  to  measure  upwash  angles  from 
flight  data  and  construct  a  calibration  of  upwash  angle, 

as  a  function  of  lift  coefficient. 

Dynamic  Response 

An  angle  of  attack  vane  system*  constitutes  a  torsional 
spring-mass-damper  mechanical  system  having  an  undamped  natural 

*The  vane  system  includes  any  internal  mass  balancing,  transducer 
elements,  accelerometer  package,  etc.,  in  addition  to  the  aero¬ 
dynamic  lifting  surface. 
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frequency  and  damping  rauio  that  describe  its  dynamic  response 
characteristics.  The  output  of  such  a  system  will  have  errors 
that  are  due  to  signal  amplitude  and  phase- lag  characteristics 
which  are  functions  of  the  frequency  relationship  of  the  exciting 
signal  and  the  natural  frequency  of  the  vane  system  as  well  as 
its  damping  ratio. 

Equations  are  derived  below  which  may  be  used  to  correct 
for  errors  in  sensed  angle  of  attack.  The  same  equations  may 
also  be  used  to  find  errors  in  the  position  of  vane-mounted 
accelerometers . 


Response  to  Sinusoidal  Inputs 

An  angle  of  attack  vane  system  is  shown  schematically 


About  the  vane  pivot,  taking  cos  A a  =  1 


u 


Lyi  -  I. 


•y“ic 


(92) 


which  may  be  rewritten  as 


■  ^®lc 

a  V. 


(93) 
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Rearranging 


Equation  (94)  is  a  second  order  differential  equation 
describing  a  system  whose  dynamic  characteristics  are 
given  by 


■y 


2  46)n 


and 


It  should  be  noted  that  the  natural  frequency  and 
damping  of  the  sytem  are  not  unique  properties  but  are 
dependent  on  Mach  number  and  altitude.  Thus  these  prop¬ 
erties  must  be  computed  for  each  flight  condition  at 
which  data  are  to  be  collected. 

Since  the  angle  of  attack  vane  system  is  approximately 
represented  by  a  second-order  dynamic  system,  the  first 
method  for  correcting  angle  of  attack  for  dynamic  lag 
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will  assume  a  sinusoidally  varying  angle  of  attack.  This 
is  a  fairly  realistic  assumption  for  an  aircraft  engaged 
in  roller  coaster  maneuvers.  Under  these  circumstances , 
then,  in  Laplace  notation  considering  only  the  error  due 
to  lag 


.2 

s 


+  24^<t>QS 


+ 


(98) 


where  K  is  a  fixed  gain,  usually  1.0. 

Rearranging  this  equation  and  introducing  ju  for  the  Lap¬ 
lace  operator,  s,  yields 


(99) 


This  relationship  defines  the  amplitude  and  phase  relation¬ 
ship  between  the  two  sinusoidal  oscillations  in  a  and 
Solving  equation  (99)  for  a  we  have 


{[1  -  +  {2^^)2F 

I  ‘“n  “n  )  j  tan 

«ic  -  e 


2C 


0) 


-1  <“] 


(100) 


Introducing  for  the  expression 


«ic  “  «ic„sinwt 


(101) 


where  a.  is  the  maximum  amplitude  of  a.  ,  it  is  possible 

ICq  1C 

to  rearrange  the  above  equations  so  that  a  solution  is 
given  for 


-  at-  rsin(<ut  +  tj) 
a  iCq 


(102) 


where 


[1-(£1_)2  + 


|[1  -  (^)2]^+  (2Cf-)2]p 
I  <“11  _ * 


-  2  cos  [tan' 


-1  ^  <^n  A 

1  _  (SL.)^  / 


(103) 


2^£i- 
1  -(^)2 


cos  [  tan" 


1  - 

‘"n 


(104) 


There  exists  another  method  for  finding  Aa„  for  the 

V 

case  where  sinusoidal  inputs  may  be  assumed  and  which 
involves  fewer  calculations.  Given ^  as  before,  that  the 
angle  of  attack  vane  is  a  second  order  system  with  an 
undamped  natural  frequency  of  oscillation,  and  a 
damping  ratio,  the  indicated  angle  of  attack  is  given 


aic(t) 


KaQSin(<at-(^) 


([1  -  (^] 


)2]2  +  (2C^) 


(105) 


where  K  is>  as  before,  a  constant  gain  (usually  1.0)  and 
is  the  maximum  amplitude  of  a. 

Making  use  of  trigonometric  identities 


aic(t) 


K 


«o 


j[l  -(^)2]2  +(2^^)2f 


.(sineutcos^  ->  coscutsin^) 


(106) 


where 


4> 


tan 


-1 


2C^ 


I .  <1^,^ 


(107) 


By  rearranging  the  terms  in  this  equation  and  making  small 
angle  approximations  (cos  (|)  =  ls  in  (!>  =  (})) 


a(t) 


{[l-(SL)2f  +  (2C^y^f 

'•  ‘“n  ^  r 

— — - — -  la 


K 


ic(t)  + 


ajp  ^coswt] 


(108) 


In  the  above  equation  wt  describes  the  sinusoidally 
varying  a.  This  is  generally  not  known  since  a  is  not 
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known;  however,  it  can  be  inferred  through  knowledge  of 
aic(t)  and  a  small  correction  to  t.  When  =  0,  it 

can  be  shown  that  tan  wt  =  <|)/  where  wt  describes  a(t) , 
rather  than 

Since  t  =  (1/w)  tan“^<J)  when  oii_(t)  =  0,  it  can  be 

i  C 

used  to  correct  the  time  associated  with  the  observed 
angle  of  attack,  aj_(t)  so  that  the  expression  for  a(t) 
becomes 


/[I  +  (^)2]2  +  . 


(109) 


where 

a.  (t  )  =  indicated  angle  of  attack  at  any  time,  t 

“ic  “ic 

®lc  *  maximum  amplitude  of  observed  angle  of  attack 


t  =  time  from  a  reference  value  of  the  indicated 

®ic  angle  of  attack 

At  =  incremental  time  correction  to  t 

(l/w)tan“‘(j)  ic 


=  phase  lag  angle  between  a(t  ) 


and 


Response  to  Reindom  Inputs 

The  previous  discussion  was  concerned  with  correcting 
indicated  angle  of  attack  when  a  sinusoidal  variation  of 
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the  input  had  been  assumed.  In  the  more  general  case,  as 
in  a  climb,  the  angle  of  attack  is  changed  in  a  series  of 
quasi-step  inputs  which  are  aimed  at  keeping  the  aircraft 
on  a  desired  climb  schedule.  In  correcting  the  indicated 
angle  of  attack  in  this  case,  a  somewhat  more  complicated 
process  is  used.  Since  the  angle  of  attack  vane  system 
is  a  second  order  system,  its  input-output  relationship 
is  given  by 


2  2 
s  +  2  4^6)S  + 


(110) 


If  the  operations  indicated  by  the  above  Laplace- 
transform  expression  are  performed,  the  following  expres¬ 
sion  in  the  time  domain  results: 

*•  *2 

a(t)  »  — i 

In  general  it  will  be  necessary  to  operate  on  the 
indicated  angle  of  attack  time  history,  a.  (t) ,  with  a 
computer  program  designed  to  give  first  and  second 
derivative  information  about  a  variable.  One  such  program 
is  called  DIRSIT  and  is  described  in  reference  12.  In 
general,  both  first  and  second  derivative  corrections  to 
indicated  angle  of  attack  time  histories  should  be  mads 
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since  omitting  the  correction  for  the  second  derivative 
can,  in  some  cases,  cause  the  end  result  to  be  as  much  in 
error,  but  in  a  different  sense,  as  having  made  no  connec¬ 
tion  at  all. 

Pitch  Rate 

The  aircraft  airspeed  vector  expressed  in  the  wind- 
axes  is,  referring  to  figure  5 

Vj  -  vj  (112) 


Fjguri  S  OHentatian  of  Vano 

To  compute  the  vane  airspeed  vector,  V^,  it  is  assumed 
that  the  Xj^  axis  passes  through  the  vane.  The  following 
equation  may  then  be  written  (reference  10,  page  440) ; 
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Vj  +  Wfa  X  iy 


(113) 


V 


When  an  aircraft  is  subjected  to  a  pitch  rate  there 
will  be  a  difference  between  the  true  airspeed  of  the 
aircraft's  center  of  gravity  and  the  local  airspeed  of 
a  vane  on  a  nose  boom.  Since  the  vane  alines  itself  with 
the  local  flow^  the  angular  alinement  error  is  the  angle 
between  the  aircraft  airspeed  and  the  vane  airspeed 
vectors . 

To  impute  the  vane  airspeed  vector,  the  follow¬ 
ing  equation  is  used  (reference  10,  page  440) 
where  the  body  angular  velocity  is 

6)jj  ■  pijj  +  qjf,  +  rkjj  (114) 

Recognising  that 


we  may  write 


V 


V 


•lb  I'b 

q  r 

0  0 


(115) 


(116) 
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Expanding  the  determinant 


Vj.  +  [^yTjjj  - 


(117) 


To  express  the  velocity  in  the  wind-axes  system  a  trans¬ 
formation  is  made  through  a  in  the  negative  direction. 
Making  use  of  the  matrix 


cosa  0  sin  a 

0  1  0 

-sina  0  cosa 


taken  from  the  section.  Coordinate  Systems  and  Transfor¬ 
mations  ,  tt.e  vector  in  brackets  in  equation  (117)  becomes 


cosa  0  sina 

0 

-  ^yqsina 

0  1  0 

Ct 

y.yr 

-sina  0  cosa_ 

’v'J  - 

-  ’yqcosa 

(118) 


Tn  t^^  /ind-axes,  equation  (117)  becomes 

Vy  ”  VjT  +  [  -  ^yqsinal  +  c^yrj  -  p.yqcosak]  (119) 
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Adding  coinponent<i  in  equation  (119) 


Vy  -  (Vj  -  J-yasina)!  +  f.yrj  -  iiyQcosak  '  (120) 

Remembering  that  the  unit  vector  Ic  is  positive  downward 
and  the  angle ,  is  positive  upward 


tan  Aa 

q 


-(-  iiyqcosa) 
Vj  -  iyqsina 


(121) 


or 


Aa„ 


1  P.yQcosa 

tan”^ — I - ; — 

Vj-  O.yqsino 


(122) 


Since  a  is  not  known ^  an  iteration  procedure  has  to  be 
used  to  compute  Aa^.  In  the  first  iteration  the  correc¬ 
tions  described  above  are  added  to  the  indicated  angle  of 
attack  so  that  the  angle 

“ic  ^®boom  bending 

is  substituted  in  place  of  a  in  equation  (122) . 

True  Angle  of  Attack 

True  angle  of  attack  is  determined  by  adding  corrections 
for  the  errors  described  above  (boom  bendihg,  upwash,  dynamic 


response,  and  pitch  rate) .  I£  a  pressure  sensing  device  is 
used,  dynamic  corrections  may  be  ignored  since  there  are  no 
moving  parts.  Upwash  corrections  are  generally  the  largest, 
peucticularly  at  low  speeds.  Corrections  for  dynamic  response 
and  pitch  rate  are  frequently  not  required  but  may  be  sizable 
during  roller  coaster  maneuvers,  and  during  climbs, 
especially  for  high  performance  aircraft.  If  values  of  pitch 
rate  and  pitch  acceleration  are  needed,  they  are  best  obtained 
by  direct  measiurement  with  on-board  instrumentation.  In  the 
event  that  the  instrumentation  is  not  available,  the  required 
data  may  be  confuted  using  equations  derived  under  the 
heading,  Attitude  Rates,  in  the  paragraph.  Errors  in  Measured 
Accelerations . 

ERRORS  IN  MEASURED  ACCELERATIONS 

Two  types  of  accelerometer  installations  have  been  con¬ 
sidered:  one  in  which  the  accelerometers  are  attached  to  a 
vane  located  on  a  nose  boom  and  the  other  with  the  acceler¬ 
ometers  hard  mounted  in  the  aircraft  near  its  center  of 
gravity.  In  both  cases  a  mechanical  misalinement  should  be 
considered.  For  accelerometers  on  a  vane  the  sensitive  axis 
of  the  longitudinal  accelerometer  may  not  coincide  with  the 
axis  of  the  vane.  Similarly,  for  eg  mounted  accelerometers, 
the  sensitive  axis  of  the  accelerometer  may  not  coincide  with 
the  airplane's  body  axis,  Xj^. 
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In  addition  to  accounting  for  mechanical  misalinementr 
axis  transformations  must  be  made  to  find  load  factors  in 
the  wind-axes  system.  For  vane-mounted  accelerometers  the 
rotation  is  made  through  the  net  angle  which  results  from 
combining  errors  due  to  boom  bending,  upwash,  dynamic  lag, 
and  pitch  rate.  In  the  case  of  fixed  eg  mounted  acceler¬ 
ometers  the  rotation  is  made  through  the  angle  of  attack, 
but  the  same  errors  must  be  evaluated  since  they  are  used  to 
compute  angle  of  attack  (excepting  dynamic  lag  when  angle  of 
attack  is  determined  from  a  differential  pressure  instrument) . 

With  accelerometers  on  a  veuie,  accelerations  are  induced 
at  the  vane  by  emgulau:  attitude  rates  and  accelerations  which 
are  not  experienced  by  the  accelerometers  when  located  at 
the  center  of  gravity.  If  the  accelerometers  are  on  board 
the  aircraft  near  the  center  of  gravity  similar  accelera¬ 
tions  are  generated,  but  the  distemce  from  the  center  of 
gravity  to  the  accelerometer  should  be  small  enough  to 
make  the'^induced  accelerations  negligible.  However,  they 
should  not  be  eliminated  arbitrarily  but  only  after  inves¬ 
tigation  of  the  particular  situation  shows  these  induced 
accelerations  to  be  trivial. 

An  additional  source  of  error  is  incurred  with  acceler¬ 
ometers  on  a  vane.  Dynamic  lag  resulting  from  rapid  con¬ 
trol  inputs  affects  the  alinement  of  the  vane  (as  well  as 
an  angle  of  attack  vane)  relative  to  the  local  flow. 
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Mechanical  Misalinement 


It  is  inevitable  that  there  will  be  some  angular  mis- 
alinement  between  the  axes  of  a  vane  and  the  sensitive  axes 
of  accelerometers  (both  normal  and  longitudinal)  attached 
to  the  vane.  Considering  the  construction  of  the  acceler~ 
ometer  system/  it  is  generally  secured  as  .  ^wo-acceler- 
ometer  package  with  their  orthogonality  well  enough  con¬ 
trolled  so  there  is  little  error.  Also/  since  acceleration 
in  the  y-direction  is  not  measured/  misalinement  from  rotar 
tion  about  the  z-axis  is  not  taken  into  account.  The  prin¬ 
cipal  source  of  misalinement  is,  then,  an  angular  rotation 
of  the  accelerometer  case  about  the  vane  axis  as  illustrated 
in  figure  6. 


accelerometer 


Figure  6  Vane  Misalinement  Angle 
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Components  of  aeceleratlon  alon^^ahd :  nomal  to  the 
vane  are  found  by  transforming  through  the  angle  to 

obtain 

\  -  (123) 

and 

"zy  "  (124) 

where  the  subscript  m  refers  to  measured  values.  In  order 

t(.  retain  the  accuracy  inherent  in  sensitive  accelerometers, 

it  is  necessary  to  know  the  magnitude  of  quite  well.  To 

illustrate,  suppose  the  n  =  Og,  n_  =  Ig,  and  e  ^  ~  1 

m  m 

degree,  n  ,  arising  from  the  component  of  n  is  0.01745g. 
^v 

Since  this  is  several  times  tiie  accuracy  which  may  be 

achieved  in  the  measured  quantities,  the  misalinement  angle 

should  be  known  to  within  about  0.1  degrees. 

Use  of  a  sensitive  accelerometer  to  measure  normal 

load  factor  is  desirable  in  computing  n  ,  (and  hence 

•'v  ^  particularly  for  climbs  or  other  maneuvers  during 
*w 

which  there  are  significant  changes  in  normal  load  factor. 

A  similar  misalinement  angle  should, be  expected  when 
the  accelerometer  system  is  mounted  near  the  center  of 
gravity  of  the  airplane.  In  this  case  the  misalinement 
angle  is  the  angular  displacement  of  the  sensitive  axis 
of  the  longitudinal  accelerometer  from  the  airplane's  body 
axis,  Xj^.  An  axis  transformation  similar  to  that  indicated 


V-6B 


by  equations  (123)  and  (124)  is  made  to  obtain  load  factors 
in  the  body  axes. 


(125) 


and 


(126) 


Observed  Load  Factors  in  the  Wind-Axes  System 

It  is  necessary  to  mcdce  axes  transformations  through 
the  angle  between  the  sensitive  axes  of  the  accelerometers 
and  the  wind  axes  for  both  types  of  accelerometer  systems 
considered.  This  might  be  done  by  making  the  transforma¬ 
tion  through  a  single  angle  in  both  cases.  The  resultant 
load  factors  have  been  found,  however,  by  first  transform¬ 
ing  through  the  misalinement,  as  indicated  in  equations 

(123)  through  (126) .  In  the  case  of  vane-mounted  acceler¬ 
ometers,  then,  an  additional  transformation  through  the 
angle 


net 


^*boom  bending  ^ 


Aa  +  Aa  .  +  Aa 


<t> 


(127) 


is  required,  and  observed  load  factors  in  wind  axes  become 


(128) 


and 


«z,  °o®^“net  * 


(129) 


The  load  factors  in  equations  (128)  and  (129)  are  in 
the  wind-axes  system  but  have  not  been  corrected  for 
induced  accelerations  derived  in  the  following  paragraph. 

For  accelerometer  systems  hard  mounted  near  the  eg., 
the  desired  load  factors  may  be  found  immediately  by  making 
a  transformation  through  the  angle  of  attack ,  which  is 
computed  from 

a  ■  CL,  +  Aa.  u  jj  (130) 

1C  boom  bending  “u  <a 

We  have ,  then 


n_  «  Hy  cosa  —  n_  slna 
*w  *b  ^b 


and 


n 


Zw 


n,  slna  +  cosa 
*b 


(131) 


(132) 


Accelerations  Induced  from  Attitude  Rates 

The  acceleration  at  the  vane  referenced  to  that  at  the 
eg  is,  from  reference  10,  page  443 


^cg 


+  ^ 

V  0  *^0  r.\ .  V 

1 _ 

dt^ 

dt 

b  L  J 

A  Xf  y  • 

h 

dt  J 

+  SijXC^b*  *’'v)  (133) 


where  the  body -axes  system  is  moving  and  the  terms  in  the 
brackets  represent  derivatives  within  the  moving  system. 
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Since  the  axis  of  the  vane  is  assiuned  to  lie  on  the 
axisr  the  first  and  second  derivatives  of  are  zero. 


[€  ■  [Si 


(134) 


Evaluating  the  remaining  terms 


dt 


(135) 


&1-.  ■ 
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p  a 
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•  . 


r 
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(0)lb  +  (  Vr)jb  -  (^q)kb  (136) 


and 


‘b  Jb  kb 

•  ♦  • 

P  q  r 

0  -Hyq 


-.f.y(q2  +r2)Tb  +  ^^qpjb  +  JlyP^^b  (I37) 


Substituting  equations  (135)  and  (137)  in  equation  (133) 


ay  *  ^Cg  +  (P**  ^^^^b  +  (P'^“^^kbl  (138) 


Transforming  the  components  to  the  wind- axes  system,  as 
in  equation  (118) 
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cos  a  0 


—  sina  0 


Sinai  +  r2) 


0  yap  4  f) 

cos  a  *>y(pr  ~  a) 


+r^)oosa  +  <pr— Q)sina] 
£v<«P  +  f> 

JlvUa^ +r^)8i“«  +  (pr-Q';cosa] 


(139) 


The  final  equation  for  acceleration  in  the  wind-axes  is 


Sy  =  a^jg  +  y|[ -<q2 +r^)cosa  +  (pr-q)sina]i  +  (qp+f)T 

+  [(q^  +  r2)sina  +  (pr— q)cosa]k| 


(140) 


The  induced  acceleration  is  made  up  of  the  terms  in 
equation  (114)  representing  the  difference  between  the 
acceleration  at  the  vane  and  that  at  the  center  of  gravity r 
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Corrected  Load  Factors  in  the  Wind-Axes  System 

When  induced  accelerations  as  defined  in  equation  (140) 
become  large  enough  to  have  significant  effects  on  load 
factors,  it  becomes  necessary  to  calculate  load  factors 
at  the  eg  in  terms  of  the  load  factors  measured  by  the  accel¬ 
erometers  located  on  a  vane. 

The  correction  to  load  factor  is 


•  ( ftgg  ) 


(141) 
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so  that,  from  equations  (128) ,  (129)  and  (140)  the  load 
factors  along  and  normal  to  the  airplane's  velocity  vector 
are 

%  “  tij^cosAonet  -  n^^sInAonet  +  ^[(q2  +  r2)cosa  -  (pr-q)sina]  ^43) 


and 


X, 

Uj^cosAonet  +  “xySiaAonet  +  ^[{q^  +  ;^>sina  +  (pr-a)cosa]  (143) 


Since  these  corrections  are  not  usually  needed  for 
accelerometers  located  near  the  eg,  load  factors  for  this 
sort  of  installation  may  be  found  from  equations  (131)  and 
(132) . 

If  the  corrections  to  load  factors  described  above  (and 
to  angle  of  attack  for  pitch  rate)  are  made,  the  best  way 
to  obtain  the  attitude  rates  and  attitude  accelerations  is 
to  measure  them  directly  with  rate  gyros  and  angular  accel¬ 
erometers.  Any  other  method  involves  differentiation  which 
will  magnify  any  errors  arising  from  inadequate  instrumen¬ 
tation,  calibration,  data  reduction  technique,  etc.  Despite 
these  limitations,  corrections  can  be  made  by  use  of  equa¬ 
tions  presented  in  the  following  paragraph  although  it  is 
iterated  that  the  results  can  be  expected  to  be  inferior  to 
direct  measurements. 
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Computed  Attitude  Rates 
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Equating  the  vector  components  of  equations  (114)  and  (146) 


p  =  (B  -  CTSin>')cosa  -  (-ysinE  a cO!sycosB)sina  (147) 


q  =  a  +  ycosB  +  acosysinB 


(148) 


r  =  (B  —  asiny)sina  +  ( —  ysinB  +  (7COsycosB)cosa  (149) 


The  parameters  in  these  equations  can  all  be  computed 
from  airspeed,  altitude,  and  angle  of  attack  except  for  B, 
B,  and  o.  To  evaluate  these  unknowns  B  is  assumed  to  be 
zero;  then  the  other  two  parameters  can  be  obtained  by 
reference  to  the  equations  of  motion  along  the  wind  y-  and 
z-axes.  The  accelerations  along  these  axes  are,  neglect¬ 
ing  smaller  order  terms,  from  equations  (47)  and  (48) 

"  Vj(ffcosycosB  —  ysinB) 

and 

a„  =  —  V,(ycosB  +  ffCosysinB) 

Zw  '' 

The  components  of  the  local  gravity  can  be  obtained  by 
transforming  the  component 

(152) 


(150) 


(151) 
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in  the  local  geocentric-axes  system  into  the  wind-axes 
system.  Applying  the  transformation  matrix  from  the  section. 
Coordinate  Systems  and  Transformations ,  results  in 

®L  “  gL(-sin}'r  +  cosysinBj  +  oosycosBk)  (153) 

Recalling  that  the  normal  load  factor  is  positive 
upwards  opposite  to  the  positive  z-axis,  the  total  load 
factor  at  the  eg  is 

"eg  =  nJ  +  nyl-n^lT  (154) 

The  vector  equation  for  load  factor  in  terms  of  air- 
craf':  acceleration  is 

rigg  =  -^(^cg  ”■  (155) 

Substituting  the  components  from  equations  (150) ,  (151) , 

(153)  and  (154)  in  equation  (155)  and  picking  out  the  y 
and  z  components 

n  =  V,( ‘jcc.sycosB  —  ysi^B)  —  g.  copysInB]  (156) 

gj.  ^ 

and 

n  =  JL[V.(ycosB  +  ffCosysinB)  +  g.  cosycosB]  (157) 

gj. 

With  the  assumption  of  zero  sideslip  the  side  load 
factor  must  be  zero.  Therefore,  from  equation  (156) 
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Yj(ac's)'cosB  —  /sinB) 


cos  y  sir;  3 


(158) 


Equ«ic;.on  (158)  can  be  solved  for  tan  B  to  give 


tan  3 


Vj,acosy 
gj^cosy  +  v^y 


(159) 


Rearranging  equr-.  .wi.  (157) 

g,n_  — '  gr  cosycos  B  •  '■/-'''osB  =  ( V%acosy)sinB  (160) 

Solving  equation  vI59)  for  (V^ocosy) ,  substituting  into 
equation  (160) ,  and  multiplying  by  cos  B 

gj^Hg  cosB  ~  (gjjcosy  +  Vjy)cos^B  »=  (gj  cosy  +  Vjy)sin^'B  (161) 
w 

Since  cos^B  t  sin^B  =  1,  equation  (161)  reduces  to 


cosB 


gLcosy  +  Vj^y 

i  **\3J 


(162) 


Multiplying  cosB  of  equation  (162)  times  tanB  of  equation 
(159)  the  following  results 


slnB 


Vjcosy 


e.n 
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or  solvina  for  o 


Kr^y, 

^sinB  (163) 

cos  y 


From  triaonometry 


sin  B 


(  1  —  cos“  B ) 


'/2 


(164) 


A 

so  that 


sin  B 


gr  cosy  +  v,y  2 ‘/^ 

( 1  -  (_t - L)^l 


w 


(165) 


Substituting  into  equation  (163)  and  rearranging 


.  0  'A 

-  (gj^cosy  +  v^y)") 
cos  y 


(166) 


Equations  (162)  and  (166)  provide  the  required  values  of 

B  and  0  to  sub.'titute  into  the  attitude  rate  equations, 

((147)  through  (149)'  .  Tiie  pitch  angular  acceleration, 

q,  required  in  equations  (1421  and  (143)  can  be  obtained 

from  time  differentiation  of  q  in  equation  (148) .  These 

*  .  ■ 

sin  B  a?id  o  become  negative  when 


1 


g,  COST  +  V /y  2 

- L_)  < 


Rr”j 


0. 
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equations  Include  the  true  normal  load  factor,  which  can 
be  approximated  satisfactorily  by  the  measured  value, 

INERTIAL  NAVIGATION  SYSTLMc- 

.-nertial  naviqation  system  from  v.'hich  are  obtained 
the  output  from  a  three-axj?  accelerometer  system  ana 
sufficient  orientation  angres  to  find  components  in  the 
wind-axis  system  has  a  very  inij^ortant  advantage  over  the 
accelerometer  installations  previously  described.  That  is, 
the  need  for  making  corrections  for  boom  bending,  upwash, 
dynamic  lag,  pitch  rate  or  for  induced  accelerations  in 
values  sensed  by  accelerometers  mounted  on  a  vane  are 
eliminated.  Since  the  inherent  accuracy  of  the  orientation 
angles  defining  the  position  of  the  three-axis  accelercmeter 
system  is  superior  to  conventional  angle  of  attack  measure¬ 
ments  and  the  uncertainties  associated  witr  above  correc¬ 
tions  are  eliminated,  ccmp>’.ted  load  factors  in  the  wind-axes 
from  an  inertial  navigatioi'.  system  can  be  expected  to  be 
significantly  more  accurate  than  those  from  the  accelerometer; 
systems  previously  described. 

Comparative  disadvantages  of  an  inertial  navigation 
sy.^tem  may  be;  increased  preflight  checkout  time,  poorer 
relicibility  with  increasea  maintenance,  and  higher  cost 
(presuming  that  the  test  aircraft  has  not  been  equipped 
with  an  operable  system) . 
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Because  of  the  variety  in  desiqn  of  inertial  nayi-:- 

(■ 

gation  systems,  the  orientation  of  the  accelerometers 
depends  on  the  specific  system  being  used.  Hence,  the 
axes  transformations  to  obtain  load  factors  in  wind-axes 
are  not  presented. 

RADAR  METHOD 

Equations  are  presented  which  may  be  used  with  data 
from  radar  (or  position  data  from  soiue  other  source  such 
as  Askania  cameras)  to  compute  load  factors  and  hence 
excess  thrust. 

Figura  7  shows  the  radar  system  located  with  respect 
to  the  geocentric  coordinate  system  in  terrm  of  oto  • 
centric  latitude,  6^^,  the  height  above  the  reference 
ellipsoid,  and  the  angular  deviation  of  the  x^.  axis  from 
true  north. 


Figure  7  Radar  System 
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The  origin  of  the  radar  system  in  geocentric  woordinates 


is  given  by 


Xg*  -  (r  +  hj.)cosSL 


(167) 


Yg'  ~  0 


(168) 


and 

V  “  -(r  +  (169) 

r,  repeating  equation  (5)  from  the  section  Geophysical 
Properties ,  is 

r  «  ro(. 99832172  +  .00167615  oos2Sl  +  .00000211  cos  4  (170) 

As  may  be  seen  by  figure  7,  rotations  of  (.;;j.+180®)  and 
6^  a.v.e  necessary  to  aline  the  radar  system  with  the  geocentric 
systsm.  These  rotations  result  in  a  rotation  matrix  given  by 


--  sinSLCOs 

—  cosSjJ 

M(  i  =* 

1 

rAni>j 

1 

-  cos^^ 

0 

(171) 

! 

1 

L—  cos5j^cos0j. 

cosS^sin^j. 

sinS^^  - 
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The  position  of  an  aircraft  in  geocentric  coordinates 
is  obtained  by  addition  of  the  rotated  radar  coordinates 
to  the  components  of  the  radar  origin  in  geocentric 
coordiijates : 


V 

“r 

^V] 

^  1  M( o,.  f  .7.  <’l)1 

>r 

+ 

•'’o' 

-'\r 

V 

-'V. 

(172) 


Since  X  y_' »  z  and  (M  ($  .6  )  ]  are  constants, 

6  6  6  IT  i>» 

the  transformations  of  aircraft  velocities  and  accelera¬ 
tions  are  given  by 


( M  ( Oj.  +  1  r 


(173) 


The  components  of  velocity  and  acceleration  ere  deter¬ 
mined  from  a  smoothing  ana  differentiating  process  of  the 
position  data. 

Since  the  y^,  system  is  attached  to  the  earth 
and  rotates  with  it  at  an  angular  rate  of  w»,  accelerations 
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due  to  are  gen^erated.  In  order  to  find  inertial  accel¬ 
erations  which  may  be  related  to  forces  acting  on  the  air¬ 
craft  and  to  lo.'id  factors,  equations  (174)  through  (176) 
are  used ; 


•• 


(174) 


•• 


2 

1  iJ 

m 


2x„  « 


a> 


(175; 


and 


•0 


(176) 


Next,  the  above  accelerations  are  transformed  to  local 
geocentric  su'.es.  This  is  accomplished  by  use  of  the  matrix 
found  in  ttie  section.  Coordinate  Systems  and  Transformations, 
and  results  in 


I 


f.  H  «« 

"g 

cos  AAlsIoSl  —  sin  AA^sin^  --  cos^l” 

Xp 

yg 

m 

sinAA^  —  cos.AAl  0 

•• 

L^'gJ 

_  -  cosAAlCosSl  •-  sIhAAlCOsSl  sinSj^  _ 

•• 

L  ZeJ 

(177) 
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AX-  is  the  difference  between  the  aircraft  longitude  and  the 
Xj 

longitude  of  the  radar  coordinate  origin  as  shown  in  figure 
3  of  that  section. 

Gravitational  attractions  as  computed  from  equation  (15) 
of  the  same  section  may  be  readily  added  to  the  accelerations 
from  equation  (177)  and  load  factors  computed: 


n  =  -  gx  )  (178) 

g  gf  ®  g 


Sr 


and 


g 


(179) 


(180) 


Transformation  into  wind  axes  but  at  zero  bank  angle 
is  accomplished  by  means  of  the  matrix  [M(a,Y)l  expressed 
in  equation  (14)  of  the  section.  Coordinate  Systems  and 
Transformations ; 


- 

"x 

*g 

=  [M(a,  y)] 

"y 

-  "v- 

L  Zg  J 

(181) 


The  angles  o  and  y  when  computed  from  radar  data  for 
use  in  eqmtion  (181)  require  a  knowledge  of  wind  speed 
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and  direction  as  a  function  of  time.  Velocities  in  local 
geocentric  axes  are  found  through  the  same  rotation  matrix 
that  was  used  with  accelerations  in  equation  (177): 


r 

% 

• 

«s 

1 _ 

cosAAlsIhSl 

sIhAAl 

cosAAlCosSl 


sinAALsinSL 
—  cosAAl 
-  sinAALCCsSL 


—  cosS^ 
0 

sin^T 


(182) 


Velocities  in  the  wind  axes  become 

=  x_  +  V,,cos^^^  (183) 

«  yg^  +  V^cos^  (184) 

=  K  (185) 

with  the  assumption  that  winds  are  in  the  local  geocentric 
system. 

Aircraft  velocity  with  respect  to  the  airmass  is 

(186) 

The  two  angles,  y,  and  o,  for  use  in  equation  (181) 
become 


'w 


and 


w 
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-1  ~ 


(ir-) 


Bank  angle  may  be  defined  as 


(188) 


—1 


(189) 


A  rotation  may  then  be  made  through  B  to  obtain  load  fac¬ 
tors  in  banked  flight.  Performing  the  transformation  the 
load  factors  become/  with  the  aid  of  equation  (189) 


(190) 


(191) 


,  2  2^ 
("v  .  +  ^7.  .  ) 


(192) 
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STANDARD  CLIMB  SCHEDULES 
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SUMMARY 


Various  standard  climb  schedules  are  considered  (e.g., 
a  segment  at  constant  calibrated  airspeed  followed  by  one  at 
constant  Mach  number) .  Equations  are  given  to  compute  the 
break  altitude  or  juncture  of  these  segments.  Climb  speed 
derivatives  from  which  acceleration  factors  may  be  computed 
are  derived  for  each  schedule. 
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SYMBOLS  UGLD  IN  THIS  SECTION 


Symbol 

Definition 

Units 

a 

speed  of  sound 

knots  or 
ft/sec 

acceleration  factor 

dimensionless 

^fE 

acceleration  factor  associated 
with  energy  height 

dimensionless 

reference  acceleration  due  to 
gravity 

ft/sec^ 

a 

L 

local  effective  acceleration  due 
tc  gravity 

ft/sec^ 

H 

c 

pressure  altitude 

ft 

«E 

energy  height 

ft 

M 

flight  Mach  number 

dimensionless 

P 

a 

ambient  pressure 

in,  Hg 

impact  pressure 

in.  Hg 

calibrated  airspeed 

knots 

equivalent  airspeed 

knots 

true  airspeed 

knots 

«a 

ambient  pressure  ratio 

dimensionless 

^a 

ambient  temperature  ratio 

dimensionless 

a 

air  density  ratio 

dimension.less 
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INTRODUCTION 


It  is  desirable  to  present  aircraft  performance  during  a 
continuous  climb  for  an  arbitrary  (standard)  climb  schedule; 
generally  the  one  which  the  pilot  was  attempting  to  follow. 
Since  climb  schedules  cannot  be  followed  precisely,  li’^rge 
deviations  in  rates  of  change  of  potential  and  kinetic  energies 
occur,  although  the  sum  of  the  two  can  usually  be  expected  to 
change  quite  uniformly.  In  order  to  find  rates  of  climb  along 
an  arbitrary  climb  schedule,  acceleration  factors  together  with 
climb  potential  (rate  of  climb  for  zero  acceleration  along  the 
flight  path)  are  used  as  described  in  the  section.  Standardiza¬ 
tion  of  Performance  Parameters.  The  information  required  to 
define  standard  schedules  and  to  compute  acceleration  factors 
is  presented  in  the  following  paragraphs. 


STANDARD  CLIfb  SChLDULLS 


In  standarcUzinr  .'‘ircraft  performance  uurino  continuous 
Clints,  it  is  desiral^lo  tc  r.aio  corrtetiens  to  an  arbitrary 
(stanuaru)  cliir.l'<  sci  eoulo;  generally  t';e  one  which  the  pilot 
was  attemptinn  to  follcv:.  Since  cliri''  schedules  cannot  be 
followed  precisely  and  clinb  potential  may  vary  rapidly  with 
speed,  large  corrections  may  result.  The  largest  corrections 
should  be  expected  vrhen  climbing  at  a  hinh  subsonic  speea 
(typical  of  best  clinib  speed  for  I.ioh  performance  aircraft) 
v'here  drag  increases  very  rapidly  wlien  tiacii  nmrbers  exceed 
those  for  best  clir.b. 

Four  standard  climb  schecules  arc  considered.  These  are; 

(1)  Initial  searent  at  constant  v/ith  the  remainder 
at  constant  P 

(2)  True  speed  specified  as  a  function  of  altitude  - 
for  any  arbitrary  schedule  vdth  altitude  montonic 
increasing  (test  and  standard  altitudes  taken  to 
be  ecjual) 

(3)  Initial  seement  at  constant  V  with  the  remainder 

e 

at  constant  f' 

(4)  True  speed  specified  as  a  function  of  altitude  - 
for  any  arl^itrary  schedule  with  altitude  periuitted 
to  decrease  but  with  Hg  monotonic  increasing  (test 
and  standard  energy  heights  take;,  to  be  equal) 


VM 


wwmpwJwcwjKUitcatafja  ' 


Constant  Mach  nuinber,  constant  V^,  and  constant  Vg  climbs  are 
not  listed  since  they  can  be  considered  as  special  cases  of  (1)  or 
(3) .  Equations  pertaining  to  the  standardization  of  climb  perfor¬ 
mance  along  each  of  the  above  schedules  are  derived  in  the  follow¬ 
ing  paragraphs. 

BREAK  ALTITUDE 

For  schedules  il)  or  (3)  it  is  convenient,  when  using  a 
digital  computer,  to  load  in  desired  standard  values  of  speed 
and  Mach  nun^ber  and  have  the  intersection  of  the  two  segments 
of  the  schedule  (break  altitude)  computed. 


Figure  1  Break  Altitude 
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To  do  this,  combinations  of  V’  less  than  or  m.ore  than 

and  subsonic  and  supersonic  numbers  must  be  considered. 

(Reference  figure  2.)  If  then,  of  course,  M  >  1.0 

so  that  the  combination  V  >  a,.,  and  M  <  1.0  cannot  exist. 

nreak  altitudes  for  the  three  remainino  combinations  of  V  and 

c 


I'i  are  derived  as  follov^s: 


M 

Figure  2  CembinatUns  of  and  M 

CASL  I  (V  <  a„T ,  M  <  1.0) 

Substituting  y  =  1.40  in  equation  (21)  of  the  section, 

Flight  Parameters  from  Sensed  Environment. 


P 


a 


[l  +  1 


(1) 


To  find  the  conm.ion  value  of  a  /P  and  hence  the  Ireak  altitude 

”C  a 


Vi~8 


the  above  equation  is  equated  to  the  corresponding  relation  con¬ 
taining  calibrated  airspeed 


-!-0.2( 


\^)213.5 

®SL 


-  1 


} 


which  follows  from  equation  (9)  of  the  section 
from  Sensed  Environment. 


1  |fl*0^(^)T'-l|  - 

^SL  J  > 

From  which 

fl  +  0.2(^)2P-  1 

;!  ®SL  ^ 

. . I 

(1  +  0.2M2'^‘^-1 


-  1 


5^  can  then  be  easily  converted  into  altitude, 
from  equation  (14;  in  the  section.  Atmospheric 
the  troposphere  in  the  U.S.  1962  atmosphere 


H 


c 


1/5.2559 

1  -8a 


6.87558x10“^ 


Vl-S 


(2) 

,  Flight  Parameters 

(3) 


(4) 

For  example. 
Environment,  for 


(5) 


CASE  II  (V^  <  M  >1.0) 

C  —  T’L 

Following  the  same  procedure,  from  equations  (2)  and  (8) 
from  the  section,  Flight  Parameters  from  Sensed  Environment. 


CASE  III  (V^  >  a,.,,  M  >  1.0) 

C  k.  ij 

As  before,  from  equations  (8)  and  (10)  of  the  same  section 


The  break  altitude  for  a  constant  -  constant  n  schedule 
can  be  found  more  readily.  Equation  (13) ,  also  from  the  same 
section,  ..lay  ■e  res  La  Led  as 

Vg  =  aM\/7'  (8) 

The  velocity  of  sound  in  a  perfect  gas  is  proportiondl  to  the 
square  root  of  tb.e  teiuperature  so  that 

-  /T  (9) 


Vl-lfl 


From  the  perfect  gas  lav:  6^  =  06^.  Substituting  this  relation¬ 
ship  and  equation  (9)  in  equation  (8) 

Vg  =>  ^gL  ^  (10) 

Solving  for  6 

3i 


V 


(11) 


The  break  altitude  may  then  be  computed,  as  for  the  con¬ 
stant  V  -  constant  M  schedule,  from  equation  (5)  or  other 
appropriate  6  -  K  relation  for  a  given  atmospheric  layer. 

Oi 

Equation  (11)  is  valid  ;for  all  values  of  M  and  Vg. 

AIRSPEED-MACn  NUMBER  RELATIOHShIPS 

The  airspeed  -  Mach  number  relationships  which  are  needed 
for  climbs  are  presented  belov/.  For  exaraple,  v;hen  using  sched¬ 
ule  (1) ,  Mach  numbers  at  the  test  altitudes  along  the  constant 
portion  of  the  climb  are  desired;  conversely,  when  a  con¬ 
stant  Mach  number  is  specified  the  corresponding  calibrated 
airspeeds  are  desired. 

The  equations  that  were  used  to  determine  break  altitudes 

for  the  same  combinations  of  and  M  are  used  to  calculate 

c 

both  V  and  M. 
c 
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CALIBRATED  AIRSPEED  AND  MACH  NUMBER  (V  <  M  <  1.0) 

C  —  SL 

V  is,  solving  equation  (3) 


(12) 


and  M  is 


X  ,2,3.5 


M 


v's 


’({ _ "SL 

V  5, 

—  +  ij  ij 

(13) 


CALIBRATED  AIRSPEED  AND  MACH  NUMBER  (V^  <  agj^,  tl  >  1.0) 
Solving  equation  (6)  is 


/ 

r  ,  ..7  1 

2/7  \ 

S„ 

166.921  M'  _i 

+  1 

"a 

(7m2-  1)^-^ 

(14) 


and  M  is 


(15) 


In  this  case  an  iterative  solurion  for  M  must  be  made  since  equa¬ 
tion  (7)  cannot  be  solved  explicitly. 
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M  >  1.0) 


CALIBRATED  AIRSPEED  AND  MACH  NUMBER  (V^  > 
Solving,,  eqyatipii  {7,)  ;for 


166.921  M‘ 

(Tpri)^ 


(16)  .  _) 


In  this  case  also,  must  be  found  throuqj^i  an  iterative  procedure. 
M  from  equation  (7)  is 


(17) 


here  again  an  iterative  procedure  is  needed  tofji.o  li.  ' 

V  may  be  found  for  the  remaining  schedules  by, first  con- 
verting  or  V.  to  V  and  then  using  equation  (12)  ,  (14)  , 
or  (16)  as  appropriate. 

CLIIIE  .SPEED  DEPIVA/riVL,  dV^/dli^ 

The  acceleration  factors  discussed  in  the  section,  Standardiza¬ 
tion  of  Performance  Parameters,  and  defined  by  the  equations 


Af 


gf  ^  gf 


and 


AfE 


£l 

«L 


(  1 


!tl!L 

«r 


(18a) 


(18h) 


are  of  importance  in  establishina  climb  uer-^ormance .  Witti  mem 
rate  of  climb  along  tlie  standard  clini>  schedule  may  be  found 

from  standardized  excess  thrust.  Values  of  dv,j./dl  ^  for  tiie 
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/ 


various  combinations  o£  and  H  are  derived  in  the  following 
paragraphs  for  a  constant  V^. 

CLIMB  SPEED  DERIVATIVE,  dV./dH  ,  (V  <  a^. /  M  <  1*0) 

From  differentiation  of  the  relation  V^  =  aM 


dH,  “  ‘  dH/  'dH, 


Since  constant  with  respect  to  the  differentiation 

indicated,  the  term 


r  V  9i3.5 


in  equation  (13)  may  be  replaced  to  give 


Differentiation  with  respect  to  produces 


5kA 


7  5a^M(^  + 

®  S3 
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Equation  (21)  might  be  used,  substituting  it  in  equation  (19); 
a  more  simple  form  may  be  derived,  however.  From  the  definition 
of 


(22) 

5  S„ 


Substituting  in  equation  (21) 


5[(1.0.2m2)C-^-i]^^ 
TS^Md  +0.2M^)'‘’^ 


(2?^ 


But 

Lifa  .  ddnSg) 

so  that  substitution  in  equation  (19)  produces 


dHc 


9  3i  5 

M  5a[(l  +  0-2Md  -1]  d(ln5a) 

7M(1 +0.2M^)'^’'’ 


(24) 


The  derivatives  da/dH  and  d()ln  6  )/dH  are  presented  in  subse- 

O  4  w 

quent  paragraphs. 
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CLIMB  SPEED  DERIVATIVE,  dV./dli  ,  (V  <  ,  M  >1.0) 

u  C  O  ojj 


Equation  (l,i)  may  be  rewritten  as 


l'^6.921M^  =  (^  +  1){7M^  -  1)^*^ 


(25) 


Differentiating  and  collecting  terms 


dM 

dH 


. [7 X  166.921  M^-  35M(|l  +  1)(7M^-  1)^'^] 


,  ^2.5  ddngg) 


- i(7M^-l) 


dH, 


(26) 


Substituting  for  from  equation  (25)  and  rearranging 


dM 

dH, 


1^166.921  M^- 


(7M^-1)^*^ 


l7M^--l)i.i"i9i 

dH,, 


(27) 


lie8.45M^(2M2  -  1) 


Substituting,  as  before,  in  equation  (19) 


11.  Mil. 

dH„  dH, 


_  a[l66.921M'^-(7M^-l)^*^](7M^-  1)  d(lt^) 


1168.45 M^(2M"  -  1) 


dH, 


(28) 


Vl-M 


CLIMB  SPEED  DERIVATIVE,  dV^/dH^,  (V^  >  agj^,  M  >  1.0) 
Rewriting  equation  (17) 


166.921M‘  -  t 


1  ( 


,M[7&)^-  n“ 

aSL 


1  UiK7M^-l)^^ 


where 


166.921M^  .  (^  +  1X7m2-1)^ 
Sa 


166.92U^)^ 

_ !sl_ 

V  2 

[7(^)2-  11 
®SL 


Since  the  form  of  equation  (30)  is  the  same  as  that  of  equa¬ 
tion  (25)  the  same  final  equation  for  dV^/dH^,  will  also  apply 
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i 

t 


in  this  case.  Repeating  equation  (28) 


a[  166.921  -  ( 7 ~  1  )^^1{ 7M^  -  1 )  ddnSg) 

1168.45M®(2M2  -  1)  dH^ 


(32) 


CLIMB  SPEED  DERIVATIVE,  dV./dli^,  (V^  not  constant) 

c  C  w 

Values  of  dV. /dh  during  a  climb  at  constant  V  (as  during 

w  w  w 

schedule  (3))  may  be  found  from  differentiation  of  =  v^//3- 
to  obtain 


liVt  Ve  d.  (33) 

■IMc  ' 

Tor  the  cases  where  V^  xs  specified  as  a  function  of 
(schedules  (2)  and  (4))  dV^/dH^  or  dV^/dKg  must  be  found  by 
numerical  differentiation  since  they  cannot  be  calculated  by 
analytic  means. 

Vvhen  Mach  niimber  is  maintained  constant  equation  (19) 
reduces  to 

<34) 

dHe  dH^ 

and  values  of  da/dli  ,  as  for  constant  V_  climbs,  are  found 

c  c 

from  the  following  equations. 
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0.0100234 


(39) 


d  a 
dH„ 


[216.63  +  0.0  0  0  30  <l8(Hj.- Hj,)] 


for  65, 61b. 80  <  ii^  <  104,986.88  feet  and 


da  ^  0.0280641 _ 

[228.65  +  0.0008534(H^-Hb)]^ 

for  104,986.88  <  L  <  154,199.48  feet 

^  Mi. 


(40) 


P’^rseURL  RATIO  ULRIVATIVL,  d(J,n  6  ) /dh 

_  3.  C 

The  follov'inc  equations  are  aiven  as  an  aid  in  coitputing 
dV^/oii^  and  are  to  be  used  in  equations  (24),  (28),  and  (32). 


'  nation  (16)  from  the  section.  Atmospheric  Environment, 
may  .  ?  rewritten  as 


(Pa)b 

Pg 

®SL 


1  +  w) 

■  ■ 


(41) 


or  taking  the  logarithm 


1„5,  .  -.SL|ni^ri,i!l(H  - 


RL 


M 


'SL 


(TJ 


Hb) 


a'b 


(42) 
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Differentiating  v;ith  respect  to 


d(lnS») 
- ^ 

dHc 


_!sl_ 


(TJv 


(Hc-Hb) 


(43) 


when  L,,  ^  0. 

M 

Taking  logs  of  equation  (17)  from  the  section,  Atmospheric. 
Environment 


In  8a 


In 


(PJ; 


a'b 


^SL  ^  ^ 


(44) 


Differentiating  v;ith  respect  to  li^ 


gS'L 

dHc  "  “R(Ta)b 


(45) 


when  L„  =  0. 

M 

Substituting  constants  appropriate  to  the  U.S.  1962  atiuosphere 


d(InSa) 


-  0.361374x10"'^ 


_ 1 _ 

1-6.875586x10“^  He 


(46) 


in  units  of  feet”^  for  -16,404.20  <  <  36,nRQ.?4  feet 


ddnSg) 

dHc 


-0.480637x10“'^ 


for  36,089.24  <  ii  < 


65,616.80 


feet 


(47) 
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=  -  0.480637x10'^  - 1 - (48) 

[l  +  1.40688  xlO~^(H^-  65,616.80) 

for  65,616.80  <  1.  <  104,986.88  feet  and 

—  c  ~ 


-  -  0.455412x10"'^  - -  (49) 

flHc  [l  +3.73252xl0"^(Hc-  104,986.88 )_ 

for  104,986.88  <  <  154,199.48  feet 

DENSITY  RATIO  DERIVATIVE,  do/dH„ 

The  following  equations  are  listed  to  provide  assistance 
in  computing  dV./dH  when  V  is  constant,  as  in  equation  33. 

t  C  0 

After  differentiating  equations  (19)  and  (20)  from  the  section, 
Atmospheric  Environment,  we  have ,  after  substitution  of 
appropriate  constants  from  the  1962  atmosphere 

~  -  -2.92618x10“^  1  -  6.87558xl0“®H„)^'^^^®  (50' 

dH,  C 

in  units  of  feet  ^  for  -16,404.20  <  H  ^  36,089.24  feet 

c 

_1.4278.';x  10“®exp[ -4, 80637x10“^  H„  -  36089.24)]  (51) 

for  36,089.24  <  65,616.80  feet 
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*  -3  55517xl0“®[l  +  1.40688xl0"^(Hg-  65616.8)]“^®'  (52) 

IHc 

for  65,616.8  5  5  104,986.88  feet  and 


for  104986.88 


-5.31944x10  '^[i  +  3.73252 X 10 “®(  Hg  -  104986.88)] 

5  H  -  154,199.48  feet, 

c 


(53) 
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SECTION  VII 

STANDARDIZATION  OF 
EXCESS  THRUST 
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SUMMARY 


Standard  excess  thrust  is  found  by  extrapolating  to  a  set 
of  standard  conditions  using  a  Taylor  series  expansion.  From 
a  functional  statement  defining  excess  thrust  and  the  equations 
of  motion,  partial  derivatives  of  excess  thrust  with  respect 
to  each  of  the  independent  variables  have  been  found.  The 
terms  are  collected  into  an  equation  which  may  be  used  for 
both  climbs  and  level  accelerations.  This  has  been  done  for 
first  order  terms,  which,  it  is  expected,  will  generally  be 
sufficient;  second  order  terms  have  been  derived  so  that  they 
may  be  included  if  their  magnitude  warrants  it.  The  standardi¬ 
zation  equation  cannot  be  solved  directly,  and  an  iteration 
procedure  is  needed.  A  description  of  it  is  included. 
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SYMBOLS  USED  IN  THIS  SECTION 


Symbol 

Definition 

Units 

"•D 

airplane  total  drag  coefficient 

dimensionless 

airplane  lift  coefficient 

dimensionless 

eg 

center  of  gravity 

pet  MAC 

engine  ram  drag 

lb 

^ex 

excess  thrust 

lb 

gross  thrust 

lb 

net  thrust 

lb 

g 

acceleration  due  to  gravity 

ft/sec^ 

reference  acceleration  due  to 
gravity 

ft/sec^ 

h 

tapeline  altitude 

ft 

pressure  altitude 

ft 

ip 

thrust  angle  of  incidence 

rad 

M 

flight  Mach  Number 

dimensionless 

^*z 

w 

load  factor  along  negative 
wind  z-axis 

dimensionless 

p 

a 

ambient  pressure 

in.  hg 

q 

dynamic  pressure 

Ib/ft^ 

r 

local  radius  of  the  earth 

ft 

s 

wing  area 

ft^ 

t 

time 

sec 

T 

a 

ambient  temperature 

deg  K 

true  airspeed 

knots 

w 

airplane  gross  weight 

lb 

VIM 


emgle  of  attack 

» 

flightpath  climb  angle  measured 
from  the  geocentric  horizontal 
plane 

flightpath  climb  angle  measured 
from  the  geodetic  horizontal 
plane 

aircraft  geocentric  latitude 
flightpath  heading  angle 
angular  velocity  of  the  earth 


rad/sec 


standard  day  conditions 

test  day  conditions 

dot  denotes  first  derivative  of  a 
quantity  with  respect  to  time 
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INTRODUCTION 


The  test  excess  thrust,  computed  by  one  of  the  methods  out¬ 
lined  in  the  section,  Determination  of  Excess  Thrust,  is  extrapo¬ 
lated  to  a  set  of  standard  conditions.  In  the  following  analysis, 
excess  thrust  is  taken  to  be  sensitive  to  the  following  parameters; 
(1)  temperature,  (2)  weight,  (3)  Mach  number,  (4)  pressure  or  alti¬ 
tude,  (5)  normal  load  factor,  (6)  power  setting,  and  (7)  center  of 
gravity.  Of  these  variables  the  first  five  obviously  have  a  bear¬ 
ing  on  thrust  and/or  drag  and  hence  on  excess  thrust.  Adjustments 
to  excess  thrust  for  changes  in  power  setting  would  not  be  made  in 
most  cases  but  might  be  desirable  for  a  variety  of  reasons;  for 
example,  a  correction  may  be  needed  to  account  for  a  power  lever 
v;hich  has  been  improperly  riaged.  Inlet  total  pressure  ratio 
was  not  included  as  an  independent  variable  since  it  was  consid¬ 
ered  to  be  a  function  solely  of  Mach  nuniber,  which  has  been 
treated  as  an  independent  variable.  This  assur:.ption  can,  in 
general,  be  made  quite  satisfactorily.  As  speed  is  increased 
and  a  Mach  number  of  perhaps  three  is  reached,  total  pressure 
ratio  may  become  a  function  of  other  variables  as  v/ell  as  Mach 
number.  In  this  event,  appropriate  functional  statements, 
depending  on  the  particular  installation,  need  to  be  written 
and  corrections  to  excess  thrust  made  for  variations  in  total 
pressure  ratio. 

Standardization  of  excess  tfirust  is  accomplished  by  first 
expanding  excess  thrust  in  a  Taylor  series  about  the  test  day 
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point  through  the  first  order  terms.  Then  the  development  of 
expressions  for  the  partial  derivatives  of  excess  thrust  with 
respect  to  each  of  the  variable  test  conditions  is  presented. 
Next,  the  equation  for  the  standard  excess  thrust  in  terms  of 
flight  test  data  is  developed.  Finally,  a  mere  precise  equa¬ 
tion  for  standard  excess  thrust  is  derived  by  following  the 
above  steps  but  carrying  the  Taylor  series  expansion  through 
the  second  order  terms. 

TAYLOR  SERIES  EXPANSION  OF  EXCIiSS  TERUST-FIRST  ORDER  TERMS 

The  two  equations  of  motion  resulting  from  force  balances 
along  and  normal  to  the  airplane's  velocity  vector  are  basic 
to  derivatives  in  the  following  paragraphs.  Excess  thrust, 
from  the  longitudinal  equation  of  motion  is* 

^ex  “  F,,cos(a+  ip)  -  qSC|3  (1) 

and  the  equation  of  motion  normal  to  the  flightpath  is 

qSCL  -  n^W  -  FgSin(a  +  ip)  (2) 

Equation  (1)  is  examined  to  determine  a  functional  rela¬ 
tionship  for  excess  thrust.  It  is  well  known  that  net  thrust 
depends  on  ambient  temperature,  Mach  number,  and  pressure  (alti¬ 
tude)  .  Power  setting,  tt,  is  also  included  so  that  changes  in 
*The  excess  thrust  is,  properly,  F  cos(a+i  )-F  -qSCp, .  This  is 

r  6  u 

quite  a  satisfactory  approximation,  hov/ever,  and  <=!imDlifies  the 
following  derivations. 
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power  lever  angle  or  other  measure  of  power  setting  can  be  made. 
Hence,  =  f  M,  P^,  tt)  .  Both  and  may  be  expressed 

as  functions  of  a,  M,  and  eg  (C^=f{a,  M,  eg)  and  C^=f(a,  M,  eg)) 
with  the  center  of  gravity  included  to  account  for  changes  in 
trim  drag.  It  is  more  convenient,  however,  to  express  drag 
coefficient  as  Cj^=f (Cj  ,  M,  eg).  From  these  relationships 

together  with  q=0.7P^M^  and  C|_  =  — from  equation 

Q  S 

(2) ,  excess  thrust  may  be  completely  defined  by  the  functional 
statement 

^ex  =  ^a'  "z'  (3) 

A  Taylor  serie;^  expansion  of  the  excess  thrust  about  the  test 
day  point  results  in  the  expression  for  standard  excess  thrust: 


Magnitudes  of  the  higher  order  terms  are  less  than  those  of 
the  first  order  terms  but  will,  at  some  test  conditions,  be 
significant.  The  partial  derivatives  in  equation  (4)  are 
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evaluated  first,  and  then  attention  is  given  to  the  second 
order  terms. 

VARIATION  OF  EXCESS  THRUST  WITH  TEST  CONDITIONS 

The  partial  derivatives  in  equation  (4)  are  found  by 
differentiating  equation  (1)  v;ith  respect  to  each  of  the 
independent  variables.  It  is  important  to  note  that,  in  each 
case,  the  normal  equation  of  motion  (eauation  (2))  provides 
a  constraint  on  the  equation  for  excess  thrust.  The  angle  of 
attack  must  vary  in  such  a  manner  that  both  equations  are  sat¬ 
isfied  simultaneously. 

Temperature 

To  determine  how  excess  thrust  varies  with  temperature, 
the  relationship  q=0.7?  is  used  and  the  oartial  derivative 

ci 

of  equation  (1)  taken  with  respect  to  temperature: 


^  -  cos(a  +  iF)^-  Fnsin(a  +  lp)||.  -  0.7 


The  reason  thrust  varies  with  airtbient  temperature  is  well 
known;  however,  the  dependence  of  angle  of  attack  and  drag 
coefficient  on  temperature  requires  an  explanation. 

Figure  1  shows  a  component  of  thrust  in  the  direction 
of  the  lift  vector.  If  the  thrust  is  varied,  then  the  angle 
of  attack  must  vary  to  adjust  the  lift  to  the  new  condition. 
The  changes  in  angle  of  attack  and  lift  cause  a  change  in  the 
induced  drag  similar  to  that  caused  by  a  weight  change. 
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Chain  rule  differentiation  of  a  =  M,  eg)  with  M  and  ca 

held  constant  produces 


=  da 


dCi 


dC^dJ^ 


and  similarly 
(9Cq 


(6) 


(7) 


Equation 


df. 


ex 


dJ, 


COS(a  +  ip  ) 


dFf 

W, 


da 

F„sin(a  +  ip  )  — 


'  +  qS-~ 
L 


(?C| 

FT. 


(8) 


The  partial  of  the  lift  coefficient  with  respect  to  temperature 
is  found  by  differentiating  the  equation  of  motion  normal  to 
the  flight  path: 


qS 


-sin(a+  ip)^  ~  FnCOs(a+  ip)~ 


(9) 


Again  using  the  transformation  expressed  by  equation  (6)  : 


dC  |_ 

K 


sin(a  +  ip)-gr^ 
_ iLa- 


da 


qS  +  F^cosCa  +  ip)r7r 

c/V/  1 


(10) 


The  equation  for  the  variation  of  excess  thrust  with  temperature 
is  found  by  substituting  equation  (10)  in  equation  (8), 


Vil-ll 


(11) 


COS(a  +  ip 


sin(  a  + 


<9Cr\ 

(?Fg  F„sin(a+ip)  +  qS^ 
*^^3  Fg  cos(  a  +  ip)  +  qS^^ 


Weight 

The  partial  derivative  of  excess  thrust  with  respect  to  weight 
using  equation  (1)  as  the  expression  for  excess  thrust  is 


-  F^  sin ( a  +  ip 


.  da 

'm 


(12) 


Again  transforming  variables  through  chain  rule  differentiation 
with  M  and  co  constant; 


da  da  \_ 

"iW  “  5Cl<9W 


and 


<9Cp  <9Cq  (9Cl 

W  ”  d^W 


(13) 


(14) 


Substituting  these  transformations  into  equation  (12)  and 
rearranging; 


fex 

aw 


-fil.sin(a  +  ip)  +  qS^ 
aCL 


ac  |_ 
W 


(15) 


The  partial  of  the  lift  coefficient  with  respect  to  weight  is 
found  by  differentiating  equation  (2). 


qS 


aCi 

W 


Fg  cos(  a  +  ip  ) 


da 


(16) 
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Substituting  equation  (13)  in  equation  (16)  we  have 


5C|_ 

W  ‘ 


cos{a  +  ip  )  +  qS 


The  general  expression  for  the  variation  of  excess  thrust 
with  weight  from  equations  (17)  and  (15)  is 


F„sm(atlF)  +  qS  — 

- — —  n_ 

^(9C|  ^ 

Fg  COS(a  4-  ip  )  + 


Mach  Number 

Inserting  the  relationship  q=0.7P  in  equation  (1) 

Ci 

and  differentiating,  we  have 


^  c  •  /  .  ZqSCq 

_"oos(a»,F)-  F„s,n(.  Vsir- 


Holding  eg  constant,  is  a  function  of  M  and  a,  but  a  is 
dependent  on  M  so  that 


Cn  *  f 


[m,  a(M)j 


Differentiation  of  equation  (20)  with  respect  to  M  yields 


5Cnl  5Cnl  da 


<9M  dm  du  IM  (9M 


which  reduces  to 


(9Cq  5Cp 

5M  5M  a  da  iMcJM 
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Pressure 


Taking  the  derivative  of  equation  (1)  with  respect  to 
pressure 


<?Pa 


dFf, 

COs(a+ 


-  Ff,sin(a  +  ip  ) 


da  qSCp 

"Pa  ^Pa 


(27) 


Substituting  the  transformation  of  variables  produced  by  the 
same  method  as  previously  used 


da  da 

and 

5Cn  5C[)  ^Ci 

^  ■  5c?ip7 


results  in 


(28) 


(29) 


df. 


ex 


dP. 


dfn 

COS(a  +  ip 

dr- 


gSCp 

p. 


-  PpSi 


a  da 


sin(a  +  ip)'^  + 

0\j  L 


qS 


aCi 


dCj  dP,_ 


(30) 


The  partial  derivative  of  lift  coefficient  with  respect  to 
pressure  is  evaluated  from  the  equation  of  motion  normal  to 
the  flightpath.  Differentiating  equation  (2) 


q 


qSC, 


<5Pa  .  .  .a  r-  .  ,  ^da 

-  ~sinva  +  ip;-rgCOS\a  +  ip/  — 


/  *1  \ 
V  JX/ 
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Using  the  transformation 


da  da 

Wl  ■  acLi>p, 


(32) 


equation  (31)  becomes,  after  rearranging 
dCi  ^sin(a  +  iF)+ 


dP. 


.da 


Fg  cos{a  +■  +  ‘1^ 


(33) 


Substituting  equation  (33)  in  equation  (30)  we  have  the 
general  expression  for  the  change  in  thrust  v;ith  pressure. 


dP. 


ex 


dP, 


cos  (  a  +  ic  )—:r-  ~ 


qSCp 


«Fa  Fa 


F_sin(a  +  ic  )  +  qSr — 

"  da 


Fg  COS  (a  +  ip  )  +  qS 


da 


dPn  ,  ,  .  .  qSCL 

^sin(a+  if)+  -y~ 


(34) 


It  may  be  more  convenient  to  make  the  trcUiti  format  ion 


dP, 


£!n^Jlc 

^Hc  aPg 


(35) 


an 

and  substitute  in  equation  (34)  with  being  evaluated 

from  the  appropriate  altitude  -  pressure  relationships. 
For  example,  differentiating 


Pg  =  Pg^^d  -  6.87558x10'^ 


(36) 
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taken  from  the  U.S.  Standard  Atmosphere,  1962,  for  altitudes 
belov/  the  tropopause 


aHj. _ 13.07b3 _ 

^**3  (1  -  6.87552x10'^ 


(37) 


in  the  units  of  feet/ lb  per  square  foot. 

Load  Factor,  n 

z 

Differentiating  equation  (1)  with  resnect  to  normal  load 
factor  yields 


dn^ 


-  F^sin(a  +  ip 


Proceeding  as  before 

da  da 

dC  L 

and 

3Cq  5Cq 

5Cl  dn^ 


equation  (38)  becomes 


(38) 


(39) 


(40) 


On, 


-  sin(a  +  Ip  )— - -  -  qS  ——  — - 


(41) 


9  Cl 

The  tern;  is  evaluated  by  takinc  the  partial  derivati'  c  of 
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equation  (2) 


qS 


dC  I 
(9n, 


W  -  F  cos(a  +  Ip ) 


da 

5n, 


(42) 


FroK  which 


aC] 

5n, 


W 


da 


cosla  f  ip)  +  qS 


(43) 


Substituting  equation  (43)  in  equation  (41) 


5Cn”^ 

sin(a  +  ip  )  +  qS  — = 
da 


r.  5C| 

Fg  cos(a  +  ip}+qS^j 


W 


(44) 


Power  Setting 

The  expression  for  the  partial  of  excess  thrust  with  respect 
to  power  setting  is  developed  by  taking  the  partial  of  equation 

(1)  with  respect  to  tt/ 


dn 


ex 


.<3F 


V  I  p 

coi;(a  +  !p  -  F„  sin(a  +  ip) 

on  ^ 


da 

dn 


-qS^ 

dn 


(45) 


Expressions  for  the  follov/ing  chance  of  variables  may  be  written 


following  the  methods  used  in  the  preceding  derivations: 


da  da 

dn  dC^dn 


(46) 
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3Cp  5C[j  3C|_ 

dn  (9C|_  dn 


equation  (45)  becomes 


3Fav  i3Fn 

“•—  »  cos( a  +  Ip  ) 

diT  ^  dn 


The  variation  of  lift  coefficient  with  power  setting  is  deter¬ 
mined  by  taking  the  partial  derivative  of  equation  (2)  with 
respect  to  power  setting. 


-~sin(a+lp)-  FgCos(a+ip)^ 


Substituting  equation  (46)  and  rearranging 


'sin(a  +  Ip ) 


Fg  cos(a  +  Ip )  +-qS 


The  variation  of  excess  thrust  with  power  setting  becomes, 
after  substitution  of  equation  (50)  in  equation  (48) 

5Cn 

,  .  ,«Fn  +  jf,  .  ,  .  , 

IT  ■  f“<«-'fv  + - 5c:  17^'"'“^''^’ 

FgCos(a  +  Ip)  -t-  qS  — = 
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Trim  Drag 

Differentiating  eauation  (1)  with  respect  to  eg  yields 


ex 

(}cq 


r-  .  ,  .  V  (9a  Q 

-  F  sin(a  +  Ip)- - qS— y 

"  deg  deg 


(52) 


As  previously  stated  =  f(a,  M,  eg).  Since  trim  drag  cor¬ 
rections  are  made  at  constant  liach  nun*ber  the  functional  rela¬ 
tionship 


C 


D 


(53) 


may  .'e  set  dov;n.  The  partial  derivative  of  with  respect  to  eg 
then  becomes 


deg  deg  'a  da  'eg  deg 


(54) 


Substituting  in  equation  (52) 


Jcg 


-  FpSin(a+ 'iU) 
''  ^  deg  deg  a  Oa  'eg. /eg 


dr. 


dc 


(55) 


From  differentiation  of  equation  (2) 


,  dC| 


dCi 


qS(“-  + 

deg  'a  da  'eg  dcg 


-) 


Fg  COS(a  +  ip) 


deg 


Solving 


for 


9a 
9  CO 


(56) 
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^a. 

deg 


dCit 

Ff,  cos(a  +  ie)  +  qS— ^1 
«  ^  da  leg 


(57) 


Substituting  equation  (57)  in  equation  (55) 


fex 

deg 


F„sl„(a*iF>+<.s'|e| 

- - - qS  — t 

dCi  deg 

eos(a  +  ip)  +  q 
9  da 


-  qs^ 

a  deg  a 
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(58) 


An  examination  of  the  preceding  derivations  shows  that  the 


term 


dCr 


F-sin(a  +  ip)  +  qS— ,, 

^  da  'M,eg 


dCi  I 

F^  eos(a  +  ip)  +  qSr-= 

8  ^  da  !M,eg 

occurs  in  each  of  the  equations  for  the  partial  derivatives. 
Setting  this  term  equal  to  A  and  summing  up  the  corrections  as 
indicated  by  equation  (2)  we  have 


f'exs  “  ^ext  +  [cos(a  +  ip)  ^  +  Asin(a  +  ip)  j^l{  -  Tg^) 

3  3*^ 

dFg  ,  .  ^  ZC^qS 

^cos(a+iF) - y - qS~ 


~  An^eWg- Wt)  +  ^ 

I  V  II 

2  CLqS 


5F, 


+  — 9sin(a  +  iF)|^(  Mg  ~  M^) 


]} 


{•' 


.  ,5F„3Hj  qSCp 


iFg  aHg 


SP, 


sin(a  +  ip) 


'  a 
A  Wt( 


-  >(Pa  -  F3,)  +  [cOS(a  +  ip)-^  +  sin(a  +  ip^(jF5  -  n, 

r  *.  L  an  djT  J  S  t 

Hz  )  +  [AqS^I  - 


■  <^Co 
*  deg 


( egg  -  eg,) 


(59) 
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Equation  (59)  can  be  simplified  by  considering  the  total 
change  in  net  thrust  caused  by  changes  in  temperature,  Mach  num¬ 
ber,  and  pressure.  The  total  change  in  net  thrust  is 


F 


(60) 


or 


"t 


- 

JTj  as 


T-  )  + 


il(Ms  -  Mt)  +  T:r^(  P, 


dM 


(5P. 


‘s. 


^a,> 


(61) 


The  change  in  net  thrust  may  be  determined  from  the  engine  manu¬ 
facturer's  specifications,  from  computer  programs  constructed  by 
the  engine  manufacturer,  or  from  fliglit  test  data. 

Substitution  of  equation  (61)  in  equation  (59)  results 


in 


ex. 


Fext  +  cos(a+  ip)(F„^-Fn^)  +  Asin(  a  +  ijr)(  Fg^  -  Fg^) 

r 

-  An2(Ws  -  Wj)  + 


.  ,S,A§ 

m  M  «9M 


(  Mg  -  Nip  +  (  -  ^4^  +  A  544^  )(  P;,  "  P:,.) 


P,  ' '  ■ 


-  A  W  ( tm  _  n-  )  + 

4s 


A  qS 


.  dCi 


deg 


-  qS— ^ 
a  deg  Kz 


(egg  -  cgp 


(62) 


TAYLOR  SERIES  EXPANSIOIT  OF  EXCESS  TliRUST-SLCOND  ORDER  TERMS 
Restating  equation  (4)  but  carrying  the  exoansion  through 
the  second  order  terms,  we  have,  omitting  pov/er  settina 


V 11-22 


ex. 


F  +  ^ 


*^*^ex/ B  PI  X  ‘^^ex/  X  <^f^py, 

+  ^3  “  P3  J  n,  -  n,  )  +  — ^(  eg,  ~  cat) 

5Pa  h  ®t  dn^  H  ^t  aeg  '-yt' 


dh, 


dJ 


?'\- V  *  v‘*--*t> 

j  3 


2  2 

+  2^=%.(Ta  -T.)(Ms-Mt)  +  2^^ 

dJ^dn  h  h  ^  '  aT^ap 


hF 

r=r-||-(Ta  -  T.  )(  P.  -  P.  ) 

2  2 

2  2 

-^<W.-Wt)'.ai^(«,-W,)(M3-M,) 

2  2 

*  V  -  <-5^<»s-W,><"z3-n.,) 

2p  2 

2  2 

*  V  *  'sfe;“s-"t""z3-"z,) 

*  ’’as  -  -  "V  ^  -  ”31’'  -  eg,) 

!>X(„  .  ,2 


— ~  ■*■  }(  cg^  -  eg*) 

dn^2  ^s  H  dn^deg  ^s  ^ 


d% 


5cg' 


■|^(cgs  -  cgj)^ 
r 


(63) 


Before  evaluating  the  parti.al  derivatives ,  it  should  be 
noted  that  the  total  change  in  net  thrust,  as  identified  in 
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equations  (<50)  and  (61)  tor  first  order  terms,  should 
also  include  the  follov/imi  terras 


eh 


!(T. 


2 


,)Pa2  5 

e^F, 


\  ~  '"as  Pat’'  and  -  MtXP^^  -  Pg^) 


Evaluation  of  the  partial  derivatives  using  equations  (1) 

and  (2)  would  become  quite  involved.  Since  the  magnitudes 
of  the  second  order  terns  are  less  than  those  of  the  first 
order  terms ,  the  equations  of  raotion  can  be  reduced  to 


P^ex  ~  Pn  “  *^D 


(64) 


and 


qSCL  =  n^W  (65) 

and  the  partial  derivatives  found  much  raore  easily  and  with 
little  loss  in  accuracy.  U'ith  these  simplifications  the  second 
partial  derivatives  are  (derivations  have  not  been  included)  as 
listed  belovj. 

list  .  ^ 
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^ex  ^  +  _J^LJL(  -  i5i=)  + 

(?M^  aM^  LM'JM  Cl  m2  aCLM  5m2  Cl  “  M^j 


<9MaPa 


Au  aife  +f!£fi(i£u^ 

msPi  p,l3M  Cl.  5Ci_2  m  . 

sc  ■ 

+  .7MS  -2Cn  +  2^  Cl  +  M - ^  Cl 

°  ^Cl  '-  aw  Cl 


NOTE:  This  correction  can  be  icnored  for  cliirJDs  and  level 

accelerations  where  either  AM  or  AP  =0.  It  will  take  on  some 

a 

value  when  both  M  and  P  are  correctea  as  during  an  ' optimum 
climb. " 


.2p  -P  ,2p 

^  ^ex  ^  _  i£L  1£d 

(9M_(?n2  <9M  Cl  ^ 


*SX. 

5M^cg 


'^^''ex  .  qSCL^^^Cc 

7  “  -  +  2  2 

‘’Pa  PPa  Pa  '’^L 

^  .  Cl«  '’^Cp 

^Pgf^n^  Pg  SC^ 
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2^ 

^  )  ■ 

^  p« 

-  qS 

1  I  Cl 

^Pg^cg 

Pa  Icl  Pa  <3Cl 

^  Pex 

a^Cp 

~~~r  “ 

an/ 

aC|_^ 

acp 

2 

a(— fi) 

^  Pex  . 

-w 

an^aeg 

aeg 

dh. 


deg 


-qS 


Standard  excess  thrust  may  be  computed  by  substitutina  the 
partial  derivatives  listed  above  in  aquation  (63)  to  find 
the  second  order  terms.  It  should  be  remembered  that  the  total 
change  in  net  thrust  is 


F  -  F 


“  Ta*)  +  M*)  +  ^(Pa  "  ^a.) 

alg  h  aw  5  ^  apg  ®s  H 

+  ^  T  -  Tg  i  Mg  - 

2aTg2  at  2^m2  s  t 

2  2 

+  Pa  -  Pa  +  ?Ar.(Ta  -  T.  )(M-  -  M*) 

2^p2  H  H  aig^M  as  at  s  t 

;)2p 

4r(T,  -T,  )(P,  -  P,  )  +  ' 


algaPa  as  at 


n 


dM 


(66) 


Hence  only  the  remaining  second  order  terms  in  equation  (63) 
need  be  added  to  equation  (62)  co  find  standard  excess  thrust 
for  an  expansion  through  the  second  order  terms. 
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The  macinitudes  of  the  second  order  terms  .can  be  expected  to 
be  less  than  those  of  the  first  order  terms,  tut  they  v/ill  be 
of  consequence  in  some  instances.  It  is  suggested  that  for 
installations  having  an  accuracy  of  0.0 lo  or  less  in  n  consid- 
eration  be  given  to  inclucinn  seconu  order  terms  in  the  stan- 
cardization  equation.  They  can  be  safely  ignored  when  using 
the  airspeed-altitude  method  since  errors  introduced  by  this 
method  are  larger  than  the  maanitudc  of  th.e  second  order  terms. 
ITERATION  PROCEDURES 

Equation  (^2)  cannot  be  solved  directly  since  (during 

climbs)  standard  normal  load  factor  depends  on  y^,  and 

standard  weight,  which  must  be  determined  from  the  fuel  consumed. 

'f’he  fuel  consumed,  in  turn,  depends  on  i  \e  normal  load  factor. 

Procedures  are  simplified  for  level  accelerations  since  Yg=0. 

Further,  trim  drag  corrections  are  related  to  center  of  Gravity 

position  v;hich  may  be  a  function  of  weinht,  as  v/hen  automatic 

fuel  seauencina  is  used.  Lenco,  iteration  procedures  are 

reouired  to  find  standard,  values  of  normal  load  fcv-tor,  v/eight, 

ano  perhaps  center  of  Gravity  position.  It  should  3:e  noted  that 

several  test  points  at  least  r>;ust  he  operated  on  at  one  time  in 

order  that  nuruerical  differentiation  to  obtain  y  and  consequently 

n  liiav  bo  carrioc.  out. 
z 

For  purposes  of  illustration,  the  hasJe,  sto].s  in  the  iteration 
procedure  for  stt nc  arc'ization  of  a  continuous  climb  are  laid  out 
in  the  following  chart. 
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SECTION  VIII 

STi^NDARDIZATION  OF 
PERFORIEANCE  PARAMETfRS 
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SUMMARY 


Onca  standard  excess  thrust  has  been  found,  the  remaining 
calculations  to  compute  other  performance  parameters  of  interest 
are  relatively  simple.  Equations  to  determine  the  basic 
parameters  for  climbs  and  level  accelerations  have  been  included 
in  this  section. 
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SYMBOLS  USED  IN  THIS  SLCTIOW 


Symbol 

Definition 

Units 

acceleration  factor 

dimensionless 

^fE 

acceleration  factor  associated 
with  energy  height 

dimensionless 

ex 

excess  tlirust 

lb 

FU 

fuel  used 

lb 

9l 

local  effective  acceleration 
due  to  gravity 

ft/sec^ 

reference  acceleration  of  gravity 

ft/sec^ 

h 

geometric  altitude 

ft 

pressure  altitude 

ft 

energy  height 

ft 

load  factor  along  the  x-axis 

dimensionless 

NAMT 

nautical  air  miles  traveled  ' 

-  _  _ 

UGMT 

nautical  ground  miles  traveled 

_  _ 

r 

local  rodius  of  the  earth 

ft 

t 

time 

sec 

true  airspeed 

ft/sec 

w 

airplane  gross  weight 

IL 

Vv 

es 

airplane  gross  weiebt  at 
engine  start 

lb 

fuel  flOV7 

Ib/sec 

flightpath  cliri)  angle  measured 
from  the  geodetic  horizontal 
plane 

rad 

Sxibscript 

s 

standard  day  conditions 
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STANDARDIZATION  OP  PERFORMRaCE  PARAMETERS 

Once  excess  thru  jt  has  been  corrected  to  standard 
conditions  an  iteration  procedure  is  used  to  find  the  para¬ 
meters  of  concern  for  both  clirbs  and  level  accelerations. 
The  basic  steps  in  the  iteration  procedure  are  shown  for 
a  clirb  in  the  section.  Standardization  of  bxcess  Thrust. 

A  sinllar  prc>cedure  woulc  be  usea  for  i*..  vcl  accelerations 
althouqh  fewer  steps  are  required  since  v  .  o  and  excess 
.nrust  goes  into  increasing  kinetic  energy  only. 

CLIMB  PARAMETERS 

In  addition  to  the  need  for  iteration  proceciures  noteu 
above  for  both  clirbs  and  accelerations,  an  additional 
complication  arises  when  generatino  cliiib  paraiiveters .  It 
has  been  accepted  practice,  historically,  to  aeterndne  dim}; 
parameters  as  a  function  of  aitituce,  beginning  at  sea 
level.  An  aircraft  may  be  at  several  tnousand  feet, 
however,  before  a  pilot  is  able  to  intercept  the  aesirea 
climb  schedule,  hence,  it  is  not  possible  to  collect 
test  aata  throuah  a  range  of  altitudes  irmeaiately  aoove 
sea  level.  Since  cliiib  perfonuance  at  sea  level  is  useful 
for  making  comparisons  v/ith  other  aircraft  and  to  Stanuard 
Aircraft  Characteristics  Charts,  it  must  be  extrapolateu 
to  sea  level,  recognizing  that  the  flight  condition  is 
impossible  and  that  the  data  is  fictitious. 
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lo  obtain  cliinb  perforr.ance  at  sta  level  tiie  rate  of 
climb  is  extrapolated  to  sea  level;  other  paraiaeters  are 
then  deduced  from  rate  of  climb.  The  steps  necessary  for 
finding  standardized  rate  of  climb  might  be  carried  out 
entirely  in  a  digital  computer;  perhaps  using  a  curve  fit 
of  available  rate  of  climb  data  and  extending  it  to  sea 
level  with  the  aia  of  thrust  and  drag  data.  Due  to  un¬ 
certainties  in  measurea  rates  of  climb,  particularly  at 
the  lowest  test  altitudes,  the  extrapolation  iray  not 
produce  valid  results.  Hence,  it  may  be  preferable  to 
interrupt  the  computer  program  to  establish  rate  of  climb 
manually  and  then  restart  th«;  prog.:aiu,  continuing  v/ith 
the  sta?.dardization  process. 

Rate  of  Climb 

An  aircraft's  climb  potential  is  defineu  by 

p  V. 

*  ex^'t 

W 

and  represents  the  rate  of  climb  which  would  be  achieved  by 
ar.  i^i’.-plane  climbing  at  a  constant  true  speed.  If  a  climb 
is  made  at  other  than  a  constant  true  rnoed,  the  rate  of 
change  of  kinetic  energy  must  be  accounted  for  to  fine  rate 
of  climb.  This  is  done  through  the  acceleration  factor,  A^, 
and  the  rate  of  climb  is,  at  standard  conditions 


p  V. 

AjWg 


"x 

*s  ‘'s 


(1) 
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where 


kl 


(1  -t 


''t. 


8r  Cl, 


(2) 


'ine  acceleration  factor  is  deri\'ed  as  follows.  Rate 
of  change  of  energy  height  may  be  stated  as 


n 


Vt  dVt 


E 


=  177h^ 


(3) 


and  also  as 


H 


E  “ 


W 


(4) 


I’rom  equation  (10)  of  the  section.  Atmospheric  hnyiron/aent 


« 

h 


(5) 


Combining  equations  (3) ,  (4) ,  and  (5) 


r  '^K 


(6) 


from  which  the  acceleration  factor  is  taken  to  be  as  defined 
by  equation  (2) .  This  acceleration  factor  is  useful  for 
climbs  described  in  the  section,  Standara  Climb  Schedules, 
during  which  altitude  is  raonotonic  increasing.  For  the  case 
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where  altitude  is  permix.ted  to  decrease,  corrections  are  made 
St  cojAstant  values  of  II,,  and  the  acceleration  factor  takes  the 

b 

form 


A 


fE 


(7) 


is  found  by  defining  rate  of  change  of  energy  height  as 


Vt  dVj  . 


He  =  Hj.  +  ~  ^  Hg 


(8) 


Solving  for  H 


c 


Vt  dVt 

«c  = 


(9) 


Substituting  equation  (5) 


gr  Vj  dVj  . 

h  =  -  V  7ir>HE 

gL  gf  tSHg  E 


(10) 


or 


h  =  AfgHg 


(11) 
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Time  to  Climb 


Once  h  has  been  established  time  to  climb  follows 
s 


easily  from 


t  = 


ft 


(12) 


>1 


for  continuous  climbs  and 


t  = 


(13) 


for  optimum  climb  schedules  which  use  energy  height  as  an 
independ^rtv.  variable. 

Distance 

Distance,  conventionally  shown  in  terms  of  nautical  air 
miles  traveled,  is 


N/.V’T.  = 


V,  cos>'dt 


(14) 


The  NAI-IT  differs  s3.ightJ.v  from  the  distance  traveled  along 
the  earth's  surface,  particularly  for  aircraft  flying  at 
extreme  altitudes  (e.g,,  aoout  0.4%  at  an  altitude  of  80,000 
feet) .  The  ground  distance  is 


NGMT 


Vt-cos)’, 


S  '  S  r  +  h, 


(U 


(15) 
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Fuel  Used  and  Weight 


The  quantity  of  fuel  useii  is  computeu  simply  by  inte¬ 
grating  the  standardizes  fuel  flov, : 


Weight  during  a  clirii)  is  based  on  the  aircraft  gross  v;eight 
at  engine  start  and  the  fuel  allowance  for  taxi,  takeoff, 
and  acceleration  at  sea  level  to  the  cliiib  schedule.  Standard 
weight  during  tlie  course  of  a  cliirb  is,  then 


Wc  =  -  (fuel  allowance)  -  PU 

S  6S 


(17) 


LEVEL  ACCELERATION  PARAI-ETERS 

The  same  parai.ieters  (tire,  distance,  and  fuel  used)  are 
desired  in  the  presentation  of  ac  '-eleration  performance*  After 
standardized  excess  thrust  as  a  function  of  i'tnch  number  or 
true  speed  has  been  established  the  acceleration  can  be  found 
from 


and  the  tire  to  accelerate  from 


(18) 


(19) 
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The  ristarjce  traveled  is  computed  ?s  for  climbs.  With 
y  -=0  distance  over  the  sartij  is 

/** 

Fuel  used,  as  for  climbs,  is  determined  during  the  standardiz¬ 
ation  process.  Repeating  equation  (16) 
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